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APPLICATION  OF  NONDESTRUCTIVE  TESTING 
TECHNIQUES  TO  MATERIALS  TESTING 

G.  S.  Kino 

SUMMARY 

A  set  of  new  techniques  in  scanning  acoustic  and  optical  mi¬ 
croscopy  is  described.  Starting  with  an  acoustic  microscope  that  can 
directly  measure  both  phase  and  amplitude,  similar  techniques  have 
been  developed  for  the  scanning  optical  microscope.  These  make  it 
possible  to  measure  range  to  a  thousandth  of  a  wavelength.  Other  tech¬ 
niques  involving  scanning  optical  microscopy  have  been  demonstrated 
which  make  it  possible  to  carry  out  profiling  of  semiconductor  circuits. 
Developments  of  these  methods  are  now  being  actively  pursued  for  use 
in  the  semiconductor  industry.  Developments  of  the  acoustic  technol¬ 
ogy  pioneered  on  this  program  have  proven  extremely  important  for 
measuring  internal  defects  in  composites  and  surface  cracks  on 
ceramics. 
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INTRODUCTION 

The  purpose  of  this  program  has  been  to  develop  new  techniques  of  nondestructive 
testing  as  in  the  forerunning  AFOSR  contract.  Contract  No.  F49620-79-C-0217.  During 
the  earlier  contract  period  we  developed  acoustic  wave  techniques  to  examine  materials  of 
the  type  typically  used  in  aircraft  and  in  various  types  of  structural  components.  During  that 
time,  we  pioneered  the  development  of  new  types  of  electronically  scanned  acoustic  imaging 
systems  and  highly  successful  new  theoretical  approaches  for  determining  the  attenuation  of 
waves  in  polycrystalline  materials.  Toward  the  end  of  that  contract,  surrounded  as  we  are 
by  Silicon  Valley,  our  interests  began  to  change  toward  the  examination  of  small  devices 
such  as  integrated  circuits  and  other  electronic  components.  We  therefore  began  to  develop 
acoustic  microscopy  techniques  suitable  for  examining  composite  materials  and  structural 
ceramics,  as  well  as  electronic  components.  This  work  was  continued  on  the  contract  being 
reviewed  here. 

A  second  program  begun  on  the  first  contract  was  concerned  with  using  fiber-optic 
techniques  to  measure  the  displacement  of  the  surface  of  an  object  by  surface  acoustic  waves 
propagating  along  it  This  noncontacting  technique  was  used  to  measure  the  depth  and  size 
of  surface  cracks  and  other  defects  in  metals  and  ceramics,  and  proved  to  be  very  powerful 
for  measuring  the  properties  of  surface  acoustic  wave  devices.  It  is  now  in  regular  use  at 
Hewlett-Packard,  Santa  Rosa,  where  it  is  employed  for  measuring  wavefronts  in  their  sur¬ 
face  acoustic  wave  resonators. 

The  work  on  this  contract  led  to  major  new  developments  in  acoustic  and  optical 
microscopy.  One  example  is  our  development  of  the  first  acoustic  microscope  capable  of 
measuring  phase  directly.  With  this  microscope,  we  were  able  to  use  Fourier  transform 
techniques  to  directly  measure  the  reflectivity,  as  a  function  of  angle  of  plane  waves,  inci¬ 
dent  on  the  surface  of  a  material,  and  to  measure  the  thickness  of  metal  films  of  the  order  of 
1000  A  thick  using  waves  with  a  wavelength  of  30  pm  .W  The  devices  were  therefore 
very  sensitive  to  minor  changes  in  material  properties,  to  the  presence  of  small  defects,  and 
particularly  cracks.  As  an  example,  we  were  able  to  measure  the  stress  variation  across  a 
disc  of  glass  with  this  noncontacting  method.  The  work  on  acoustic  microscopy  is  now 
being  continued  by  Professor  Khuri-Yakub  on  another  program  where  the  techniques  de¬ 
veloped  out  of  the  AFOSR  program  have  proved  to  be  particularly  useful  for  measuring  the 
width,  depth,  and  features  of  defects  inside  epoxy  composite  materials,  and  for  the  mea¬ 
surement  of  cracks  in  ceramics  associated  with  surface  machining. 


The  work  on  fiber-optic  sensors  combined  with  our  background  in  acoustic  mi¬ 
croscopy,  which  gave  us  a  different  viewpoint  from  standard  optical  microscopists,  led  to  a 
range  of  new  types  of  scanning  optical  microscopes.  The  first  of  these,  an  electronically- 
scanned  scanning  optical  microscope,  was  developed  on  this  contract  and  was  the  first  opti¬ 
cal  microscope  capable  of  directly  measuring  phase  and  amplitude.^  One  example  of  the 
use  of  this  technique  is  that  of  a  precision  range  sensor  suitable  for  use  in  robotics.^ 

Later  developments  arising  from  this  work  have  included  yet  other  types  of  real-time 
scanning  optical  microscope,  and  a  device  called  the  correlation  microscope.  The  first  of 
these  devices  is  now  being  developed  commercially  from  our  design.  The  major  advantage 
of  such  scanning  optical  microscopes  is  that  they  have  excellent  depth  resolution  along  with 
transverse  resolution,  and  contrast  is  somewhat  better  than  with  standard  microscopes.  This 
makes  it  possible  to  measure  the  profiles  of  integrated  circuits  with  low  submicron 
resolutions,  and  to  determine  the  thickness  of  films  with  resolutions  of  the  order  of  10  A  . 

Three  theses  were  written  during  the  time  of  this  contract.  One  was  on  precision 
phase  measurements  in  acoustic  microscopy.  5  The  second  was  on  the  fiber-optic  sensor 
and  the  electronically  scanned  scanning  optical  microscope.^  Both  of  these  devices  resulted 
in  patents.  The  third  thesis  concerned  itself  with  the  development  of  the  electronically- 
scanned  scanning  optical  microscope.^ 

We  shall  discuss  this  work  in  more  detail  below,  and  refer  the  reader  to  the  Ap¬ 
pendix  and  our  published  papers  for  still  fuller  descriptions  of  the  research. 


DESCRIPTION  OF  PROJECTS 
Phase  Measuring  Acoustic  Microscope 


A  50  MHz  phase  measuring  acoustic  microscope  was  constructed  as  shown  in 
Fig  1.8  A  ceramic  transducer,  with  a  20%  bandwidth,  was  indium  bonded  to  a  quartz 
buffer  rod  in  which  was  ground  a  a  3.2  mm  spherical  depression  which  formed  a  lens 
with  a  focal  length  of  4.24  mm  in  water.  When  the  surface  of  a  solid  was  placed  nearer  to 
the  lens  than  the  focal  plane,  two  waves  could  be  set  up  by  a  short  tone  burst  exciting  the 
transducer.  The  first  corresponded  to  direct  specular  reflection  from  the  surface,  and  the 
second  to  the  excitation  of  a  Rayleigh  wave  which  propagated  along  the  surface.  The  re¬ 
ceived  tone  burst  corresponding  to  the  Rayleigh  wave  is  delayed  relative  to  the  specular  re¬ 
flection.  Thus,  it  was  possible  to  gate  the  two  tone  bursts  and  switch  them  into  separate 
circuits.  By  this  means,  we  were  able  to  design  a  circuit  which  could  measure  the  phase 
difference  of  the  carriers  corresponding  to  the  two  tone  bursts,  and  hence  the  surface  wave 
velocity  along  the  substrate.  By  using,  instead,  a  reference  corresponding  to  a  pulse  re¬ 
flected  from  the  flat  sides  of  the  buffer  rod,  it  was  also  possible  to  measure  the  relative 
phase  or  time  delay  of  the  specularly-reflected  pulse  from  the  object,  and  hence  the  distance 
from  the  object  with  good  accuracy. 
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Fig.  1.  Schematic  of  the  phase  measurement  system. 


We  developed  the  phase  measurement  system  described  in  more  detail  in  reference 
8,  which  is  part  of  the  Appendix.  Suffice  it  to  say  that  with  pulses  only  two  or  three  cycles 


long,  it  is  possible  to  measure  the  phase  to  an  accuracy  of  the  order  O.lO  ;  hence  we  were 
able  to  measure  changes  in  velocity  along  the  Rayleigh  wave  path,  which  is  of  the  order  of 
1-2  mm  long,  to  accuracies  of  one  part  in  10^. 

In  the  accompanying  papers  in  the  Appendix,  1  >2,8  examples  are  given  of  the  power 
of  this  device.  We  were  able  to  measure  the  thickness  of  thin  indium  overlays  of  the  order 
of  1000  A  thick,  with  good  accuracy  and  also  the  stress  distribution  in  a  glass  plate.  By 
using  Fourier  transforms  with  the  phase  and  amplitude  information  obtained,  we  were  able 
to  measure  the  reflectance  function  as  a  function  of  angle  for  a  plane  wave  incident  on  vari¬ 
ous  types  of  substrates. 

These  techniques  required  considerable  development  of  fundamental  new  phase 
measurement  methods  for  short  tone  bursts.  At  the  same  time,  we  developed  a  considerable 
body  of  theory  to  interpret  the  measurements.  This  theory  has  been  fundamental  to  our  un¬ 
derstanding  of  both  scanning  acoustic,  and  the  later  scanning  optical,  microscopes. 

Electronically-Scanned  Scanning  Optical  Microscope 

This  device  has  the  capability  of  measuring  both  optical  amplitude  and  phase.  The 
ideas  behind  it  are  basic  to  much  of  the  present  work  we  are  carrying  out  on  new  programs 
on  phase  contrast  microscopy,  and  the  work  developed  towards  the  end  of  the  AFOSR 
program  on  mechanically-scanned  scanning  optical  microscopes.3>6,7 ,9,io 

In  the  electronically-scanned  microscope,  shown  in  Fig.  2,  a  collimated  laser  beam 
passes  through  a  Bragg  cell  and  is  diffracted  into  two  beams,  one  of  which  may  be  scanned 
by  varying  the  input  frequency  fb  to  the  Bragg  cell.  The  other  stationary  undiffracted  beam 
is  used  as  a  reference.  Both  beams  are  passed  through  transfer  lenses  into  the  microscope 
objective  and  are  focused  to  two  spots  on  the  object.  The  light  from  the  two  spots  passes 
back  through  the  lenses  and  the  Bragg  cell.  The  reference  beam  is  diffracted  into  a  beam 
which  travels  in  the  same  direction  as  the  reflected,  scanned  beam.  The  two  coincident 
beams,  one  of  which  is  upshifted  by  a  frequency  fb  and  the  other,  which  is  downshifted 
by  fb ,  pass  into  a  detector  and  an  output  signal  is  obtained  at  a  frequency  2fb . 
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Fig.  2.  Schematic  of  the  electronically-scanned  microscope. 

As  the  beam  is  scanned  across  the  object,  the  output  signal  from  the  detector  at  a 
frequency  2fb  in  the  100-200  MHz  range  has  an  amplitude  which  varies  with  the  reflec¬ 
tivity  of  the  sample,  and  a  phase  which  depends  on  the  phase  difference  between  the  refer¬ 
ence  beam  and  the  moving  beam.  Since  phase  can  be  measured  to  an  accuracy  of  0.1°  , 
this  implies  that  variations  in  height  can  be  measured  to  an  accuracy  of  the  order  of  1  A  . 
In  our  experimental  system,  phase  and  amplitude  can  be  measured  to  this  accuracy  at  a  rate 
of  20  jxs/spot .  As  far  as  we  are  aware,  this  system  is  one  of  the  fastest  for  measuring 
phase  and  amplitude  simultaneously  that  has  been  developed.  The  results  obtained  with  it, 
and  the  method  of  digitization  of  the  data,  are  described  in  the  accompanying  articles  in  the 
Appendix.^’ll 

The  advantages  of  phase  contrast  microscopy  are  already  familiar.  In  reflection  mi¬ 
croscopy,  measurement  of  phase  makes  it  possible  to  measure  profiles,  or  height,  indepen¬ 
dently  of  changes  in  reflectivity.  In  addition,  contrast  between  regions  of  different  heights, 
but  with  the  same  reflectivity,  such  as  a  groove  in  metal,  are  greatly  enhanced  by  the  use  of 
a  phase  image.  As  an  example,  a  line  scan  of  a  titanium-gold  line  on  silicon  is  shown  in 
Fig.  3.  It  will  be  seen  that  phase  measurements  give  better  results  than  amplitude  mea¬ 
surements  in  this  case.  Using  a  0.95  N.A.  objective  lens,  with  argon  laser  illumination  of 
a  wavelength  of  510  nm  ,  the  edge  definition  of  the  phase  scan  of  a  step  of  aluminum  on 
aluminum,  from  the  10-90%  height  points,  is  250  nm  .  Furthermore,  because  both 
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amplitude  and  phase  information  is  available,  a  Fourier  transform  of  the  line  scan  can  be 
taken,  and  inverse  filtering  processes  carried  out.  In  one  example,  the  10-90%  definition 
was  shown  to  be  improved  to  130  nm ,  with  similar  results  for  gratings  and  for  finite 
width  strips. 

The  Scanning  Optical  Microscope 

We  have  mainly  been  concentrating  on  optical  methods  in  the  last  few  years,  because 
the  developments  in  our  own  and  other  laboratories  of  the  Confocal  Scanning  Optical 
Microscope  (CSOM),  in  addition  to  providing  noncontacting,  nondestructive  methods  of 
measurement,  make  it  possible  to  measure  depth  and  width  at  the  same  time  with  submicron 
resolutions.  The  transverse  resolution  is  better  than  that  of  a  standard  optical  microscope 
using  the  same  lens.  Optical  techniques  do  not  damage  the  material,  they  offer  great 
flexibility,  and  are  easy  to  use.  Unlike  the  electron  microscope,  the  optical  techniques  do 
not  require  a  vacuum  and  it  is  easy  to  carry  out  sectioning,  nondestructively,  to  obtain  in¬ 
formation  on  features  at  different  depths  in  the  interior  of  a  translucent  material  or  to  find  the 
profile  of  surface  features.  Confocal  scanning  microscopy,  because  of  its  excellent  range 
and  transverse  definitions,  can  provide  a  facility  for  optical  "cross-sectioning"  that  allows 
glare-free  fluorescence  imaging  of  fine  biological  structures  and  imaging  of  internal  features 
in  translucent  materials. 

New  tools,  with  submicron  resolution,  are  much  needed  for  nondestructive  measurements 
of  small  structures  in  semiconductors,  integrated  circuits,  magnetic  recording  heads,  and 
micromachined  samples.  The  same  tools  are  of  great  importance  to  manufacturing  tech¬ 
nology  because  they  will  make  it  possible  to  carry  out  precision  metrology  and  detect  faults 
in  parts  during  the  manufacturing  process.  Since  they  employ  noncontacting,  nondestructive 
methods,  they  can  perform  a  large  number  of  measurements  in  a  short  time.  The  basic 
principle  of  the  scanning  optical  microscope  is  shown  in  Fig.  4.  Light  passes  through  a 
pinhole  to  a  microscope  objective  lens  and  the  pinhole  is  imaged  on  the  object.  The  light 
reflected  from  the  object  passes  through  the  objective  lens  and  back  through  the  pinhole.  If 
the  object  is  moved  out  of  focus,  the  reflected  light  reaching  the  pinhole  is  defocused  and  the 
light  does  not  pass  through  it  Consequently,  if  a  detector  is  placed  behind  the  pinhole  and 
the  amplitude  of  the  signal  received  varies  rapidly  with  the  defocus  distance  z ,  the  distance 
dz(3  dB)  between  3  dB  points  can  be  shown  to  be  given  by  the  approximate  formula:  ^ 


.  _  0.45*. 

2  1  —  cos  0  Q 

where  sin  00  =  N.A.  is  the  numerical  aperture  of  the  lens. 

Typically,  with  a  wide  aperture  lens,  the  3  dB  points  of  this  response  are  only 
0.5  fim  apart.  At  the  same  time,  the  transverse  response  of  this  microscope  tends  to  be 
better  than  that  of  a  standard  optical  microscope  because  the  lens  is  used  both  on  transmit 
and  on  receive.  Consequently,  the  point  spread  function  is  squared,  the  spatial  frequency 
response  is  doubled,  and  the  3  dB  width  of  the  transverse  response  is  decreased  by  a 
factor  of  0.7  from  that  of  a  standard  microscope  using  the  same  lens.  Either  a  laser  source 
or  an  incoherent  source  can  be  used,  since  the  sidelobe  level  of  the  response  is  so  low  that 
speckle  caused  by  interference  from  regions  outside  the  focused  spot  is  no  longer  a  prob¬ 
lem. 

We  have  constructed  a  scanning  optical  microscope  based  on  this  principle.  We  mechani¬ 
cally  scan  the  object  and  construct  a  raster  image  which  is  stored  in  a  computer  and  then 
displayed  as  a  video  image.  It  takes  a  few  seconds  to  form  one  frame. 1 1  We  have  also 
constructed  an  optical  range  finder  using  the  same  technique.4  Results  taken  with  both  de¬ 
vices  are  described  in  the  Appendix. 

Towards  the  end  of  the  program  we  devised  new  techniques  for  scanning  differential 
contrast  microscopy.10  These  techniques  are  described  in  the  Appendix.  The  basic  purpose 
behind  these  techniques  is  to  get  a  better  definition  of  edges  and  of  heights  by  using  differ¬ 
ential  methods.  When  the  object  is  moved  in  a  direction  z ,  the  output  in  amplitude  passes 
through  a  maximum.  Obviously,  it  would  be  easier  to  find  where  this  maximum  is  and  on 
which  side  of  the  maximum  the  focus  is  by  differentiating  the  curve.  A  very  simple  way  of 
doing  this  is  to  move  the  object  back  and  forth  periodically  at  a  frequency  Q  .  The  signal 
output  at  this  frequency  Q  will  essentially  be  die  differential  of  the  intensity  as  a  function 
of  distance.  A  better  way  to  do  this,  which  does  not  require  mechanical  movement,  is  to  use 
an  electro-optic  cell  made  of  PLZT  and  modulate  the  phase  of  the  wave  passing  through  the 
cell  using  transparent  electrodes.  This  technique  is  somewhat  equivalent  to  a  Zemike 
microscope  but  has  better  contrast  and  better  definition.  Results  obtained  with  this  method 
are  described  in  reference  10  in  the  Appendix. 
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Fig.  4.  An  illustration  of  the  principles  of  confocal  microscopy. 

A  second  technique  involves  using  a  Bragg  cell  with  two  spots  placed  side  by  side 
and  moved  with  each  other.  This  can  give  a  differential  output  with  respect  to  transverse 
position.  Such  a  technique  is  very  useful  for  measuring  slope  or  the  position  or  sharp 
edges. 

Since  the  AFOSR  program  was  on  a  cost-free  extension  from  December  1986  to 
September  1987,  this  work  was  placed  on  another  program  in  December  1986. 

CONCLUSIONS 

The  work  described  here  describes  major  advances  in  both  acoustic  and  optical  mi¬ 
croscopy.  New  techniques  for  directly  measuring  phase  and  its  range  have  been  demon¬ 
strated.  The  use  of  these  techniques  gives  excellent  transverse  definition,  along  with  good 
range  definition.  The  work  is  continuing  on  other  programs  and  the  new  developments  in 
optical  microscopy  have  excited  a  great  deal  of  interest  in  the  semiconductor  industry. 
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Abstract— The  measurement  and  the  use  of  phase  in  acoustic  micros¬ 
copy  are  discussed.  It  is  demonstrated  that  in  many  applications  phase 
can  be  used  to  provide  sensitivity  and  information  unparalleled  by  am¬ 
plitude-only  measurement  methods.  A  technique  capable  of  high- 
accuracy  measurement  of  the  phase  of  short  RF  acoustic  pulses  is  de¬ 
scribed.  The  power  of  this  phase  measurement  technique  is  illustrated 
in  a  number  of  applications.  Surface  material  property  measurements 
such  as  the  Rayleigh-wave  velocity  and  the  inversion  of  the  complex  F(;) 
to  obtain  the  reflectance  function  of  a  liquid-solid  interface  are  consid¬ 
ered.  Surface  topography  mapping  based  on  phase  measurement  is  ex¬ 
amined.  A  Fourier  transform  approach  for  precision  determination  of 
linewidlhs  comparable  to  the  resolution  spot  size  is  also  presented. 

I.  Introduction 

HE  SCANNING  acoustic  microscope  is  a  high-reso¬ 
lution  imaging  system  in  which,  unlike  conventional 
optical  microscopes,  it  is  relatively  straightforward  to 
measure  the  phase  of  the  return  signal.  With  only  a  few 
exceptions  [l]-[3],  attention  has  centered  only  on  the  use 
of  intensity  information.  In  some  of  the  main  areas  of  ap¬ 
plication  of  the  acoustic  microscope,  the  phase  of  the  re¬ 
ceived  signal  plays  an  important  role.  For  example,  in  the 
so-called  F(z)  measurements  [4]-[7],  the  amplitude  of  the 
received  signal  V(z)  as  a  function  of  the  separation  be¬ 
tween  the  lens  and  the  substrate  exhibits  periodic  peaks 
and  nulls.  This  phenomenon  is  due  to  the  beating  between 
a  specularly  reflected  signal  from  the  substrate  and  a  de¬ 
layed  leaky  Rayleigh-wave  signal,  which  reemits  to  the  lens 
while  propagating  along  the  surface  of  the  substrate.  The 
phase  difference  between  these  two  signal  components  is 
therefore  of  great  importance,  and  in  fact  it  controls  the 
contrast  of  reflectionai  images.  In  addition,  it  has  been 
shown  that  the  independent  measurement  of  phase  and 
amplitude  can  be  very  useful  in  the  determination  of  the 
elastic  constants  of  tissue  [8],  [9],  However,  little  effort 
has  been  made  to  extract  the  phase  of  the  return  signal 
separately. 

One  reason  for  the  reluctance  to  make  use  of  the  phase 
is  that  it  is  not  generally  trivial  to  measure  the  phase  of  a 
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high-frequency  tone  burst  with  sufficient  accuracy.  We  will 
describe  here  an  approach  that  is  relatively  easy  to  imple¬ 
ment  and  capable  of  yielding  high-precision  phase  data, 
even  from  very  short  tone  bursts.  The  details  of  the  phase 
measurement  technique  will  be  published  elsewhere,  and 
we  will  concentrate  here  on  some  of  the  applications  of 
the  system  to  acoustic  microscopy. 

There  are  many  motivations  for  making  measurements 
of  this  type.  Firstly,  they  offer  a  direct  indication  of  ma¬ 
terial  properties  as  in  the  measurement  of  Rayleigh-wave 
velocity.  Also,  the  results  of  the  phase  measurement  can 
be  used  to  infer  the  width  and  the  height  of  surface  fea¬ 
tures.  Finally,  the  combination  of  amplitude  and  phase 
measurements  can  be  used  in  the  reconstruction  processes 
in  which  complete  information  about  the  interaction  be¬ 
tween  acoustic  field  and  material  is  required,  as  in  the  in¬ 
version  of  V(z)  data  to  find  the  reflection  coefficient  as  a 
function  of  angle  [10]. 

II.  Phase  Measurement  Scheme 
A.  Acoustic  Lens  Configuration 

In  making  precise  phase  measurements  that  are  related 
to  physical  properties,  it  is  generally  important  to  ensure 
that  the  reference  signal,  against  which  the  phase  of  the 
probing  signal  is  to  be  compared,  and  the  probe  itself  share 
as  many  of  the  instrumental  and  environmental  phase  dis¬ 
turbances  as  possible.  In  other  words,  precise  measure¬ 
ments  are  best  done  in  an  interferometer,  where  the  two 
arms  are  closely  matched. 

In  this  work  we  have  made  use  of  two  acoustic  mea¬ 
surement  configurations  that  largely  satisfy  this  general 
condition.  The  first  (Fig.  1(a))  is  reminiscent  of  the  de- 
focused  condition  used  in  V(z)  measurements,  except  that 
here  a  fixed  separation  between  the  lens  and  the  sample 
surface  is  maintained.  Two  components  of  the  acoustic 
field  returned  to  the  lens  contribute  most  significantly  to 
the  output  signal:  the  on-axis  specular  reflection  of  the 
longitudinal  wave  in  the  water  (L  in  Fig.  la);  and  the  off- 
axis  rays  (/?),  which  satisfy  the  condition  vw/vr  =  sin  $R 
for  the  conversion  of  longitudinal  waves  with  velocity  vw 
in  the  water  to  leaky  Rayleigh  waves  with  velocity  vR  on 
the  surface  of  the  sample  [4],  [5], 

These  two  signals,  L  and  R,  differ  slightly  in  path  length 
in  water  but  otherwise  experience  much  the  same  environ- 
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Fig.  I.  Acoustic  lens  configurations,  (a)  Rayleigh-wave  velocity  perturba¬ 
tion  measurement,  (b)  Topography  mapping. 


mental  disturbances.  However,  the  Rayleigh-wave  com¬ 
ponent  acquires  substantial  phase  delay  at  the  surface  of 
the  specimen,  hence  the  effective  path  length  for  the  R 
signal  is  longer  than  that  of  the  specular  reflection.  Indeed 
in  our  experiments  we  use  the  different  arrival  time  of 
these  two  return  signals  to  facilitate  their  separation  and 
subsequent  phase  comparison.  The  excitation  signal  is  two 
to  three  cycles  in  duration,  and  the  defocus  distance  is 
sufficient  that  there  is  no  temporal  overlap  of  the  signals. 
It  is  clear  that  changes  in  the  Rayleigh-wave  velocity  as  a 
function  of  position  along  the  surface  of  the  specimen  may 
be  sensed  in  this  way.  By  combining  a  high-accuracy  phase 
measurement  scheme  with  a  suitably  extended  Rayleigh- 
wave  path,  remarkably  sensitive  measurements  can  be 
made. 

The  second  configurauon  that  concerns  us  is  illustrated 
in  Fig.  1(b).  In  this  case  the  lens  is  positioned  so  that  its 
focus  is  at  the  surface  of  the  specimen  or  slightly  above  it 
so  that  no  Rayleigh  wave  of  importance  is  excited.  A  ref¬ 
erence  path  is  provided  by  an  annular  beam  that  propa¬ 
gates  through  the  flat  outer  periphery  of  the  lens.  By  ex¬ 
citing  the  lens  with  a  short  pulse,  and  once  again  using 
time  discrimination  to  separate  the  signals  from  the  two 
different  paths,  an  interferometer  is  formed.  Now  small 
local  changes  in  the  surface  topography  result  in  large 
changes  in  the  phase  of  the  focused  beam  relative  to  the 
phase  of  the  large  diameter  reference  beam. 

In  both  measurement  configurations  it  is  apparent  that 
the  signals  of  interest  arrive  at  different  times,  and  this 
would  ordinarily  make  phase  comparison  impossible. 
However,  the  electronic  system  described  in  the  next  sec¬ 
tion  essentially  reconstructs  two  phase-coherent  continu¬ 
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ous  signals  that  are  directly  related  to  the  amplitudes  and 
phases  of  the  two  time-distinct  return  signals. 

The  first  configuration  is  also  used  where  the  signal  of 
interest  is  the  transducer  output  voltage,  when  the  system 
is  operated  in  the  F(z)  mode.  In  this  case  the  exciting 
signal  has  a  much  longer  duration  so  that  the  return  signals 
(L  and  R)  overlap  and  interfere  at  the  transducer  even  for 
large  defocus  distances.  The  idea  here  is  to  measure  both 
the  amplitude  and  the  phase  of  the  signal  relative  to  some 
reference  as  a  function  of  the  axial  separation  z  between 
the  acoustic  lens  and  the  specimen.  The  resulting  complex 
V(z)  can  be  inverted  [10]  to  obtain  the  angular  depen¬ 
dence  of  the  reflectance  of  the  liquid-specimen  interface. 
From  the  reflectance  function,  various  material  property 
parameters  can  readily  be  extracted. 

In  the  experiments  described  here,  a  center  frequency 
of  50  MHz  was  used  with  an  acoustic  transducer  having  a 
bandwidth  of  20  percent.  The  lens  material  was  fused 
quartz  and  had  a  radius  of  curvature  of  3.2  mm,  giving  a 
focal  length  in  water  of  4.24  mm.  The  opening  aperture 
of  the  lens  was  5.0  mm  in  diameter,  corresponding  to  a 
maximum  half  angle  of  36°  or  an  f-number  of  0.85. 

Providing  the  same  fractional  bandwidth  can  be  main¬ 
tained,  there  should  be  little  difficulty  in  applying  the  same 
techniques  to  a  system  operating  at  much  higher  fre¬ 
quency,  perhaps  up  to  2  GHz.  Beyond  that  frequency,  lim¬ 
itations  in  the  electronic-switching  components  presently 
available  may  present  difficulties. 

B.  Electronics 

The  essential  elements  of  the  phase  measurement  elec¬ 
tronics  are  shown  schematically  in  Fig.  2.  The  system  will 
be  described  in  detail  in  a  forthcoming  paper  [18],  where 
various  sources  of  error  and  an  analysis  of  the  ultimate 
performance  are  discussed.  A  more  general  review  of  the 
concepts  will  suffice  here. 
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The  basic  scheme  involves  signal  recovery  by  synchro¬ 
nous  detection  and  subsequent  phase  measurement  with  a 
lock-in  amplifier.  The  key  component  (Fig.  2)  in  the  signal 
source  is  the  single  sideband  generator  (SSB)  with  syn¬ 
chronous  outputs  at  100  kHz,  10.6  MHz,  and  10.7  MHz. 
The  10.7-MHz  output  is  in  fact  the  upper  sideband  of  the 
product  of  the  100  kHz  and  10.6-MHz  signals.  The  lower 
sideband  at  10.5  MHz  is  suppressed  by  at  least  50  dB  using 
a  standard  FM  radio  system  IF  filter.  The  translation  os¬ 
cillator  shifts  the  operating  frequency  up  to  the  desired 
center  frequency  of  the  acoustic  system,  which  is  50  MHz 
in  this  case.  A  reference  signal  at  49.9  MHz  is  also  gen¬ 
erated.  Again  the  lower  sideband  components  of  the  prod¬ 
uct  signals  have  to  be  removed.  At  this  point  however  the 
filtering  requirements  are  not  as  stringent  as  in  the  SSB 
generator  and  tunable  bandpass  filters  (F3  and  F4)  with 
five-percent  bandwidth  are  adequate,  since  the  sidebands 
are  now  widely  separated.  This  signal  generation  scheme 
is  flexible  in  that  the  operating  frequency  is  tunable  over  a 
fairly  wide  range,  which  is  limited  in  this  case  to  100  MHz 
by  the  bandpass  filters.  For  much  higher  frequency  of  op¬ 
eration,  filters  with  greater  selectivity  would  be  needed  or, 
alternatively,  successive  stages  of  heterodyning  could  be 
used  with  filtering  at  each  stage  to  step  up  incrementally 
to  the  desired  operating  frequency. 

As  illustrated  in  Fig.  2  with  a  50-MHz  acoustic  micro¬ 
scope,  the  50-MHz  continuous  wave  (CW)  signal  is  time¬ 
gated  to  produce  a  short  tone  burst  that  excites  the  acous¬ 
tic  transducer.  Switch  SW3  is  used  as  a  time  gate  to  pass 
the  low-level  acoustic  return  signals  of  interest  and  to  block 
high-level  extraneous  ones,  which  may  damage  the  pream¬ 
plifier.  The  acoustic  return  signals  are  time-separated  RF 
pulses,  and  they  are  electronically  separated  through  time¬ 
gating  into  two  channels  as  shown.  Each  of  the  resulting 
signals  is  mixed  with  the  49.9  MHz  reference  signal,  and 
the  product  is  narrow-band  filtered  to  extract  the  100-kHz 
component. 

It  can  be  shown  [18]  that  the  100-kHz  signals  are  es¬ 
sentially  low-frequency  CW  replicas  of  the  pulse  modu¬ 
lated  RF  acoustic  signals,  bearing  identical  phase  and  am¬ 
plitude  information.  With  a  biphase  lock-in  amplifier  tuned 
to  100  KHz,  one  can  readily  measure  the  phase  difference 
between  the  two  signals  and  also  the  amplitude  of  the  sig¬ 
nal  being  fed  to  the  “signal”  channel  of  the  lock-in,  if  it 
is  required.  Phase  sensitivity  is  basically  limited  by  the 
lock-in  amplifier  because  the  problem  of  system  noise  can 
in  most  cases  be  overcome  by  increasing  the  integration 
time  in  the  lock-in  output  stage.  This  effectively  reduces 
the  system  noise  bandwidth,  but  it  also  includes  the  ob¬ 
vious  disadvantage  of  longer  data  acquisition  time.  Good- 
quality  lock-in  amplifiers  routinely  have  phase  resolution 
of  the  order  of  0.1°,  which  corresponds  to  an  overall  sys¬ 
tem  phase  sensitivity  of  1/3600  of  a  wavelength. 

In  velocity  perturbation  measurements  the  measured 
phase  is  dependent  on  the  distance  between  the  lens  and 
the  specimen.  This  distance  changes  with  the  sample  sur¬ 
face  topography  and  is  susceptible  to  thermal  drift  as  well. 


thus  introducing  phase  error  into  the  measurement.  To 
minimize  this  error,  a  feedback  mechanism  is  employed 
to  keep  the  lens-to-specimen  spacing  constant.  The  acous¬ 
tic  lens  is  mounted  on  a  piezoelectric  (PZT)  stack  so  that 
its  vertical  position  can  be  adjusted  continuously  by  an  '■ 

electronic  control  signal.  The  acoustic  on-axis  reflection  * 

pulse,  whose  phase  is  a  direct  measure  of  the  lens  to  sam-  * 

pie  distance,  is  applied  to  the  reference  channel  of  the  lock-  * 

in  amplifier  as  shown  in  Fig.  2.  The  reference  channel  f 

generates  a  100-kHz  constant-amplitude  phase-locked  ] 

replica  of  the  reference  input,  which  is  then  compared  with  1 

the  100-kHz  output  of  the  SSB  generator  to  produce  the  * 

control  signal  for  the  PZT  stack.  The  acoustic  lens  auto-  ! 

maticaily  tracks  the  surface  topography  of  the  sample  dur-  } 

ing  scanning  to  ensure  that  the  measured  phase  change  is  ] 

due  to  material  property  variation  alone.  This  feedback  j 

mechanism  is  also  used  effectively  in  such  modes  of  op-  j 

eration  as  topography  measurement  to  compensate  for  \ 

thermal  effects. 

III.  Applications 

The  power  of  this  measurement  scheme  is  illustrated  in 
the  following  with  a  number  of  examples.  We  will  consider 
both  material  property  measurements,  such  as  Rayleigh- 
wave  velocity  and  complex  reflectivity,  and  surface  topog¬ 
raphy  measurements. 

A.  Velocity  Perturbation  Measurements 

In  this  section  we  will  specifically  deal  with  the  pertur¬ 
bation  to  the  Rayleigh-wave  propagation  velocity  as  a  re¬ 
sult  of  material  property  change  and  also  the  presence  of 
surface  residual  stress. 

1)  Velocity  Perturbation  Due  to  Thin-Film  Over¬ 
lay:  The  sample  in  this  example  is  a  multiple-thickness 
indium  film  deposited  on  glass.  The  thicknesses  are  240 
A  and  620  A ,  respectively.  The  objective  was  to  measure 
the  perturbation  of  the  Rayleigh-wave  velocity  caused  by 
the  indium  film.  A  line  scan  over  the  surface  of  the  sample 
(Fig.  3)  exhibits  phase  changes  of  7°  and  11°  for  the  240- 
A  and  380- A  step  changes  in  film  thickness.  The  spatial 
resolution  of  the  system  is  defined  by  the  Rayleigh-wave 
path  length  on  the  substrate  and  is  determined  to  be  about 
0.8  mm  from  the  step  transition  width  in  the  line  scan.  It 
can  be  calculated  from  first  order  perturbation  theory  [11] 
that  the  velocity  perturbation  due  to  the  indium  film  is  0.18 
percent  and  0.46  percent  for  the  240- A  and  620- A  layers, 
respectively.  Based  on  these  estimated  parameters,  one 
would  expect  phase  changes  of  7.6°  and  12°  for  the  240- 
A  and  380- A  step  transitions.  Hence  there  is  fairly  good 
agreement  between  the  experimentally  obtained  and  the¬ 
oretically  predicted  phase  changes.  It  should  also  be  noted 
that  in  the  line  scan,  the  small  phase  variations  in  the  sup¬ 
posedly  flat  regions  of  the  indium  film  are  real  and  re¬ 
peatable.  The  fluctuations  are  less  than  0.5°  and  are  due 
to  nonuniformity  in  the  thickness  of  the  indium  film.  Since 
the  phase  sensitivity  of  the  system  is  limited  by  the  lock- 
in  amplifier  to  0.1°,  this  measurement  technique  can  po- 
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Fig.  3.  Measurement  of  velocity  perturbation  due  to  indium  film  on  glass. 

tentially  detect  velocity  perturbations  on  the  order  of  one 
part  in  10s. 

2)  Residual  Stress  Measurement:  In  nondestructive 
testing,  knowledge  of  the  surface  residual  stress  distribu¬ 
tion  in  a  component  plays  an  important  role  in  the  predic¬ 
tion  of  failure  modes.  One  way  of  characterizing  residual, 
stress  is  by  measuring  the  acoustic  wave  propagation  ve¬ 
locity,  which  varies  linearly  with  the  local  stress  [12].  With 
the  measurement  configuration  shown  in  Fig.  1(a)  the 
measured  change  in  the  relative  phase  between  the  L  and 
R  signal  pulses  as  the  lens  is  scanned  can  be  shown  to  be 
directly  proportional  to  the  residual  stress  on  the  object 
surface  [13]. 

The  sample  used  in  the  experiment  is  a  glass  disk  of  two 
inches  in  diameter  that  was  heated  and  then  thermally 
quenched  with  air  jets  to  introduce  a  radial  distribution  of 
residual  stress.  Fig.  4(b)  shows  the  radial  variation  of  the 
measured  phase  perturbation.  For  comparison  a  destruc¬ 
tive  test  was  carried  out  on  a  similar  sample  to  determine 
the  actual  residual  stress  distribution.  A  Vicker’s  indenter 
was  used  at  a  prescribed  load  to  produce  median  cracks  at 
different  points  on  the  glass  disk.  By  measuring  the  crack 
lengths,  the  residual  stress  as  a  function  of  radial  distance 
was  calculated  [14],  The  result  of  the  destructive  test  is 
shown  in  Fig.  4(a).  The  variation  of  the  principal  stresses 
are  essentially  the  same,  with  the  radial  and  tangential 
components  tracking  each  other  to  within  20  MPa.  Since 
the  phase  measurement  is  omnidirectional  in  that  a  spher¬ 
ical  lens  was  used  to  launch  Rayleigh  waves  propagating 
in  all  directions  along  the  surface,  the  resulting  phase  per¬ 
turbation  is  a  measure  of  the  sum  of  the  principal  stresses. 
Comparison  of  Figs.  4(a)  and  4(b)  shows  that  there  is 
fairly  good  corroboration  between  the  phase  perturbation 
curve  and  the  actual  residual  stress  distribution.  From 
these  results  we  empirically  deduce  that  40  MPa  of  stress, 
which  is  the  change  from  the  center  of  the  sample  to  the 


-25  0  ‘25 

RADIAL  DISTANCE  (mm) 

Fig.  4.  Residual  stress  measurements  on  the  surface  of  a  tempered  glass 
disk,  (a)  Variation  of  radial  (open  symbols)  and  tangential  stresses  (closed 
symbols)  across  a  diameter  of  the  glass  disk  estimated  from  indentation 
fracture  tests,  (b)  Radial  variation  of  measured  phase  perturbation. 

limit  of  the  line  scan,  corresponds  to  a  phase  variation  of 
30°. 

B.  Surface  Topography  Measurements 
We  have  used  the  acoustic  microscope  in  the  measure¬ 
ment  configuration  shown  in  Fig.  1(b)  as  a  high-resolution 
noncontacting  profilometer  [13],  [15].  Besides  being  able 
to  make  use  of  phase  profile  to  map  the  depth  variation, 
we  can  also  determine  with  high  precision  the  transverse 
profile  of  surface  features  in  special  cases.  Of  particular 
interest  in  this  regard  is  the  measurement  of  the  widths  of 
long  rectangular  strips,  which  is  a  matter  of  great  concern 
in  the  fabrication  of  semiconductor  components.  We  will 
show  here  both  experimentally  and  theoretically  the  ad¬ 
vantages  of  utilizing  phase  to  determine  the  linewidth.  We 
will  also  introduce  here  a  Fourier  transform  technique  for 
linewidth  measurement  in  cases  where  the  strip  width  is 
comparable  to  the  spot  size. 

1)  Depth  Profiling:  Topography  images  of  metallized 
stripe  patterns  on  a  fused-silica  substrate  are  shown  in 
Figs.  5(a),  5(b),  and  5(d).  The  patterns  have  progressively 
finer  pitches  and  the  linewidths  are  250,  125,  and  62.5 
pm,  respectively.  The  metallization  is  gold  with  a  thick¬ 
ness  of  about  3000  A .  The  gold  stripes  show  up  as  bright 
areas  in  the  images.  At  50  MHz.  with  an  f-numberof  0.85 
and  a  uniformly  excited  aperture,  the  Rayleigh  resolution 
of  the  acoustic  lens  is  1.13  FX  =  0.96  X  or  29  pm,  while 
the  3-dB  resolution  is  0.64  FX  =  0.54  X  or  16  /*m.  The 
62.5-^m  line  image  in  Fig.  5(d)  is  clearly  resolved,  as 
would  be  expected.  The  amplitude  image  of  the  \25-pm 
stripes  is  shown  in  Fig.  5(c)  for  comparison.  There  is  vir- 
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250  f. im 

(e) 

Fig.  5.  Topography  images  of  striped  metallization  patterns  on  a  fused 
quartz  substrate.  The  metallization  is  gold  and  the  thickness  is  about  3000 
A .  (a)  Phase  image  of  250-/im  lines,  (b)  Phase  image  of  125-^m  lines, 
(c)  Amplitude  image  of  125-pm  lines,  (d)  Phase  image  of  62.5-ftm  lines, 
(e)  Perspective  plot  of  the  measured  topography  map  of  the  250- jim  line 
pattern. 

tually  no  contrast  to  surest  the  presence  of  a  striped  pat¬ 
tern,  clearly  illustrating  that  phase  is  a  far  more  sensitive 
means  of  gauging  distances.  The  measured  phase  differ¬ 
ence  between  the  surface  of  the  metallization  and  the 
surface  of  the  fused-silica  substrate  is  8.3°,  which  corre¬ 
sponds  to  a  film  thickness  of  2800  A  after  taking  into 
account  the  angular  dependence  of  the  complex  reflec¬ 
tance  function  and  the  effect  of  focusing. 

We  will  show  here  that,  based  on  theoretical  consider¬ 
ations,  indeed  the  phase  rather  than  the  magnitude  of  the 
acoustic  microscope  output  V(x)  is  a  more  sensitive  mea¬ 
sure  of  the  depth  profile.  Consider  the  simple  case  of  a 
rectangular  strip  of  width  w  on  a  substrate  of  identical 
material.  Assume  that  the  reflection  coefficient  T0  of  the 
surface  is  real.  Then  the  spatial  response  of  the  strip 
structure  is  given  by 

R(x)  =  r„  +  T0  ( e‘ *  -  1)  rect  (jc/w)  (1) 


where  <£  is  the  phase  change  due  to  the  thickness  of  the 
strip.  Let  the  response  of  the  imaging  system  to  a  line 
source  be  h(x).  Furthermore,  suppose  that  the  lens  is  uni¬ 
formly  illuminated  in  the  back  focal  plane  and  the  lens  is  ' 
aberration-free  so  that  h(x)  is  real.  The  acoustic  micro¬ 
scope  output  is  therefore  given  by 

V(x)  =  R(x)  *  h(x)  (2) 

or 

V(x)  =  r 0K  +  r0  ( e ’*  -  1)  six)  (3) 

where 

s(jc)  =  h(x)  *  rect  ( x/w )  (4) 

and 

i+oo 

h(x)  dx. 

—  OO 

The  first  term  on  the  right-hand  side  of  (3)  corresponds  to 
the  constant  background  reflection  from  the  substrate,  and 
the  second  term  corresponds  to  the  additional  spa¬ 
tially  varying  contribution  that  is  due  to  the  strip.  It  can 
easily  be  shown  that  for  a  thin  strip  where  <t>  «  1,  to 
second  order  in  <j>,  the  magnitude  and  the  phase  of  V(x) 
are  given  respectively  by  the  relations 

w*)  =  r 0K  +  [Six)  -  K]  4>2  (5) 

and 

w*)  =  £*<*>•  (6) 

The  Fphasefx)  is  directly  proportional  to  <j>,  whereas  the 
spatially  varying  part  of  Fmag(x)  depends  on  <t>  to  the  sec¬ 
ond  order.  Thus  for  thin  strips,  phase  is  a  much  more  sen¬ 
sitive  measure  of  the  depth  profile. 

2)  Linewidth  Measurement:  A  major  advantage  of 
scanned  confocal  imaging  systems  such  as  ours  is  that  be¬ 
cause  the  point  spread  function  is  always  positive  and  falls 
off  rapidly  at  large  distances,  there  is  hardly  any  ripple  in 
the  amplitude  of  the  step  response  [16],  [17].  A  similar 
smooth  transition  in  the  phase  of  the  step  response  can 
also  be  observed  in  the  perspective  plot  shown  in  Fig.  5(e). 
For  the  purpose  of  comparison,  we  have  carried  out  a  nu¬ 
merical  analysis  for  the  step  phase  responses  of  different 
imaging  configurations.  For  small  step  heights,  or  equiv¬ 
alently  small  phase  changes,  the  normalized  phase  step 
response  is  of  the  form 

/3(jc)  =  |  g(r)  cos’1  (x/r)  r  dr  (7) 

where  x  is  the  distance  between  the  center  of  the  beam 
and  the  step,  and  g(r)  is  the  point  response  function  of 
the  imaging  system.  It  is  assumed  that  regions  on  either 
side  of  the  step  are  of  equal  reflectivity.  The  step  response 
function  is  plotted  in  normalized  form  in  Fig.  6  for  two 
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distance;  in  wavelengths 

ig.  6.  Slep  phase  response  of  a  normal  phase  interference  microscope 
[solid  curve)  and  a  confocal  imaging  system  (theory  in  dashed  curve, 
experimental  result  in  closed  circles). 


afferent  cases.  The  solid  curve  corresponds  to  g(r)  = 
inc  (2  r  sin  0/X),  which  is  the  point  response  of  a  normal 
ihase  interference  microscope  where  the  illumination  is 
i  plane  wave.  Here  6  is  the  maximum  half  angle  of  the 
objective  lens;  X  is  the  wavelength  of  the  illumination;  and 
line  (x)  is  defined  as  2/,(irx)/irx.  Note  that  in  this  case 
>(r)  reverses  in  sign  and  its  corresponding  step  phase  re¬ 
sponse  does  not  change  monotonically  with  distance.  The 
dashed  curve  corresponds  to  g(r )  =  jinc2  (2r  sin  0/X), 
the  point  response  of  a  confocal  system.  The  confocal  re¬ 
sponse  shows  a  distinct  lack  of  ripples.  The  experimen¬ 
tally  measured  step  phase  response  is  also  plotted  in  closed 
circles  showing  close  agreement  with  theory.  In  both  im¬ 
aging  configurations,  the  50-percent  threshold  of  the  phase 
profile  demarcates  the  location  of  the  edge.  It  is  obvious, 
however,  that  the  smooth  transition  of  the  confocal  step 
phase  response  enables  a  straightforward  and  unambigu¬ 
ous  determination  of  not  only  the  strip  thickness  but  also 
the  strip  width. 

The  50-percent  thresh  •  i  criterion  used  to  locate  the 
edges  of  a  strip  for  estimating  the  linewidth  is  a  valid  one 
for  wide  strips.  For  strip  widths  comparble  to  the  spot  size 
of  the  imaging  system,  the  step  responses  of  the  two  edges 
are  in  close  proximity  and  tend  to  interfere  with  each 
other,  making  the  establishment  of  a  general  criterion  for 
locating  the  strip  edges  difficult.  We  will  describe  here  a 
Fourier-transform  method  that  circumvents  these  difficul¬ 
ties. 

Applying  Fourier  transform  to  (6)  and  using  the  defi¬ 
nition  of  ,v(x)  in  (4),  the  spatial  frequency  spectrum  of 
|/pha*(/J  is  given  by 

^ph **(/,)  =  ^  //(/,)  •  »v  sine  ( wf,)  (8) 

where  H(ft)  is  the  Fourier  transform  of  h(x).  The  term 
with  the  sine  function  corresponds  to  the  Fourier  trans¬ 


Fig.  7.  Spectral  responses  of  5000- A -thick  gold  lines  on  a  fused-quartz 
substrate.  The  nominal  linewidths  are  76  /jm,  43  pm.  and  33  |im.  re¬ 
spectively. 


form  of  the  spatial  variation  of  the  strip  and  its  zeroes  are 
located  at 


/,  =  nAv,  n  =  I,  2,  3.***.  (9) 

The  term  //(/,)  represents  the  angular  spectral  response 
of  the  imaging  system.  For  a  confocal  system  \H(  fx)\  is 
maximum  at/,  =  0  and  generally  falls  off  monotonically, 
until  it  becomes  zero  at  the  upper  cutoff  frequency  of 

/*..  =  1W> 

where  F  is  the  f-number  of  the  lens,  and  X  is  the  wave¬ 
length  of  the  illumination  [16],  [17).  Hence  the  zeroes  of 
Kphasc(/t)  are  those  of  the  strip  response  within  the  pass- 
band  of  the  imaging  system.  As  long  as 

vv  >  FX  (10) 

the  strip  width  w  can  be  determined  using  (8)  from  the 
locations  of  the  zeroes  of  FphaSe(/)-  In  our  case  F  =  0.85, 
so  in  principle  we  can  determine  linewidths  as  small  as 
0.85  X  or  25  p m  at  50  MHz  in  water. 

The  merit  of  this  Fourier-transform  technique  lies  in  its 
simplicity.  Given  a  priori  knowledge  of  the  strip  geome¬ 
try,  one  does  not  need  to  know  the  exact  nature  of  imaging 
system.  The  mere  location  of  the  zeroes  for  the  spectral 
response  of  the  object  suffices  for  the  accurate  determi¬ 
nation  of  the  linewidth. 

Fig.  7  illustrates  the  experimental  |Fpkl,c(/,)|  for  >v  = 
76  )zm,  43  |im,  and  33  gm,  respectively.  The  strips  are 
5000- A -thick  gold  lines  on  a  fused-quartz  substrate.  The 
acoustic  microscope  operates  at  50  MHz  (X  =  30  ^m) 
with  F  =  0.85,  and  the  acoustic  beam  is  focused  on  the 
surface  of  the  fused-quartz  substrate.  The  IW/,)I  is 
plotted  on  a  log  scale.  The  weighting  due  to  the  imaging 
system  response  |//(/,)|  ail  the  way  out  to  the  resolution 
limit  is  very  much  in  evidence.  The  linewidths  estimated 
from  the  position  of  the  first  zero  of  ^hasc (/.)  are  76.3, 
43.5,  and  31.4  /zm.  respectively.  They  agree  very  well  with 
the  optically  measured  widths.  The  minor  discrepancies 
are  mainly  due  to  the  nonuniformity  in  the  widths  of  the 
gold  lines. 

In  addition,  this  Fourier  transform  approach  is  essen¬ 
tially  insensitive  to  defocusing.  It  is  obvious  from  (8)  that 
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Fig.  8.  Spectral  responses  of  2000- A -thick  76-/tm-wide  chrome  line  on  a 
fused  quanz  substrate  at  various  defocusing  distances:  z  =  0  ^m.  +60 
/im.  and  +80  (im.  respectively. 

defocusing  only  changes  the  spectral  response  of  the  im¬ 
aging  system  H(  fx)  slightly  and  does  not  introduce  any 
new  zeroes  to  nor  does  it  affect  the  zeroes  of  |  Vvhiiffx)\. 
Fig.  8  shows  the  experimental  |Fphase(/r)l  for  a  2000- A  - 
thick  and  76-/zm-wide  chrome  line  on  fused  quartz  at  var¬ 
ious  defocusing  distances.  The  focus  was  located  on  the 
substrate,  and  60  pm  and  80  pm  above  the  substrate,  re¬ 
spectively.  The  locations  of  the  first  zero  in  all  three  cases 
are  essentially  coincident,  giving  a  width  estimate  of 
76.3  gin. 

3)  Independent  Measurement  of  Velocity  Perturbation 
and  Topography:  In  conventional  acoustic  micrographs, 
surface  topography  and  material  property  both  contribute 
to  contrast  in  the  image,  and  their  respective  effects  are 
generally  indistinguishable  from  one  another.  With  the  to¬ 
pography  tracking  mechanism  described  in  Section  II-2 
built  into  our  phase  measurement  system,  surface  topog¬ 
raphy  and  material  property  change  can  essentially  be  ob¬ 
tained  independent  of  each  other. 

Fig.  9(a)  is  the  topographic  line  scan  of  an  aluminum 
film  step  on  glass  with  a  nominal  thickness  of  5000  A . 
The  initial  and  trailing  slow  phase  changes  are  due  to  the 
large  spatial  extent  of  the  reference  signal  as  it  traverses 
the  step,  whereas  the  abrupt  phase  change  is  due  to  the 
focused  beam  as  it  crosses  the  step.  Fig.  9(b)  is  the  same 
profile  on  an  expanded  -;cale.  The  phase  change  caused  by 
the  step  is  11.25°,  which  corresponds  to  a  thickness  of 
4688  A.  With  a  limiting  phase  resolution  of  0.1°,  the  ul¬ 
timate  height  sensitivity  of  the  system  is  about  50  A . 

Fig.  9(c)  shows  a  velocity  perturbation  scan  of  the  same 
sample.  The  phase  change  measured  is  7.75°,  correspond¬ 
ing  to  a  velocity  perturbation  of  0.25  percent,  which  is  in 
fair  agreement  with  the  theoretically  calculated  value  of 
0.21  percent. 
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Fig.  9.  I. 'dependent  measurements  of  velocity  perturbation  and  topography 
profile  due  to  a  5000- A  aluminum  film  step,  (a)  Topography  profile  scan, 
(b)  Topography  profile  scan  on  an  expanded  scale,  (c)  Velocity  pertur¬ 
bation  sen. 


curve  to  determine  only  the  Rayleigh-wave  velocity.  The 
solid  curves  in  Fig.  10  show  the  magnitude  and  the  phase 
of  the  theoretical  reflectance  function  of  a  water-fused- 
silica  interface.  The  reflectance  function  contains  infor¬ 
mation  on  the  longitudinal  and  shear  critical  angles,  and 
thus  the  velocities  of  propagation  of  the  respective  modes 
in  fused  silica.  Also,  the  Rayleigh  critical  angle  corre¬ 
sponds  to  the  point  at  which  the  phase  is  x  radians  in  the 
region,  where  the  phase  curve  undergoes  a  rapid  2x  ra¬ 
dian  change. 

It  can  be  shown  that  the  relation  between  K(z)  and  the 
reflectance  function  R(6)  is  essentially  one  of  Fourier 
transformation  [10].  The  inversion  formula  is  given  by 

J  (F(m)  exp  [jitu/sm2  (a/2)]} 

R0)  ~  [  t/.+  (0]2  P'V)  [1-2  t  sin2  (a/2)] 


C.  Rejiectance  by  Inversion  of  Complex  V(z)  where 

It  has  been  demonstrated  that  the  reflectance  function 
at  a  liquid-solid  interface  can  be  obtained  by  inverting  the 
corresponding  complex  V(z)  data  [10].  The  principle  mo¬ 
tivation  behind  this  endeavor  is  that  one  can  obtain  much  an(j 
more  useful  material  property  information  from  the  re¬ 
flectance  function  than  the  usual  and  rather  limited  treat¬ 
ment  of  measuring  the  periodicity  of  the  nulls  in  the  V(z) 


_  4;  sin2  (a/2) 
X 


sin2  (6/2) 
sin2  (a/2) 
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FUSED  SILICA 
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Fig.  10.  Theoretical  and  experimental  reflectance  function  for  a  water- 

fused-quartz  interface  obtained  by  inversion  of  the  corresponding  com¬ 
plex  P(z). 

The  a  is  the  maximum  half  angle  of  the  lens;  6  is  the 
incident  angle  with  respect  to  the  liquid-solid  interface; 
and  X  is  the  wavelength  of  the  acoustic  wave  in  the  liquid 
medium.  In  addition,  U\  is  the  back  focal  plane  illumi¬ 
nation,  and  P  is  the  pupil  function  of  the  lens.  Generally 
these  two  lens  parameters  are  not  known,  but  [Uf]2  P2  in 
the  denominator  of  (11)  can  be  calibrated  by  measuring 
the  V(z)  for  a  material  whose  reflectance  function  has  uni¬ 
form  amplitude  and  phase  over  the  range  of  angular  spec¬ 
tral  components  excited  by  the  acoustic  lens. 

The  measurements  were  conducted  at  10  MHz.  A 
water-lead  interface,  in  which  the  reflectance  function  is 
uniform  in  magnitude  and  phase  out  to  6  =  40°,  was  used 
to  calibrate  the  [C/f]2  P2  term.  The  experimental  reflec¬ 
tance  function  of  fused  silica,  obtained  from  the  inversion 
procedure,  is  superposed  on  Fig.  10  in  dashed  lines.  There 
is  good  agreement  between  theory  and  experiment  in  both 
magnitude  and  phase.  The  measured  shear  critical  angle 
at  0  =  23.5°  and  the  Rayleigh  critical  angle  at  6  =  25.85° 
compare  extremely  well  with  theoretical  values. 

IV.  Conclusion 

We  have  demonstrated  that  phase  information  is  a  use¬ 
ful  asset  in  acoustic  microscopy.  Phase  can  be  used  to  pro¬ 
vide  sensitivity  and  information  in  surface  material  char¬ 
acterization  unparalleled  by  amplitude-only  measurement 
techniques.  Most  of  the  concepts  described  in  this  work 
have  complete  generality  and  can  readily  be  used  in  other 
applications.  In  fact,  many  ideas  presented  in  this  work 
have  been  applied  to  a  new  scanning  optical  microscope 
[19]  developed  by  Jungerman  and  Kino  for  surface  topog¬ 
raphy  measurement.  In  the  optical  case  both  the  amplitude 
and  the  phase  of  the  reflected  beam  can  be  measured  as 
in  the  acoustic  microscope.  Preliminary  results  from  the 
scanning  optical  microscope  give  further  confirmation  of 
most  of  the  concepts  put  forth  in  this  paper. 
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Abstract— It  is  demonstrated  that  the  reflectance  function  R(8)  of  a 
liquid-solid  interface  can  be  obtained  by  inverting  the  complex  V(z)  data 
collected  with  an  acoustic  microscope.  The  inversion  algorithm  is  based 
on  a  nonparaxial  formulation  of  the  F(j)  integral,  which  establishes  the 
Fourier  transform  relation  between  R($)  and  F(j).  Examples  are  given 
to  show  that  with  this  measurement  technique,  the  acoustic  phase  ve¬ 
locities  of  the  propagating  modes  in  the  solid  medium  can  easily  be 
determined  and  material  losses  can  be  estimated.  The  same  technique 
is  also  used  for  characterizing  imaging  performance  of  focused  sys¬ 
tems.  Applications  in  thin-fllm  measurement  are  also  discussed. 

I.  Introduction 

NEW  TECHNIQUE  for  measuring  and  interpreting 
the  acoustic  material  signature  or  the  so-called  Viz) 
in  acoustic  microscopy  will  be  presented  in  this  paper. 
V(z)  refers  to  the  variation  of  the  normalized  voltage  out¬ 
put  of  an  acoustic  microscope  as  a  function  of  the  sepa¬ 
ration  distance  z  between  the  focal  point  of  the  acoustic 
lens  and  the  plane  surface  of  a  reflecting  object  [1]— [4]. 
This  variation  of  voltage  response  with  distance  is  char¬ 
acteristic  of  the  material  being  examined  and  is  an  im¬ 
portant  source  of  contrast  in  acoustic  micrographs.  The 
Vfz)  effect  is  a  well-understood  phenomenon,  and  various 
models  have  been  put  forth  to  explain  how  it  is  related  to 
acoustic  material  properties  [1],  [2],  [5],  [6],  [9].  So  far, 
practical  application  revolves  around  the  inverse  process 
of  determining  acoustic-property-dependent  parameters  by 
interpreting  selective  features  of  the  magnitude  of  the  Viz) 
function.  For  example,  the  leaky  surface  acoustic  wave  ve¬ 
locity  can  be  obtained  directly  from  the  periodicity  of  the 
interference  peaks  and  nulls  of  the  magnitude  of  the  Viz) 
function  in  the  negative  z  region  [4]-[7],  where  the  object 
surface  is  located  between  the  focal  point  and  the  lens. 
Another  approach  invol\j->  the  Fourier  transform  of  the 
negative  z  portion  of  the  magnitude  of  the  Viz)  function  in 
order  to  identify  propagating  modes  at  the  liquid-solid  in¬ 
terface  and  to  determine  their  respective  phase  velocities 
[81. 

The  basic  theory  and  experimental  result  of  a  more  gen¬ 
eral  Fourier  transform  technique,  which  utilizes  both  the 
magnitude  and  phase  of  the  Viz)  function  for  material 
property  characterization,  will  be  given  in  this  paper.  The 
theoretical  treatment  is  a  combination  of  angular  spectra 
and  ray  optics  formalisms.  It  is  carried  out  specifically  for 
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the  case  where  the  acoustic  transducer  is  a  spherical  shell 
that  radiates  directly  into  the  liquid  coupling  medium  and 
hence  forms  a  focusing  lens  by  itself.  An  integral  expres¬ 
sion  relating  the  complex  Viz)  and  the  reflectance  R(6)  of 
the  liquid-solid  specimen  interface  as  a  function  of  the 
angle  6  between  the  direction  of  an  incident  plane  wave 
and  the  normal  to  the  interface  surface  is  derived.  It  will 
be  shown  that  by  a  suitable  transformation  of  variables, 
R(6)  can  be  obtained  by  an  inverse  Fourier  transform  of 
V{z). 

The  primary  motivation  behind  this  work  is  that  the 
complex  reflectance  function  R(0)  contains  a  wealth  of 
readily  available  acoustic  material  property  information. 
The  solid  curves  in  Fig.  7  show  the  theoretical  magnitude 
and  phase  of  the  reflectance  function  of  a  water-fused- 
silica  interface.  The  longitudinal  and  shear  critical  angles, 
9U  and  0&,  show  up  as  distinctive  features  on  the  magni¬ 
tude  curve  and  their  locations  give  the  propagation  veloc¬ 
ities  of  the  respective  modes  directly.  In  addition  the  Ray¬ 
leigh  critical  angle  dRc  corresponds  to  the  point  where  the 
phase  curve  becomes  it  in  a  region  where  it  undergoes  a 
rapid  2%  radian  change.  Also,  as  will  be  shown  in  Section 
IV,  the  effect  of  loss  in  the  solid  medium  can  readily  be 
observed  in  and  estimated  from  the  reflectance  function 
Rid).  There  are  other  potential  areas  of  application  in¬ 
volving  film  thickness  measurement  and  the  thin-film 
acoustic  matching  between  two  media  of  vastly  different 
impedances,  which  will  also  be  explored  in  Section  IV. 

II.  Nonparaxial  Formulation  of  the  V{z)  Integral 

The  derivation  of  the  relation  between  k'fz)  and  R  (9)  in 
the  literature  [1],  [2],  [6],  [9]  can  be  classified  into  two 
general  categories:  those  involving  Fourier  spectral  de¬ 
composition  [1],  [2],  [6)  and  those  based  on  ray  optic 
models  [9].  The  angular  spectrum  approach  is  more  ame¬ 
nable  to  inversion  since  Viz)  can  be  expressed  as  an  inte¬ 
gral  involving  R(0),  and  the  integral  can  easily  be  manip¬ 
ulated  to  establish  a  Fourier  transform  relation  between 
K(z)  and  Rid).  The  ray  optic  approach  essentially  regards 
Viz)  as  the  interference  between  a  family  of  specularly 
reflected  rays  and  a  family  of  laterally  displaced  leaky  sur¬ 
face  acoustic  wave  rays.  It  is  not  complete  in  that  it  ne¬ 
glects  the  contributions  from  regions  near  the  longitudinal 
and  shear  critical  angles.  Also,  since  it  does  not  take  dif¬ 
fraction  into  account,  it  is  not  valid  near  the  focal  plane. 
In  addition,  the  ray  optic  formalism  in  the  form  usually 
given  is  incompatible  with  inversion  procedures  for  find¬ 
ing  Rid).  A  common  drawback  to  both  kinds  of  formula- 
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Fig.  1.  Schematic  of  spherical  focused  transducer  used  in  t'tz)  measure¬ 
ments. 


tions  is  that  they  both  make  use  of  the  paraxial  approxi¬ 
mation. 

An  alternative  hybrid  formulation  of  the  relation  be¬ 
tween  V(£)  and  R(6)  will  be  given  in  the  following  treat¬ 
ment.  This  approach  combines  features  from  the  angular 
spectrum  as  well  as  the  ray  optic  models  and  is  valid  for 
all  z  including  the  focal  plane.  Moreover,  no  paraxial  as¬ 
sumption  has  to  be  made  in  its  development.  The  special 
case  of  a  focused  spherical  transducer  radiating  directly 
into  water  with  no  buffer  medium,  as  shown  in  Fig.  1,  is 
considered.  The  practical  significance  of  this  acoustic  lens 
configuration  will  be  discussed  in  Section  IV.  The  elimi¬ 
nation  of  the  buffer  rod  simplifies  the  derivation.  However, 
the  final  solution  is  valid  for  any  lens  configuration  that 
generates  a  converging  spherical  wave,  and  can  be  gen¬ 
eralized  for  use  with  other  systems. 

l/sing  the  notation  given  in  Fig.  1,  the  potential  or 
pressure  field  emitted  by  a  transducer  of  diameter  2  r0  is 
equivalent  to  that  produced  by  a  perfect  converging  spher¬ 
ical  wave  impinging  on  a  circular  aperture  of  radius  r0. 
All  acoustic  field  terms  will  be  assumed  to  vary  as  exp 
(jut).  Referring  to  the  aperture  plane  P  in  Fig.  1,  the  po¬ 
tential  field  distribution  in  cylindrical  coordinates  is  given 
by  the  relation 


<Mr)  =  P(r) 


exp  [jkJr 2  +  f]\ 


(i) 


where  f,  =  Jfl  -  ri  o;  =  /0  cos  0;  f0  is  the  focal 
length;  0Q  is  the  maximum  half  angle  subtended  by  the 
aperture  at  the  focus;  and  k  is  the  propagation  constant  in 
the  liquid  medium.  P(r)  is  the  pupil  function  of  the  lens, 
which  in  the  most  general  case  takes  into  account  the  non- 
uniform  excitation  of  the  transducer  and  any  aberration 
effects.  In  the  simplest  case  for  an  ideal  lens 


P(r)  = 


0  <  r  <  r0 
r  <  r0 


(2) 


Making  use  of  the  circular  symmetry  of  the  field  distri¬ 
bution,  <t>(r)  can  be  decomposed  into  cylindrical  radial 
waves  of  wave  number  kr  by  applying  the  Hankel  trans¬ 
form 


A(kr)  =  ( 
Jo 


P{r) 


exp  [jkjr  +  /,] 


Jtiiks)  r  dr.  (3) 


Here  /c2  =  Ic^  +  k?r.  Since  the  transducer  diameter  is  nor¬ 
mally  many  wavelengths  wide  and  the  integrand  is 
weighted  by  r,  the  contribution  to  the  integral  by  J0(krr ) 
comes  mainly  from  regions  where  krr  is  large.  Therefore 
A(kr)  can  be  approximated  by  replacing  J^krr)  with  its 
large  argument  asymptotic  form 


Jo(krr)  ~ 


./  2r\ 

/  7r\ 

exp 

It  _  v\ 

+  exp 

\l2irk,r 


(4) 


Thus,  it  follows  that 


A(kr) 

where 


rn,'«P  Wi(r>l  +  exP  Uh(r)]  A 
Jo  W  Jl,k,r(?  +  f\) 

h(r)  =  (kjr2  +  f\  -  krr  +  ^ 
fair)  =  (kJr  +  f\  +  krr  - 


(5) 


The  method  of  stationary  phase  will  now  be  applied  to 
evaluate  (5)  [10].  Consider  the  phase  of  the  first  term  of 
the  integrand  (r).  The  contribution  to  the  integral  comes 
from  the  region  near  where  > p\(r)  -  0 


nr) 


kr 


s  -  K  =  0 


or  where  r  =  r' 


r'  = 


kr 


/.  = 


(6) 


Note  that  r'  is  positive  in  this  case. 

Next,  consider  the  phase  of  the  second  term  of  the  in¬ 
tegrand  4/2 (r).  Setting  ^(r)  to  zero  yields  the  result 


kr' 


=  ~kr 


(7) 


which  implies  that  r'  is  negative.  Since  in  a  cylindrical 
coordinate  system  r  is  always  positive,  there  is  no  contri¬ 
bution  to  A(kr)  from  the  second  term  of  the  asymptotic 
form  of  J0(krr). 

Equation  (5)  therefore  becomes 


A{kr)  = 


(  2t  y/2n,_„  exp  [jknn]  kr 
Wro- ')!/  V27T krr'(r'2  +  f\)  k; 


(8) 


where  r'  is  given  by  (6)  and 

4/\(r')  = 


kr'2 


Jr'1  +  f\  (r'2  +  f])m 
Substituting  (6)  into  (8)  gives  the  result 

A(kr)  =  P(kr)  k  exp  \jk.ft]/kt 


k> 

&f\ 


(9) 
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where 


kz  =  V/fc2  -  T*. 


A(kr ;  0)  =  P(fc,)  k/kz 


A(kn-  z)  =  /*(*,)  k  exp  [  — y*zz]/fcz. 


z)  =  /?(£,)  /4(kr;  z). 


Viz)  =sn=J- 


1  (0V.  -  0'  I 
J  J 


0'mJ) 


j  4v; 


</S 


1  dz 

,  30r 

M,  =  — . 


OZ 


tW  =  ~jk.A(kr) 
W.(kr)  =  A-fl(^r). 


By  applying  Parseval’s  Theorem  and  substituting  in  (16), 
(14)  can  be  written  in  the  form 


Viz)  = 


Aik/,  z)  Bik/  z)  k.kr  dkr 


Aik/,  Z)  A*(kr;  z)  k,kr  dkr 


Substituting  (12)  and  (13)  into  (17),  the  final  result  is 


(10) 


To  find  the  corresponding  angular  spectrum  at  a  plane 
of  distance  z  +  /,  from  the  aperture  plane  P,  each  spectral 
component  is  multiplied  by  the  propagation  factor  exp 
[ ~jkz (z  +  /i)].  For  instance,  at  the  focal  plane  z  -  0 


Viz)  = 


PU)  «(*,)  e-^lz2Mkr<JK 


P  ikr) 


(18) 


kr  dkr 


(ID 


which  simply  means  that  with  an  ideal  pupil  function  all 
angular  spectral  components  arrive  in  phase  at  the  focal 
point.  In  general,  for  any  arbitrary  z 


By  making  the  following  substitution  of  variables,  a 
comparison  can  be  made  between  (18)  and  Atalar’s  par¬ 
axial  angular  spectrum  formulation  of  the  Viz)  integral 
[1],  Let 


(12) 


kjk  =  r/f0. 

Equation  (18)  can  be  written  in  the  form 


Suppose  the  plane  reflector  specimen  is  placed  at  z,  and 
the  reflectance  function  of  the  liquid-specimen  interface 
is  Rikr).  The  reflected  angular  spectrum  at  z  is 


/>  V/0)  Pir/fo) 


Viz)  = 


exp  [-jlkz'J  1  -  (r//0)2] 
Vl  ~  ir/fo)2 


(19) 


r  dr 


P\r/fQ) 


r  dr 


(13) 


Vl  -  ir/fo)1 


The  reciprocity  principle  of  Kino  and  Auld  given  in  the 
Appendix  is  employed  to  find  the  normalized  voltage  re¬ 
sponse  V(z)  of  the  focused  transducer  to  the  reflected  wave 
in  which  the  spectrum  is  given  by  (13).  V{z)  is  equal  to 
the  scattering  coefficient 


(14) 


where  uz  is  the  normal  displacement  at  the  liquid-solid 
interface,  and  the  integration  is  evaluated  over  the  same. 
The  superscripts  /  and  r  designate  the  incident  and  the 
reflected  fields  respectively,  and 

,  _  30' 


(20) 

which  differs  from  Atalar’s  solution  by  the  factor 
1/Vl  -  ir/f0)2  in  the  denominator  of  each  integrand.  In 
the  limit  of  the  paraxial  approximation  where  ir/f0)2  is  as¬ 
sumed  to  be  small,  this  additional  factor  tends  to  unity, 
and  (18)  becomes  the  expression  derived  by  Atalar. 

III.  Fourier  Transform  Relation  Between 
Viz)  and  Rid) 

The  possibility  of  obtaining  Ri8)  by  inverting  V{z)  has 
been  alluded  to  in  optical  microscopy  [11],  and  the  inver¬ 
sion  has  been  demonstrated  in  acoustic  microscopy  [12] 
based  on  Atalar’s  paraxial  K(z)  integral.  The  Fourier 
transform  relation  between  Viz)  and  Rid)  for  the  nonpar- 
axial  formulation  given  in  (18)  can  be  established  as  fol¬ 
lows.  Since  k?z  =  k2  -  k*r,  it  follows  that 


k.dk.  =  —krdkr 
Equation  (18)  becomes 


(15) 


■>kcos$o 


For  the  circularly  symmetric  incident  and  reflected  radial 
wave  components,  the  Hankel  transforms  of  the  «.  com¬ 
ponents  are 


K(z)  =  C  P2ikz)R(k:)  exp  i~j2kzz)  dkz 
where  C  is  the  normalization  factor  given  by 

[pkcosOo  -|  - 1 

P2ik.)  dk: 


(21) 


(22) 


(23) 


(16) 


Noting  that  kz  =  k  cos  d,  where  k  =  2i r/\,  and  making 
the  following  change  of  variables 


z 

M  =  x 


t  =  2  cos  d 


(24) 


(17) 


we  obtain  the  result 

« 2  cos  fa 


Viu) 


-j; 


P2it)  R{t)  exp  (—j2irut)  dt.  (25) 


svm  y^jooi  Asziaflac si  aasaaaa  aaa 


216 


IEEE  TRANSACTIONS  ON  SONICS  AND  ULTRASONICS,  VOL,  SU-32.  NO.  2,  MARCH  1985 


In  short-hand  notation 


IV.  Experimental  Results 


m  =  5{P2(t)  R(t)}  (26) 


where  {  }  denotes  the  Fourier  transform.  The  inversion 
to  obtain  R  ( t )  is  given  by 


R(t)  = 


P\t f 

0. 


2  >  t  >  2  cos  90 
otherwise. 


(27) 


Note  that  R(t)  is  only  determined  for  the  angular  range  of 
insonification  covered  by  the  pupil  function  of  the  lens.  In 
general  P(t)  is  not  known  a  priori.  However,  Pi(t)  can 
conveniently  be  calibrated  by  inverting  the  V(z)  function 
for  a  liquid-solid  interface  in  which  the  reflectance  has 
uniform  magnitude  and  phase  over  the  angular  extent  of 
the  lens  pupil  function.  In  particular,  if  a  perfectly  reflect¬ 
ing  surface  with  R(t)  =  1  is  used 


P\t)  =  $F-'{Ko(m)}  (28) 


where  F0(u)  is  the  V(z)  of  the  “ideal”  reflector. 

This  inversion  approach  to  material  characterization  is 
different  in  that  it  requires  the  measurement  of  V(z)  in  both 
amplitude  and  phase  whereas  previous  work  [  1  ]— [8]  uti¬ 
lized  only  the  amplitude  or  intensity  of  F(z).  It  is  generally 
not  trivial  to  conduct  high-accuracy  phase  measurement 
in  acoustic  microscopy.  A  special  phase-measurement  sys¬ 
tem  [13]  has  been  developed  to  enable  the  precise  acqui¬ 
sition  of  phase  data  and  the  subsequent  experimental  dem¬ 
onstration  of  the  inversion  technique. 

It  is  interesting  at  this  point  to  examine  the  range  res¬ 
olution  of  a  confocal  imaging  system  such  as  the  reflection 
acoustic  microscope  for  an  ideal  plane  reflector;  i.e.,  R(6) 
=  1.  Assuming  an  ideal  lens  pupil  function  of  the  form  in 
(2),  the  voltage  response  of  the  acoustic  microscope  as  a 
function  of  defocusing  distance  can  be  calculated  from 
(25)  to  be 

J  2  cos  0o 

exp  [-j2irut\  dt 
or 

V(u)  =  -  C  exp  [~jiru{  1  +  cos  0O)1 

■  sine  [2«(1  —  cos  0O)].  (29) 
The  3-dB  point  of  |  (^(m) |  is  given  by  the  relation 
2(1  -  cos  0())  uMB  =  0.457T 


or 


Z)<JB 


0.225X 
1  -  cos  dn 


Therefore  the  3-dB  range  resolution  or  the  focal  depth  for 
a  plane  reflector  is 


0.45X 

1  -  COS  Or)' 


(30) 


In  the  V(z)  measurement  a  special  acoustic  transducer 
is  used.  The  piezoelectric  element  is  a  spherical  shell  with 
an  air  backing  and  a  quarter-wave  matching  layer  on  the 
front  face.  The  transducer  radiates  directly  into  water,  and 
thus  forms  an  aberration-free  spherical  lens  that  focuses 
the  emitted  acoustic  beam  into  a  diffraction-limited  spot 
at  the  center  of  curvature  of  the  transducer.  This  acoustic 
lens  configuration  is  advantageous  for  V(z)  measurements 
because  the  elimination  of  the  buffer  rod  removes  the  sys¬ 
tematic  noise  due  to  reverberations  inside  the  buffer  rod 
that  tend  to  overlap  in  time  with  the  reflected  signal  from 
the  object.  With  such  a  construction  the  acoustic  lens  can 
be  translated  over  a  wide  range  of  z  to  collect  V(z)  data  in 
an  essentially  noise-free  environment.  The  maximum  limit 
to  the  vertical  movement  is  imposed  by  how  far  the  lens 
can  be  moved  towards  the  specimen.  The  problem  with 
multiple  echoes  in  the  water  path  between  the  transducer 
and  the  specimen  is  overcome  by  exciting  the  transducer 
with  a  short  tone  burst  and  time-gating  the  first  arrival 
return  signal. 

Two  transducers  of  this  type  with  F/\  and  F/0.1  aper¬ 
tures,  respectively,  are  used  in  the  experiment.  The  radius 
of  curvature  or  focal  length  /0  is  16  mm  in  both  cases. 
From  geometric  considerations,  the  maximum  possible 
negative  defocusing  distance  is  given  by 

zmin  =  -/.  =  -/oVl  -  1/(2F)2  (31) 

where  F  is  the  /-number  of  the  lens.  For  the  FI 0.7  lens 
|  Zmin  I  -  11-2  mm,  and  for  the  F/l  lens  |zmin|  =  13.9  mm. 
The  importance  of  zmin  will  be  discussed  in  Section  V. 

Both  transducers  have  a  fundamental  frequency  of  ap¬ 
proximately  3  MHz.  For  V(z)  measurements,  they  are  op¬ 
erated  at  the  third  harmonic  frequency  of  10.17  MHz  for 
the  F/0.1  lens  and  10.7  MHz  for  the  FI  1  lens.  Many  prac¬ 
tical  considerations  go  into  this  particular  choice  of  oper¬ 
ating  frequency  range.  The  inversion  formalism  is  most 
conveniently  applied  to  a  lossless  liquid  medium.  The  loss 
in  water  is  essentially  negligible  at  10  MHz  and  so  satis¬ 
fies  this  requirement.  Moreover,  the  wavelength  in  water 
of  about  150  ^m  is  sufficiently  large  to  render  any  irreg¬ 
ularity  and  instability  in  the  vertical  translation  mecha¬ 
nism  insignificant.  By  the  same  token  the  measurement  is 
less  subject  to  error  as  a  result  of  surface  finish  imperfec¬ 
tions  in  the  sample. 

A  lead  sample  is  used  to  calibrate  the  pupil  function  of 
the  lens.  With  a  longitudinal  velocity  vL  —  2200  m/s,  the 
longitudinal  critical  angle  for  the  water-lead  interface  oc¬ 
curs  at  6  =  43°,  as  illustrated  in  the  theoretical  plot  in 
Fig.  2.  The  reflectance  function  is  essentially  flat  in  both 
amplitude  and  phase  from  9  =  0°  to  40°.  Fig.  3  shows  the 
magnitude  and  phase  of  the  V(z)  for  lead  from  z  =  -  8 
mm  to  z  =  +8  mm.  The  phase  displayed  is  in  fact  that  of 
F(j)  with  the  linear  phase  term  2 kz  subtracted  out.  Fig.  4 
shows  P2(0)  for  the  F/0.1  lens  obtained  by  inverting  the 
experimental  V(z)  for  lead.  The  pupil  function  is  not  uni¬ 
form  in  cither  magnitude  or  phase  especially  near  the  edge 
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Fig.  2.  Magnitude  and  phase  of  the  theoretical  reflectance  function  of  a  Fig.  5.  Calibration  of  the  pupil  function  of  the  10.7  MHz  F/\  focused  trans¬ 
water-lead  interface.  ducer. 


LEfl0  FUSED  SILICA 


MHz  FI 0.7  focused  transducer.  Fig.  6.  Experimental  F'(z)  of  water-fused-silica  interface  at  10.17  MHz. 


P*(0)  OF  F/07  LENS 


tf  ICf  2Cf  g  3  tf  40“  50* 


of  the  lens  where  an  accentuated  peak  occurs  around  6  = 
38°.  Premature  roll-off  of  the  lens  angular  response  is  also 
evident,  and  the  illumination  from  the  outer  periphery  of 
the  lens  is  weak.  The  experimental  P2(d)  for  the  FI  1  lens 


shown  in  Fig.  5  exhibits  a  similar,  though  less  pro¬ 
nounced,  peak  at  d  =  23.5°.  These  anomalies  in  the  lens 
response  are  due  to  nonuniformity  in  the  thickness  of  the 
matching  layer  as  a  result  of  difficulties  encountered  in  the 
fabrication  process  [23]. 

It  should  be  pointed  out  that  this  simple  calibration  pro¬ 
cedure  can  be  utilized  as  an  alternative  method  for  char¬ 
acterizing  and  predicting  imaging  performance  of  focused 
systems.  Instead  of  measuring  the  point  spread  function 
at  the  focal  plane  as  is  the  case  in  most  conventional  tech¬ 
niques,  the  characterization  is  done  in  terms  of  the  pupil 
function  illumination.  This  approach  has  direct  relevance 
and  application  in  the  fabrication  of  spherical  acoustic 
lenses,  where  the  sphericity  of  the  lens  and  the  uniformity 
of  the  matching  layer  are  of  primary  concern. 

Fused  Silica 

The  experimental  V(z)  for  fused  silica  is  shown  in  Fig. 
6.  Fig.  7  displays  the  corresponding  R(d),  which  is  ob¬ 
tained  by  inverting  the  K(z)  and  then  dividing  out  the  cal¬ 
ibrated  value  of  P2(d)  for  the  FI 0.7  lens.  The  theoretical 
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Fig.  7.  Comparison  of  the  theoretical  and  experimental  reflectance  function  Fig.  9.  Comparison  of  the  theoretical  and  experimental  reflectance  function 
for  a  water- fused-silica  interface.  for  a  water-aluminum  interface. 

ALUMINUM  SILICON  NITRIDE 


I  °K- 

?  i  A  , 


-L.il 


Fig.  8.  Experimental  F(z)  of  a  water-aluminum  interface  at  10.17  MHz.  Fig.  10.  Experimental  F(z)  of  a  water-silicon-nitride  interface  at  10.17 

MHz. 


reflectance  is  also  shown  for  comparison.  A  value  of 
vw  =  1486  m/s  is  used  as  the  velocity  in  water  in  the 
inversion  algorithm.  There  is  good  general  agreement  be¬ 
tween  theory  and  experiment  in  both  magnitude  and  phase. 
The  experimental  shear  critical  angle  at  9  =  23.5°  and  the 
Rayleigh  critical  angle  at  9  =  25.85°  compare  extremely 
well  with  the  theoretical  values.  The  longitudinal  critical 
angle  is  not  reproduced  in  the  experimental  result  due  to 
the  problem  of  degradation  angular  resolution  caused  by 
the  spatial  truncation  of  the  V(z)  data.  This  difficulty  will 
be  elaborated  upon  in  Section  V. 

Aluminum 

Figs.  8  and  9  are  the  V(z)  and  R(9)  for  aluminum.  Again 
the  agreement  between  theory  and  experiment  is  very 
good.  The  longitudinal  and  shear  critical  angles  show  up 
clearly  in  the  experimental  result,  in  this  case  at  9  =  13.64° 
and  9  =  28.64°,  respectively.  The  Rayleigh  critical  angle 
is  located  at  9  =  30.86°,  which  gives  a  phase  velocity  of 
2897  m/s  for  the  leaky  surface  acoustic  wave. 


Silicon  Nitride 

Figs.  10  and  11  show  the  V(z )  and  R(6)  for  hot-pressed 
silicon  nitride,  which  is  a  ceramic  material  with  high 
acoustic  velocity.  All  the  critical  angles  occur  at  low-in¬ 
cidence  angles,  and  only  the  Rayleigh  critical  angle  at 
6  =  15.0°  is  reproduced  in  the  experimental  phase  curve. 
It  should  be  noted  that  the  dip  in  the  experimental  mag¬ 
nitude  curve  for  silicon  nitride  and  smaller  dips  at  the  Ray¬ 
leigh  critical  angle  in  Figs.  7  and  9  for  fused  quartz  and 
aluminum  are  only  an  artifact  introduced  by  the  spatial 
truncation  of  the  experimental  V(z)  data.  Similar  dips  can 
occur  as  a  result  of  shear  loss  in  the  solid,  an  effect  widely 
known  as  the  Rayleigh  critical-angle  phenomenon  [14], 
[15].  However,  the  acoustic  losses  in  all  three  materials 
are  known  to  be  negligible  at  10  MHz,  therefore  the  ex¬ 
perimental  dips  cannot  possibly  be  associated  with  loss. 
As  will  be  shown  in  the  following  section  on  error  analy¬ 
sis,  this  experimental  artifact  is  well  understood  and  in 
some  cases  can  be  alleviated.  Hence,  the  inversion  tech- 
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SILICON  NITRIDE 


PLEXIGLAS 


Fig.  II.  Comparison  of  the  theoretical  and  experimental  reflectance  func¬ 
tion  for  the  water-silicon-nitride  interface. 


Fig.  13.  Comparison  of  the  theoretical  and  experimental  reflectance  func¬ 
tion  for  a  water-Plexiglas  interface. 


PLEXIGLAS 


TEFLON 


Fig.  12.  Experimental  F(z)  of  a  water-Plexiglas  interface  at  10.17  MHz. 


nique  does  not  preclude  the  characterization  of  such 
acoustic  properties  as  shear  loss,  temperature  dependence 
of  the  Rayleigh-wave  velocity,  and  various  contributing 
factors  to  the  Rayleigh  critical-angle  phenomenon. 

Plexiglas 

Figs.  12  and  13  are  the  V(z)  and  R{d),  respectively,  for 
plexiglas  or  lucite,  a  low-acoustic  velocity  and  high  loss 
material.  This  example  illustrates  how  material  loss  can 
be  determined.  The  longitudinal  critical  angle  is  located 
at  0  =  32.7°,  which  corresponds  to  a  longitudinal  wave 
velocity  of  2750  m/s.  The  magnitude  of  the  reflectance 
peaks  at  0.7  rather  than  one,  as  it  would  be  for  a  lossless 
substrate.  The  amount  of  diminution  of  the  peak  level  de¬ 
pends  on  the  loss  for  the  longitudinal  mode.  The  theoret¬ 
ical  curve  is  fitted  to  the  experimental  one  by  varying  the 
longitudinal  loss  factor  Qt.  The  QL  is  found  to  be  about 
50,  which  translates  into  an  attenuation  coefficient  of  232 
Np/m  or  20  dB/cm  at  10  MHz,  which  agrees  well  with 
other  published  values  [16]. 


Fig.  15.  Comparison  of  the  theoretical  and  experimental  reflectance  func¬ 
tion  for  a  water-teflon  interface. 

Teflon 

V(z)  and  R(6)  for  teflon  are  displayed  in  Figs.  14  and 
15.  As  expected  no  critical  angle  is  observed,  since  teflon 
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5/im  GOLD  FILM  ON 
FUSED  SILICA 


Fig.  16.  Experimental  K(z)  of  a  water  and  5-jim  gold  61m  on  a  fused-silica 
interface  at  10.7  MHz. 


S/lm  GOLD  FILM  ON 
FUSED  SILICA 


Fig.  17.  Comparison  of  the  theoretical  and  experimental  reflectance  func¬ 
tion  for  a  water  and  5-fim  gold  film  on  a  fused-silica  interface. 


has  a  longitudinal  velocity  of  1400  m/s,  which  is  lower 
than  that  of  water. 

Thin-Film  Measurement 

The  following  example  demonstrates  the  potential  ap¬ 
plication  of  acoustic  microscopy  in  thin-film  character¬ 
ization  and  also  confirms  a  thin-film  matching  phenome¬ 
non  predicted  by  numerical  computation.  The  thin-film 
structure  used  in  this  example  is  a  5-jim-thick  gold  film 
deposited  on  a  fused-silica  substrate.  At  10  MHz  the  thick¬ 
ness  of  the  gold  film  corresponds  to  about  1.5  percent  of 
the  longitudinal  wavelength.  Figs.  16  and  17  are  the  plots 
of  V(z)  and  R(d),  respectively.  The  agreement  with  the 
theoretically  generated  R(d)  is  good.  R(8)  shows  a  null  at 
6  =  17.23°,  an  angle  between  the  longitudinal  and  shear 
critical  angles  of  the  fused  silica  substrate.  Physically  this 
means  that  the  incident  longitudinal  mode  in  water  cou¬ 
ples  very  strongly  into  a  bulk  propagating  mode  in  the 
fused  silica  substrate.  The  physics  of  this  coupling  is  not 
well  understood,  but  it  is  believed  that  the  longitudinal 


wave  in  water  converts  in  the  gold  film  to  a  leaky  Sezawa 
wave,  which  leaks  into  the  fused  silica  substrate  in  the 
form  of  a  propagating  shear  wave.  The  coupling  efficiency 
and  the  angle  of  incidence,  at  which  maximum  transmis¬ 
sion  occurs,  have  been  shown  to  be  functions  of  the  film 
thickness  by  numerical  analysis.  Thus  by  measuring  the 
reflectance,  one  can  get  information  about  the  film  thick¬ 
ness.  The  experimental  demonstration  of  this  thin-film 
phenomenon  also  opens  up  the  possibility  of  using  the  thin- 
film  matching  technique  in  the  design  of  acoustic  trans¬ 
ducers  to  improve  the  transmission  efficiency  between  two 
media  with  vastly  different  acoustic  impedances. 


V.  Error  Analysis 

The  introduction  of  experimental  artifacts  in  the  exper¬ 
imental  R(d)  due  to  the  spatial  truncation  of  the  V(z)  curve 
is  treated  in  this  section.  This  problem  arises  because  of 
the  finite  distance  over  which  V(z)  data  can  be  collected. 
Equation  (26)  shows  that  V(u)  and  [P2(f)/?(t)]  form  a  Fou- 
rier-transform  pair.  It  can  be  shown  by  Fourier-transform 
theory  that  since  [P2(f) /?(t)]  is  finite  in  the  t  domain,  V(u) 
has  to  be  infinite  in  the  u  domain.  Experimentally  V(z) 
can  only  be  obtained  for  some  finite-width  data  window. 
Therefore  the  actual  waveform  V'(z)  used  in  the  inversion 
is  a  truncated  version  of  V(u).  Thus 

V'(u )  =  V(u)  rect  (u/D)  (32) 

where  u  =  z/X,  and  D  is  the  width  of  the  data  window 
defined  in  terms  of  the  number  of  wavelengths  in  water. 
Noting  that  multiplication  in  the  u  domain  corresponds  to 
convolution  in  the  transform  domain  t,  the  reflectance 
function  obtained  by  inverting  V '(u)  is 

[P2(t)P(t)]'  =  [P\t)  P(f)]*  D  sine  (Dt)  (33) 

where  *  denotes  convolution. 

The  effects  of  the  convolution  between  the  reflectance 
function  and  the  waveform  sine  (Dt)  are  twofold.  First,  the 
angular  resolution  in  the  t  domain  is  degraded.  Second, 
because  of  the  oscillatory  nature  of  sine  (Dt),  sharp  fea¬ 
tures  in  P2(t)R(t)  tend  to  generate  ripples  in  the  resulting 
[P2(t)/?(f)]'.  The  resolution  degradation  effect  can  be  es¬ 
timated  as  follows.  The  full  width  between  zeroes  of  the 
main  lobe  of  sine  (Dt)  is 

At  =  2  ID.  (34) 

The  At  is  essentially  the  transition  width  of  the  response 
to  a  sharp  step  in  P2(t)  R(t)  and  hence  can  be  regarded  as 
the  worst  case  resolution  in  t.  Using  the  relation  t  =  2  cos 
6  given  in  (24),  we  get  the  following  expression  for  the 
angular  resolution  in  terms  of  the  angle  of  incidence  6 


A6 


At 

2  sin  6 


(35) 


Combining  (34)  and  (35),  the  resolution  in  0  in  terms  of 
the  width  of  the  data  window  D  is  found  to  be 


A8 


_ 1_ 

D  sin  6 


(36) 
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Fig.  18.  Variation  of  angular  resolution  degradation  as  a  function  of  the 
incident  angle  and  the  normalized  data  window  width  D. 


SILICON  NITRIDE 


Fig.  19.  Simulated  reflectance  function  of  a  water-silicon  nitride  interface 
obtained  by  inverting  spatially  truncation  F(z)  data. 


Equation  (36)  is  plotted  in  Fig.  18  with  D  as  a  param¬ 
eter.  The  data  window  D  used  in  the  experiment  is  about 
100.  This  implies  angular  resolution  of  about  3.3°  at  0  = 
10°,  1.7°  at  0  =  20°,  and  1.1°  at  0  =  30°.  The  experi¬ 
mental  reflectance  of  silicon  nitride  in  Fig.  11  is  a  patho¬ 
logical  example  of  the  degradation  effect  of  this  angular 
resolution.  The  Rayleigh  critical  angle  occurs  at  a  low 
value  of  0  =  15.0°  and  in  the  vicinity  of  this  angle,  the 
magnitude  of  the  reflectance  is  unity  but  the  phase  goes 
through  a  rapid  2 r  radian  change  over  a  2°  angular  range. 
The  net  result  of  the  convolution  with  sine  (Df)  is  a  sharp 
dip  in  the  magnitude  of  the  experimental  reflectance  and 
a  smoothing  of  the  phase  curve  near  the  Rayleigh  critical 
angle.  The  effect  of  the  convolution  is  much  less  serious 
at  angles  of  high  incidence  as  shown  in  Fig.  18  and  as 
evidenced  by  the  experimental  R(0)’s  of  fused  silica  and 
aluminum,  where  only  a  slight  dip  occurs  in  the  reflec¬ 
tance  magnitude. 

A  computer  simulation  has  been  carried  out  where  the 
theoretical  R(0)  for  silicon  nitride  is  used  to  generate  the 
V{z),  which  is  truncated  and  then  inverted  in  an  identical 
manner  to  the  experimental  data.  The  result  is  shown  in 
Fig.  19.  The  simulated  R(0)  exhibits  exactly  the  same  be¬ 
havior  in  both  amplitude  and  phase  as  the  experimental 
one  in  Fig.  11. 

The  locations  of  the  critical  angles  are  generally  used  to 
determine  the  phase  velocities  of  propagating  modes, 
which  are  given  by 


o.  =  iv/sin  0 

where  vw  is  the  velocity  in  water;  the  asterisk  stands  for 
L  (longitudinal),  S  (shear),  or  R  (Rayleigh);  and  the  sub¬ 
script  c  denotes  critical  angle.  The  percentage  error  in  the 
estimation  of  v.r  as  a  consequence  of  the  truncation  of 
V(z)  can  be  shown  to  be 


Av. 

v. 


1  cos  0.  „  1  cos  0 

- c-  Ad.  = - 5—. 

2  sin  0.r  ‘ID  sin2  6 


(37) 


Equation  (37)  is  plotted  in  Fig.  20  with  D  as  a  parameter. 


Fig.  20.  Phase  velocity  estimation  error  as  a  function  of  angle  of  incidence 
due  to  degradation  of  the  angular  resolution. 


The  error  decreases  drastically  with  increasing  incidence 
and  increasing  D.  For£>  =  100,  |  AvJv.\  is  four  percent 
at  0  =  20°  and  1.75  percent  at  0  =  30°.  One  should  bear 
in  mind  that  (37)  represents  the  worst-case  estimation,  and 
the  actual  error  could  well  be  substantially  smaller,  espe¬ 
cially  for  the  Rayleigh  critical  angle  0Rc.  The  determina¬ 
tion  of  6Rc  involves  locating  the  point,  where  the  phase  of 
R(6)  is  t  radians.  Around  0Rc  the  magnitude  of  the  reflec¬ 
tance  function  for  a  lossless  material  is  constant  while  the 
phase  ¥  can  be  shown  to  have  the  form  [9] 


tan  ¥  =  2 


K 


k  sin  dRc 
a 


(38) 


where  a  is  the  leak  rate  of  the  Rayleigh  wave.  The  phase 
function  ¥  is  the  antisymmetric  in  kr  about  'If  =  ir  around 
the  point  kr  =  k  sin  $Rc.  It  is  similarly  antisymmetric  in  k. 
or  t  in  the  neighborhood  of  kz  —  k  cos  dRc.  Since  the  sine 
function  with  which  [P2(/) /?(f)l  is  convolved  is  symmet¬ 
ric,  the  phase  function  Sk  around  the  Rayleigh  critical  an¬ 
gle  remains  essentially  unchanged  by  the  convolution. 
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provided  that  the  width  of  the  main  lobe  of  the  sine  func¬ 
tion  At  is  smaller  than  the  transition  width  of  ¥  through 
the  Rayleigh  critical-angle  region.  For  example,  in  the  case 
of  fused  silica  and  aluminum,  the  respective  leaky  sur¬ 
face-acoustic-wave  velocities  calculated  from  the  experi¬ 
mental  phase  curve  are  almost  identical  to  those  predicted 
by  theory. 

The  second  problem  associated  with  the  generation  of 
ripples  can  be  partly  remedied,  though  with  some  further 
loss  in  angular  resolution,  by  applying  a  smooth  apodi- 
zation  function  for  the  V(z)  data  before  inversion  takes 
place.  The  apodization  function  used  in  processing  the 
V(z)  data  is 


a(u)  = 


sin2  ( itu/D ) 
n  +  sin2  ( ttu/D ) 


(38) 


where  n  is  a  free  parameter  for  adjusting  the  tapering 
characteristics  of  the  apodization.  The  n  =  0.1  is  used  to 
process  the  experimental  data  shown  here.  The  resulting 
apodization  essentially  leaves  the  data  in  the  center  of  the 
window  unchanged,  but  it  behaves  much  like  Hamming 
weighting  at  the  edges  of  the  data  window. 

The  obvious  solution  to  both  the  resolution  degradation 
and  ripple  problems  is  to  increase  the  data  window  width 
D  of  V'(u).  Equation  (31)  shows  that  this  can  be  realized 
by  increasing  the  radius  of  curvature  /0  or  reducing  the 
aperture  size  or  both.  Also  D  can  be  increased  by  increas¬ 
ing  the  frequency  of  operation.  For  silicon  nitride,  R(6) 
can  be  reproduced  much  more  faithfully  if  A0  is  reduced 
to  j",  which  can  be  achieved  by  using  an  FI  1.5  lens  with 
a  focal  length  of  32  mm  operating  at  20  MHz. 

Another  important  source  of  error  is  in  the  estimation 
of  the  velocity  of  water.  Since 

u  =  z/\  =  zf/v 


then 


du  =  dv  =  -udv/v. 


(39) 


Suppose  the  wrong  velocity  is  used  in  the  inversion.  In 
(25)  this  is  equivalent  to  changing  u  to 

u'  =  u  +  Au  =  u(l  -  Aviv).  (40) 

The  resulting  inversion  is  given  by 

[/*2(0 /?(*)]'  =  |  V(u)  exp  [j2iru(l  -  Avlv)t]  du 
or 


that  the  shift  distortion  as  a  function  of  incidence  angle  is 
of  the  form 


|  Ad  |  =  cot  6  Aviv  (42) 

which  shows  that  the  most  serious  errors  occur  at  low  an¬ 
gles  of  incidence. 

The  problem  of  determining  the  true  transform 
P2(d)R(0)  from  a  finite  segment  V'(z )  of  V(z)  is  a  common 
one  in  Fourier  analysis  and  spectral  estimation.  Various 
techniques  exist  in  the  literature  for  extrapolating  V'(z)  so 
that  a  more  accurate  determination  of  P2{d)R(B)  can  be 
made.  Since  P2(0)f?(0)  is  a  bandlimited  function,  the 
maximum  spatial  frequency  being  limited  by  the  angular 
extent  of  the  pupil  function,  V(z)  is  analytic  in  the  entire 
z  axis  [17].  In  principle,  an  iterative  algorithm  proposed 
by  Papoulis  [18]  can  be  used  to  improve  the  accuracy  of 
the  estimation  of  P1(d)R(6). 

Conclusion 

We  have  demonstrated  that  the  reflectance  function  of  a 
liquid-solid  interface  can  be  determined  by  inverting  the 
corresponding  complex  V(z)  data  from  an  acoustic  micro¬ 
scope.  This  inversion  technique  represents  a  more  com¬ 
plete  approach  to  material  characterization  than  previous 
F(z)-related  work.  The  phase  velocities  of  the  various 
propagating  modes  in  the  solid  medium  can  be  obtained 
directly  from  the  reflectance  function.  The  effect  of  ma¬ 
terial  loss  can  also  be  observed  and  quantified.  Moreover, 
this  measurement  technique  provides  a  means  of  gauging 
imaging  performance  of  focused  systems  by  directly  mea¬ 
suring  the  pupil  function  illumination.  In  addition,  there 
are  useful  practical  applications  in  the  area  of  thin-film 
characterization,  and  many  interesting  possibilities  exist 
for  more  complex  structures,  such  as  multi-layered  films. 
The  nonparaxial  formulation  of  the  V(z)  integral  is  im¬ 
portant  in  that  it  lays  a  sound  theoretical  foundation  for 
the  inversion  measurement  technique.  The  excellent 
agreement  between  the  theoretically  and  experimentally 
obtained  reflectance  functions  further  supports  the  validity 
of  the  nonparaxial  theory. 

Although  this  work  has  been  carried  out  and  discussed 
in  the  context  of  acoustic  microscopy,  the  validity  and  ap¬ 
plicability  of  many  of  the  underlying  concepts  extend  to 
optical  microscopy  as  well.  Provided  one  can  accurately 
measure  the  optical  phase,  which  is  not  trivial  but  cer¬ 
tainly  realizable  [22],  the  inversion  algorithm  described 
here  can  be  used  to  obtain  the  optical  reflectance  function 
from  the  corresponding  complex  optical  P(z)  function. 


[/>2(/)/?(r)]'  =  P2fl[(l  -  Au/u)r].  (41) 

Therefore  using  an  erroneous  v  results  in  a  stretched  re¬ 
flectance  function  in  t.  Since  the  velocity  of  water  has  a 
large  temperature  coefficient,  4  m/s  per  °C,  this  error  can 
be  significant.  It  may  cause  enough  misalignment  of  the 
inversions  of  the  sample  material  of  interest  and  the  cali¬ 
bration  lead  sample  to  have  a  serious  effect  on  the  deter¬ 
mination  of  R(fi).  From  (35)  and  (41)  it  can  easily  be  shown 
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Appendix 

Auld  and  Kino  [  19]— [21  ] ,  by  using  the  reciprocity  theo¬ 
rem,  were  able  to  determine  the  normalized  reflected  sig¬ 
nal  or  reflection  coefficient  itl  from  an  object.  The  theory 
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where  the  integral  is  taken  over  the  surface  of  the  object, 
u  is  the  displacement,  and  p  is  the  pressure  fields  associ¬ 
ated  with  an  incident  wave  of  temporal  frequency  to.  Here 
exp  {jwt)  time  dependence  is  assumed.  The  superscript  i 
denotes  the  incident  or  transmitted  wave  when  the  object 
is  not  present,  and  the  unsuperscripted  terms  denote  the 
total  field  at  the  obstacle.  The  parameter  P  is  the  power 
exciting  the  transducer  for  a  given  incident  signal  u‘,  p‘. 

The  total  fields  can  be  written  in  the  form 


u  =  u‘  +  ur  (A2) 

P  =  P‘  +  Pr  (A3) 

where  the  superscript  r  denotes  the  waves  reflected  from 
the  object.  Substituting  (A2)  and  (A3)  into  (Al)  yields 


(u‘pr  -  urp‘)n  dS 


(A4) 


If  the  object  is  a  semi-infinite  plane  normal  to  the  z  direc¬ 
tion,  s,,  can  be  written  in  the  normalized  form 


\  {u[pr  -  u'p)dS 

nt)  = - = - - -  (A5) 

2  J  i/'p'  dS 

where  the  asterisk  denotes  the  complex  conjugate.  Finally 
it  is  convenient  to  write  the  pressure  in  terms  of  the  po¬ 
tential.  For  a  liquid  with  u  =  V</>,  it  can  be  shown  that 
p  =  jwp<t>.  Hence  it  follows  that 


{*(«;*'  -  «;<*>')  ds 

V{Z)  = - - - -  (A6) 

2  j  dS 

Note  that  for  a  perfect  plane  reflector  located  at  the  focal 
plane  z  =  0  of  a  lens,  <j>r  =  <j>‘  -  <t> and  urz  =  -u[; 
therefore  V(z)  =  I. 
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A  new  phase  measurement  system  based  on  a  sychronous  detection  scheme  is  described.  It  is 
capable  of  yielding  high  precision  phase  data  with  very  short  tone  bursts.  A  theoretical  analysis  of 
the  system  is  given  and  system  performance  characteristics  are  also  discussed.  This  phase 
measurement  technique  is  especially  suited  for  acoustic  microscopy  applications.  Examples  in 
material  characterization  and  surface  topography  mapping  are  presented. 


INTRODUCTION 

The  scanning  acoustic  microscope  is  a  high  resolution  coher¬ 
ent  imaging  system  in  which,  unlike  conventional  optical 
microscopes,  both  the  amplitude  and  phase  of  the  output 
signal  are  available.  With  only  a  few  exceptions,  attention 
has  centered  only  on  the  use  of  amplitude  or  intensity  infor¬ 
mation.  Little  effort  has  been  made  to  extract  the  phase  of 
the  return  signals  separately,  despite  the  fact  that  it  has  been 
shown  that  the  independent  measurement  of  phase  and  am¬ 
plitude  can  be  very  useful.1'3  The  reason  for  the  current  re¬ 
luctance  to  make  use  of  the  phase  is  that  short  tone  burst 
signals  are  generally  used  in  acoustic  microscopy  to  enable 
time  differentiation  among  multiple  echoes.  It  is  not  general¬ 
ly  trivial  to  electronically  measure  the  phase  of  a  high-fre¬ 
quency  tone  burst  with  sufficient  accuracy  to  be  viseful. 

We  will  describe  here  an  approach  which  is  relatively 
easy  to  implement  and  which  is  capable  of  yielding  high- 
precision  phase  data,  even  from  very  short  tone  bursts.  A 
Fourier  analysis  of  the  measurement  system  is  carried  out  to 
elucidate  the  underlying  principle  of  the  measurement  tech¬ 
nique.  The  phase  sensitivity  of  the  system  in  the  presence  of 
noise  is  also  analyzed. 

We  will  show  a  number  of  practical  applications  in 
which  phase  is  used  exclusively  to  carry  out  measurements 
which  otherwise  would  either  be  impossible  or  less  desirable 
to  realize  with  amplitude-only  techniques.  The  examples 
given  involve  the  measurement  of  Rayleigh  wave  velocity 
perturbation  for  film  thickness  determination,  and  the  use  of 
the  acoustic  microscope  as  a  noncontacting  profilometer. 
We  will  also  demonstrate  how  the  phase  velocities  and  at¬ 
tenuation  of  propagating  modes  in  a  solid  material  can  be 
measured  by  using  the  complex,  i.e.,  both  amplitude  and 
phase,  output  signal  of  an  acoustic  microscope. 

I.  USE  OF  PHASE  IN  ACOUSTIC  MICROSCOPY 

Phase  measurements  in  acoustic  microscopy  can  give  a 
direct  indication  of  material  properties,  such  as  the  Rayleigh 
wave  velocity,  and  can  be  used  to  infer  the  width  and  the 
height  of  surface  features. 

In  making  precise  phase  measurements  which  are  relat¬ 
ed  to  physical  properties,  it  is  important  to  ensure  that  the 
reference  signal  against  which  the  phase  of  the  probing  sig¬ 


nal  is  to  be  compared  and  the  probe  itself  share  as  many  of 
the  instrumental  and  environmental  phase  disturbances  as 
possible.  In  other  words,  precise  measurements  are  best  done 
in  an  interferometer  in  which  the  two  arms  are  closely 
matched. 

In  this  work  we  have  made  use  of  two  acoustic  operating 
configurations  where  this  general  condition  is  largely  satis¬ 
fied.  The  first  is  where  the  acoustic  lens  is  defocused  so  that 
the  sample  is  located  above  the  focal  plane,  as  shown  in  Fig. 
1(a).  The  two  components  of  the  acoustic  field  returned  to 
the  lens  which  contribute  most  significantly  to  the  output 
signal  are  ( a )  the  on-axis  specular  reflection  of  the  longitudi¬ 
nal  wave  in  the  water  [L  in  Fig.  1(a)],  and  (b)  the  off-axis 
rays  (R)  associated  with  the  leaky  Rayleigh  wave  on  the 
solid.4,5  These  rays  satisfy  the  condition  vw/vK  =  sin  0R, 
where  vw  is  the  longitudinal  wave  velocity  in  water  and  vR  is 
the  leaky  Rayleigh  wave  velocity. 


RAY  MODEL 


Fig.  1.  Acoustic  micro¬ 
scope  operating  configura¬ 
tions  for  (a)  Rayleigh  wave 
velocity  perturbation  mea¬ 
surement,  (b)  topography 
mapping. 


ACOUSTIC  MICROSCOPE  LENS 
FOR  T0P06RAPHY  PROFILING 


446  Rev.  Set.  Instrum.  57  (3),  March  IMS 


0034-6748/86/030446-07301 .30  ©  1 986  American  Institute  ot  Physics 


446 


■WHMW  IMWlUWUlftJ  WU  WUTUIRJ  iffimj  mi  w  iwmi  tm  w  vpu  rjmi  m  w  icj  »u  *  u  s  trwv  »v  wj 


These  two  signals  L  and  R  differ  slightly  in  path  length 
in  the  water,  but  otherwise  experience  much  the  same  envi¬ 
ronmental  disturbances.  However,  the  Rayleigh  wave  com¬ 
ponent  has  substantial  phase  delay  over  the  surface  of  the 
solid  material.  Hence,  the  effective  path  length  for  the  R 
signal  is  longer  than  that  for  the  specular  reflection  L.  In¬ 
deed,  in  our  experiments,  we  use  the  different  arrival  times  of 
these  two  return  signals  to  facilitate  their  separation  and  sub¬ 
sequent  phase  comparison.  The  exciting  signal  is  two  to 
three  cycles  in  duration,  and  the  defocus  distance  is  suffi¬ 
cient  so  that  there  is  no  temporal  overlap  of  the  signals.  It  is 
clear  that  changes  in  the  Rayleigh  wave  velocity  as  a  func¬ 
tion  of  position  along  the  surface  of  the  specimen  may  be 
sensed  in  this  way.  By  combining  a  high-accuracy  phase 
measurement  scheme  with  a  suitably  extended  Rayleigh 
wave  path,  remarkably  sensitive  measurements  can  be  made. 

The  second  operating  configuration  is  illustrated  in  Fig. 
1  (b).  In  this  case  the  lens  is  positioned  so  that  its  focus  is  at 
the  surface  of  the  specimen,  or  slightly  above  it,  so  that  no 
Rayleigh  wave  of  importance  is  excited.  A  reference  path  is 
provided  by  an  annular  beam  that  propagates  through  the 
flat  outer  periphery  of  the  lens.  By  exciting  the  lens  with  a 
short  pulse,  and  once  again  using  time  discrimination  to  sep¬ 
arate  the  signals  from  the  two  different  paths,  an  interferom¬ 
eter  is  formed.  Small  local  changes  in  the  surface  topography 
result  in  large  changes  in  the  phase  of  the  focused  beam  rela¬ 
tive  to  the  phase  of  the  large  diameter  reference  beam. 

In  both  measurement  configurations,  it  is  apparent  that 
the  reference  and  signal  pulses  arrive  at  different  times.  This 
would  ordinarily  make  phase  comparison  impossible;  how¬ 
ever,  the  phase  measurement  system  described  in  the  next 
section  essentially  reconstructs  two  phase-coherent,  contin¬ 
uous  signals  which  are  directly  related  to  the  amplitudes  and 
phases  of  the  two  time-distinct  return  signals. 

The  first  configuration  is  also  used  where  the  signal  of 
interest  is  the  transducer  output  voltage  when  the  system  is 
being  operated  in  the  so-called  V(z)  mode.6  In  this  case,  the 
exciting  signal  has  a  much  longer  duration  so  that  the  return 
signals  (L  and  R)  overlap  and  interfere  at  the  transducer, 
even  for  large  defocus  distances.  The  idea  here  is  to  measure 
both  the  amplitude  and  the  phase  of  the  signal  relative  to 
some  reference  as  a  function  of  the  axial  separation  z  between 
the  acoustic  lens  and  the  specimen.  The  resulting  complex 
V(z)  can  be  inverted7,8  to  obtain  the  angular  dependence  of 
the  reflectance  of  the  liquid/specimen  interface.  From  the 
reflectance  function,  various  material  property  parameters 
can  be  extracted. 

In  most  of  the  experiments  described  here,  a  center  fre¬ 
quency  of  SO  MHz  was  used  with  an  acoustic  transducer 
having  20%  bandwidth.  The  lens  material  was  fused  quartz 
and  had  a  radius  of  curvature  of  3.2  mm,  giving  a  focal  length 
in  water  of  4.24  mm.  The  opening  aperture  of  the  lens  was  5.0 
mm  in  diameter,  corresponding  to  a  maximum  half  angle  of 
36*  or  an/ number  of  0.85. 

II.  PHASE  MEASUREMENT  SYSTEM 

The  phase  of  a  short  rf  pulse  is  difficult  to  measure  di¬ 
rectly  because  the  corresponding  Fourier  spectrum  is  a  con¬ 
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Fig.  2.  Fourier  spectra  of  a  single  rf  pulse  and  a  synchronous  pulse  train. 

tinuous  spread  of  spectral  components,  as  shown  in  Fig. 
2(a).  Only  a  minute  amount  of  power  is  associated  with  the 
frequency  co0  of  the  rf  carrier  from  which  the  pulse  is  derived. 
An  acoustic  microscope,  however,  is  excited  repetitively  at  a 
constant  rate  and  the  rf  pulses  occur  as  a  pulse  train.  If  the 
pulses  are  coherent  in  the  sense  that  they  originate  from  a 
single  cw  signal,  as  depicted  in  Fig.  2(b),  the  time  waveform 
can  be  regarded  as  a  periodically  sampled  cw  signal.  The 
corresponding  Fourier  spectrum  consists  of  discrete  fre¬ 
quency  components  spaced  by  the  repetition  rate  cos .  The 
frequency  component  situated  at  <y0  can  easily  be  extracted 
by  narrow-band  filtering,  as  shown.  In  essence,  a  cw  replica 
of  the  pulse-moduated  rf  signal  is  generated  to  facilitate 
phase  measurement. 

Figure  3  shows  the  configuration  of  the  phase  measure¬ 
ment  system  used  to  carry  out  precision  phase  measurement 
with  a  50-MHz  acoustic  microscope.  The  basic  scheme  in¬ 
volves  signal  recovery  by  synchronous  detection  and  subse¬ 
quent  phase  measurement  with  a  lock-in  amplifier. 

The  signal  source  is  a  single  sideband  (SSB)  generator 
with  synchronous  cw  outputs  at  100  kHz,  10.6,  and  10.7 

50  MHt  PWR  PVN 


Fig.  3.  Schematic  for  the  phase  measurement  system- 
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MHz.  The  10.7-MHz  output  is  derived  from  the  upper  side¬ 
band  component  of  the  product  of  the  100  kHz  and  10.6 
MHz  signals.  The  companion  lower  sideband  at  10.5  MHz  is 
suppressed  by  at  least  50  dB  using  a  standard  FM  radio  sys¬ 
tem  IF  filter.  The  translating  oscillator  shifts  the  operating 
frequency  up  to  the  center  frequency  of  the  acoustic  trans¬ 
ducer,  which  is  50  MHz  in  this  example.  A  49.9-MHz  refer¬ 
ence  signal  is  also  generated.  Again  the  lower  sideband  com¬ 
ponents  of  the  product  signals  have  to  be  removed.  However, 
the  filtering  requirement  is  not  as  stringent  as  in  the  SSB 
generator  and  tunable  bandpass  filters  ( F3  and  F4 )  with  5  % 
bandwidth  are  adequate  since  the  sideband  components  are 
well  separated  in  this  case.  This  signal  generation  scheme  is 
flexible  in  that  the  operating  frequency  is  tunable  over  a  fair¬ 
ly  wide  range,  up  to  100  MHz  .  For  a  higher  frequency  of 
operation,  F3  and  F4  can  be  replaced  with  bandpass  filters 
with  better  selectivity.  An  alternative  approach  is  to  add 
successive  heterodyning  and  filtering  stages  to  step  up  incre¬ 
mentally  to  the  desired  operation  frequency. 

The  50- MHz  cw  signal  is  time  gated  (SW4)  to  produce 
a  short  tone  burst  to  excite  the  acoustic  transducer.  Switch 
SW3  is  used  as  a  time  gate  to  pass  only  the  low-level  acoustic 
return  signals  of  interest,  and  to  block  high-level  extraneous 
ones  which  might  damage  the  preamplifier.  The  gated 
acoustic  return  signal  g(t)  typically  consists  of  two  time- 
separated  rf  pulses  g,  (/)  and  g2  (/),  whose  relative  phases 
are  to  be  measured.  For  example,  in  the  case  of  surface 
acoustic  wave  velocity  perturbation  measurements,  the  rela¬ 
tive  phase  between  the  on-axis  return  pulse  and  the  convert¬ 
ed  leaky  Rayleigh  wave  pulse  is  measured.  Alternatively,  for 
topography  measurements,  the  relative  phase  between  the 
focused  beam  signal  and  the  unfocused  ring  beam  signal 
gives  the  local  height  of  the  sample.  In  either  case,  the  non- 
overlapping  pulses  are  electrically  separated  through  time 
gating  by  SW1  and  SW2  into  two  channels,  as  shown.  Each 
of  them  is  mixed  with  the  49.9  MHz  reference,  and  the  pro¬ 
duct  is  narrowband  filtered  (FI  and  F2)  to  extract  the  100- 
kHz  component. 

It  will  be  shown  in  Sec.  IV  that  the  100-kHz  signals 
q[  (f)  and  q2  (r)  are  essentially  low-frequency  cw  replicas  of 
the  pulse-modulated  rf  acoustic  signals  g,  (f )  and  g2  (/), 
bearing  identical  amplitude  and  phase  information.  With  a 
bi-phase  lock-in  amplifier  tuned  to  100  kHz,  we  can  measure 
the  phase  difference  between  the  two  signals,  and  also  the 
amplitude  of  the  signal  fed  to  the  signal  channel  of  the  lock- 
in,  if  it  is  required.  The  phase  sensitivity  of  the  system  is 
essentially  limited  by  the  lock-in  amplifier  since,  in  most 
cases,  the  problem  of  system  noise  can  be  overcome  by  in¬ 
creasing  the  integration  time  in  the  lock-in  output  stage,  thus 
effectively  reducing  the  noise  bandwidth,  with  the  obvious 
disadvantage  of  longer  data  acquisition  time.  Good  quality 
lock-in  amplifiers  routinely  have  phase  resolutions  of  0. 1°, 
which  translates  into  an  overall  system  phase  sensitivity  of 
1/3600  of  a  wavelength. 

III.  APPLICATIONS 

A.  Velocity  perturbation  due  to  thln-film  overlay 

The  operating  configuration  shown  in  Fig.  1  (a)  is  used 
to  determine  the  change  in  Rayleigh  wave  propagation  ve¬ 


locity  as  a  result  of  material  property  change.  The  sample  in 
this  example  is  a  multiple-thickness  indium  film  deposited 
on  glass.  The  thicknesses  are  240  and  620  A,  respectively. 
The  objective  was  to  measure  the  perturbation  of  the  Ray¬ 
leigh  wave  velocity  caused  by  the  indium  film.  A  line  scan 
over  the  surface  of  the  sample  (Fig.  4)  exhibits  phase 
changes  of  7°  and  1 1°  for  the  240-  and  380-A  step  changes  in 
film  thickness.  The  spatial  resolution  of  the  system  is  defined 
by  the  Rayleigh  wave  path  length  on  the  substrate  and  is 
determined  to  be  about  0.8  mm  from  the  step  transition 
widths  in  the  line  scan.  It  can  be  calculated  from  first-order 
perturbation  theory9  that  the  velocity  perturbation  due  to 
the  indium  film  is  0.18%  and  0.46%  ,  respectively,  for  the 
240-  and  620- A  layers.  Based  on  these  estimated  parameters, 
one  would  expect  phase  changes  of  7.6°  and  12°  for  the  240- 
and  380-A  step  transitions.  Hence,  there  is  fairly  good  agree¬ 
ment  between  the  experimentally  obtained  and  theoretically 
predicted  phase  changes.  It  should  also  be  noted  that  in  the 
line  scan,  the  small  phase  variations  in  the  supposedly  flat 
regions  of  the  indium  film  are  real  and  repeatable.  The  fluc¬ 
tuations  are  less  than  0.5°  and  are  due  to  nonuniformity  in 
the  thickness  of  the  indium  film.  Since  the  phase  sensitivity 
of  the  system  is  limited  by  the  lock-in  amplifier  to  0.1°,  this 
measurement  technique  can  potentially  detect  velocity  per¬ 
turbations  on  the  order  of  one  part  in  105 .  Conversely,  if  the 
material  property  parameters  are  known,  we  can  determine 
film  thicknesses  accurately  from  the  velocity  perturbation 
measurements. 

In  this  operating  configuration,  the  measured  phase  is 
dependent  on  the  distance  between  the  lens  and  the  speci¬ 
men.  This  distance  changes  with  the  sample  surface  topo¬ 
graphy  and  is  susceptible  to  thermal  drift  as  well,  thus  intro¬ 
ducing  phase  error  into  the  measurement.  To  minimize  this 
error,  a  feedback  mechanism  is  employed  to  keep  the  lens-to- 
specimen  spacing  constant.  The  acoustic  lens  is  mounted  on 
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Fig.  4.  Measurement  of  velocity  perturbation  due  to  indium  film  on  glass. 
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a  piezoelectric  (PZT)  stack  so  that  its  vertical  position  can 
be  adjusted  continuously  by  an  electronic  control  signal.  The 
acoustic  on-axis  reflection  pulse,  whose  phase  is  a  direct 
measure  of  the  lens-to-sample  distance,  is  applied  to  the  ref¬ 
erence  channel  of  the  lock-in  amplifier,  as  shown  in  Fig.  3. 
The  reference  channel  generates  a  100-kHz  constant-ampli¬ 
tude,  phase-locked  replica  of  the  reference  input,  which  is 
then  compared  with  the  100-kHz  output  of  the  SSB  gener¬ 
ator  to  produce  the  control  signal  for  the  PZT  stack.  In  oper¬ 
ation,  the  acoustic  lens  automatically  tracks  the  surface  to¬ 
pography  of  the  sample  during  scanning  to  ensure  that  the 
measured  phase  change  is  due  to  material  property  variation 
alone. 

B.  Surface  topography  measurements 

We  have  used  the  acoustic  microscope  in  the  measure¬ 
ment  configuration  shown  in  Fig.  1  (b)  as  a  high-resolution 
noncontacting  profilometer. 10,11  Topography  images  of  me¬ 
tallized  stripe  patterns  on  a  fused  silica  substrate  are  shown 
in  Figs.  5(a),  5(b),  and  5(d).  The  patterns  have  progres¬ 
sively  finer  pitches  and  the  linewidths  are  250,  1 25,  and  62.5 
pm,  respectively.  The  metallization  is  gold  with  a  thickness 
of  about  3000  A.  The  gold  stripes  show  up  as  bright  areas  in 
the  images.  At  50  MHz,  with  an  F  #  of  0. 85  and  a  uniformly 
excited  aperture,  the  Rayleigh  resolution  of  the  acoustic  lens 
is  1.13F/1  =  0.96/1  or  29  pm  ,  while  the  3-dB  resolution  is 
0.64F/1  =  0.54/1  or  16  pm.  The  62.5-/zm  line  image  in  Fig. 
5(d)  is  clearly  resolved,  as  would  be  expected.  The  ampli¬ 
tude  image  of  the  125-/zm  stripes  is  shown  in  Fig.  5(c)  for 
comparison.  There  is  virtually  no  contrast  ( less  than  1  % )  to 
suggest  the  presence  of  a  striped  pattern,  clearly  illustrating 
that  phase  is  a  far  more  sensitive  means  of  gauging  distances. 
The  measured  phase  difference  between  the  surface  of  the 
metallization  and  the  surface  of  the  fused  silica  substrate  is 
8.3*,  which  corresponds  to  a  film  thickness  of  2800  A  after 
taking  into  account  the  angular  dependence  of  the  complex 
reflectance  function  and  the  effect  of  focusing. 

Providing  the  same  fractional  bandwidth  can  be  main¬ 
tained,  there  should  be  little  difficulty  in  applying  the  same 
measurement  techniques  to  acoustic  microscopes  operating 


Fig.  S.  Topography  images  of  striped  metallization  patterns  on  a  fused 
quartz  substrate.  The  metallization  is  gold  and  the  thickness  is  at>out  3000 
A.  (a)  phase  image  of  250-/zm  lines,  (b)  phase  image  of  125-/zm  lines,  (c) 
amplitude  image  of  123-zim  lines,  (d)  phase  image  of  62.3-izm  lines. 
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at  a  much  higher  frequency,  perhaps  up  to  2  GHz.  Beyond 
that  frequency,  limitations  in  the  presently  available  elec¬ 
tronic  switching  components  may  present  difficulties. 

C.  Reflectance  by  Inversion  of  complex  V(z) 

It  has  been  demonstrated  that  the  reflectance  function 
at  a  liquid/solid  interface  can  be  obtained  by  inverting  the 
corresponding  complex  V(z)  data.7,8  The  principal  motiva¬ 
tion  behind  this  endeavor  is  that  one  can  obtain  a  lot  of  useful 
material  property  information  from  the  reflectance  function. 
The  solid  curves  in  Fig.  6  show  the  magnitude  and  the  phase 
of  the  theoretical  reflectance  function  of  a  water/fused  silica 
interface.  The  reflectance  function  contains  information  on 
the  longitudinal  and  shear  critical  angles,  and  thus  the  veloc¬ 
ities  of  propagation  of  the  respective  modes  in  fused  silica. 
Also,  the  Rayleigh  critical  angle  corresponds  to  the  point  at 
which  the  phase  is  n  radians  in  the  region  where  the  phase 
curve  undergoes  a  rapid  2rr  radian  change. 

It  can  be  shown  that  the  relation  between  V(z)  and  the 
reflectance  function  R{6)  is  essentially  one  of  Fourier  trans¬ 
formation.  The  inversion  formula  is  given  by8 

R(t)  =Jr-,[K(tz)]/P2(f),  (1) 

where  u  =  z/k  and  t  =  2k z  /k.  The  parameter  k  is  the  propa¬ 
gation  constant  in  water,  and  kr  and  kz  =  ^k2  —  k  2  are  the 
corresponding  propagation  constants  in  the  r  and  z  direc¬ 
tions,  respectively.  The  parameter  /*(/)  is  the  pupil  function 
of  the  acoustic  lens;  it  describes  the  complex  amplitude  dis¬ 
tribution  of  the  focused  acoustic  beam  as  a  function  of  angle 
0  where  sin  0  —  kr/ko. 

Generally,  the  lens  pupil  function  P  is  not  known  a 
priori.  However,  it  can  be  calibrated  by  measuring  and  then 
inverting  the  complex  F(z)  for  a  material  whose  reflectance 
function  has  uniform  amplitude  and  phase  over  the  range  of 
angular  illumination  covered  by  the  acoustic  lens. 

The  V(z)  measurements  were  conducted  at  10  MHz .  A 
water/lead  interface,  whose  reflectance  function  is  uniform 
in  magnitude  and  phase  out  to  6  —  40°,  was  used  to  calibrate 
the  lens  pupil  function  P.  The  experimental  reflectance  func¬ 
tion  of  fused  silica,  which  results  from  the  inversion  proce¬ 
dure,  is  superposed  on  Fig.  6  in  dashed  lines.  The  agreement 
is  good  between  theory  and  experiment  in  both  magnitude 
and  phase.  The  measured  shear  critical  angle  at  0  =  23.5° 
and  the  Rayleigh  critical  angle  at  0  =  25.85°  compare  ex¬ 
tremely  well  with  theoretical  values. 

IV.  FOURIER  ANALYSIS  OF  THE  PHASE 
MEASUREMENT  SYSTEM  AND  NOISE 
CONSIDERATIONS 

Consider  the  signal  recovery  portion  of  the  system  in 
Fig.  3,  starting  from  where  the  acoustic  signal  is  received  to 
the  100-kHz  crystal  filter  output.  Suppose  the  acoustic  mi¬ 
croscope  is  excited  with  a  coherent  rf  pulse  train  of  the  form 

/u)  =  cos  rect(~-y~£_)  ’  (2) 

where  <y0  is  the  rf  carrier  frequency  and/,  =  \/Tp  is  the 
repetition  rate  of  the  system.  The  width  of  the  excitation 
pulse  Tw  is  made  sufficiently  short  so  that  the  acoustic  re- 
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Fig.  6.  Theoretical  and  experimental  reflec¬ 
tance  function  for  a  water/fused  silica  inter¬ 
face  obtained  by  the  inversion  of  the  corre¬ 
sponding  complex  K(z). 


turn  pulses  of  interest  are  time  separated.  Typically,  we 
would  measure  the  relative  phase  between  two  such  pulse 
trains, g,  (0  and g2  (t).  In  general,#,  (t)  andg2  (/)  are  dif¬ 
ferent  acoustic  responses  to  the  excitation  signal /( t )  and  can 
be  modeled  by 

g,(0  =/(0*M0, 

(3) 

g2(t)  —f(0*h2(t), 

where  h ,  (0  and  A  2  (0  are  the  respective  impulse  response 
functions,  and  “*”  denotes  convolution,  g,  (f)  andg2  (0  are 
separated  by  time  gating  for  processing  in  parallel  channels, 
as  shown  in  Fig.  3. 

The  frequency  spectrum  of  the  signal  in  each  processing 
channel  is  given  by 

A  sinc[  Tw  {(o  +  <u0)/2w]/f,  («y) 

G,(<u)  =  - - - 

X  ^S(a)  +  con  —  1(0, ) 

A  sinc[ Tw(co  —  co0 )/2ir]Hi(a) 


X^\5(ft>  -  <y0-  la,),  (4) 

where  H,  (a>)  is  the  Fourier  transform  of  A,  (0,  and  i  is  the 
channel  index.  The  magnitude  spectrum  | G,  (<a)  \  is  depicted 
in  Fig.  7(a). 

The  multisideband  signal  g,  (f )  is  heterodyned  down  to 
a  center  frequency  of  <y,F  =  a>0  —  coR .  The  corresponding 
spectrum  \  P,  (&>)|  is  shown  in  Fig.  7(b).  We  have  to  be  care¬ 
ful  about  aliasing  because  of  the  overlap  between  the  discrete 


spectral  islands  centered  at  ±  (oIF .  With  a  correct  choice  of 
oIF  and  the  repetition  rate  (o, ,  as  discussed  below,  the  dis¬ 
crete  frequency  components  situated  at  ±  (olF  can  be  ex¬ 
tracted  by  a  narrow  bandpass  filter  to  produce  a  cw  output 
$/(»), 

<7,(0  =  K  \H ,  {(Oq)  |cos[<y,Ff  +  ^,(<u0)],  (5) 

where  0,((oo)  =  arg  H,  (co0 ).  The  signal  q,(t)  is  in  essence  a 
low-frequency  cw  replica  of  the  rf  pulse  g,  (t)  in  the  sense 
that  it  contains  the  same  amplitude  and  phase  information  as 
the  rf  carrier  frequency  <u0  of  the  acoustic  signal.  The  signals 
q{  (f)  and  q2  (0  are  then  used  as  the  reference  and  signal 
inputs  to  a  lock-in  amplifier  where  the  relative  phase 
9 ,  (c o0 )  —  9 2  (<u0 )  *s  measured. 

It  is  evident  from  Fig.  7(b)  that  the  choice  of  the  repeti- 


(<x>o+ wp)  -  wjr  ♦  u»o  ♦  <up 

(b) 


Fig.  7.  Magnitude  spectra  (a)  of  the  gated  periodic  rf  signal,  (b)  after  mix¬ 
ing  with  the  reference  frequency. 
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tion  rate  cos  of  the  phase  measurement  system  is  crucial  to 
the  faithful  extraction  of  the  discrete  frequency  components 
at  ± &),F  in  P,(co)  .  The  following  condition  needs  to  be 
satisfied  to  avoid  aliasing: 

2<uIFt ^na)„  n  =  1,2,3 .  (6) 

Practical  narrowband  filters  have  finite  bandwidth  that 
must  be  taken  into  account  in  aliasing  considerations.  Since 
filter  response  time  is  inversely  proportional  to  the  band¬ 
width,  it  is  desirable  to  utilize  the  largest  bandwidth  possible 
in  order  to  minimize  settling  time.  Therefore,  cos  should  be 
chosen  so  that 

co,  =  2coIF/n  +  1/2  (7) 

to  give  a  maximum  bandwidth  of  co,  for  FI  and  F2  . 

Next  we  consider  the  effect  of  additive,  Gaussian  noise 
on  the  phase  resolution  of  the  measurement  system.  A  de¬ 
tailed  noise  analysis  of  the  system  is  omitted  for  brevity.  It  is 
well  known  that  the  signal-to-noise  ratio  (S/N )  is  conserved 
in  heterodyning,  and  that  removing  the  unnecessary  band¬ 
width  from  a  narrowband  signal  improves  the  S/N  . 

Let  S,/N,  be  the  signal-to-noise  ratio  of  the  acoustic 
input  signal  to  the  time  gate  SW2  in  Fig.  3.  It  can  be  shown 
that  the  S/N  at  the  output  of  the  100-kHz  crystal  filter  F2  is 
given  by  the  relation12 

A  =  _LAA,  (8) 

Nf  2  Bf  N, 

where  BF  is  the  bandwidth  of  the  filter.  Equation  ( 8 )  clearly 
shows  the  improvement  in  S/N  by  reducing  the  noise  band¬ 
width  from  B ,  to  BF,  as  c  would  expect. 

It  can  also  be  shown  that  the  corresponding  root-mean- 
square  phase  noise  at  the  output  of  the  crystal  filter  F2  is 
given  by  the  relation12 

(MW2=1  /(Sf/Nf).  (9) 

The  low-pass  filter  at  the  output  stage  of  the  lock-in 
amplifier  can  further  reduce  the  phase  noise  if  its  equivalent 
noise  bandwidth  BL  is  such  that  2 BL  <BF.  If  the  reference 
channel  input  is  noise  free,  A0  RMS  at  the  lock-in  output  is 
given  by 

\  2  mm(BF,2BL)  l  1 

(A0RMS)  =  - n - 7T7T7-  =  ~F~7TT  ’ 

Bf  Sf/Nf  Sq/ n o 

(10) 

where  S0/N0,  defined  as  the  equivalent  output  signal-to- 
noise  ratio,  is 

A.  =  J - El - II,  (id 

N0  2  mm (Bf,2Bl)  N, 

When  the  reference  channel  input  is  also  noisy,  the 
phase  noise  power  spectrum  is  simply  the  sum  of  the  power 
spectra  of  the  signal  and  reference  channel  inputs 

( A^rms  ) TOTAL  ~  (  A^RMS  )  SIG  +  (i^RMS  1  REF 

=  _ I _  +  _ I _  =  _ I _ 

VA,„  -V(v<w  V^OIA1  ’ 


where  {Sn/N0  )REF  is  the  equivalent  output  S/N  of  the  ref¬ 


Fig.  8.  Experimentally  measured  phase  noise  as  a  function  of  the  equivalent 
input  noise. 


erence  channel  input,  and  ( S0/N0  )TOtal  is  the  total  equiva¬ 
lent  output  S/N. 

Figure  8  shows  the  experimentally  measured  phase 
noise  as  a  function  of  S,  /Nf .  The  measurement  system  is  set 
up  for  B i  =  300  MHz,  BF  =  100  Hz,  BL  —  25  Hz;  only  the 
signal  channel  input  is  noisy.  The  input-to-output  S/N  im¬ 
provement,  as  given  by  Eq.  ( 1 1 ),  is  65  dB,  which  is  within  2 
dB  of  our  measurement.  The  experimental  data  exhibit 
clearly  the  inverse  relationship  between  A0RMS  and  S,  /Ni  as 
predicted  by  theory. 

The  typical  operating  parameters  for  the  50-MHz 
acoustic  microscope  system  are  BF  —  100  Hz,  pulse  width 
Tw  =  100  ns,  and  repetition  period  Tp  =  64  /zs.  The  domi¬ 
nant  source  of  noise  is  the  thermal  noise  in  the  preamplifier 
following  SW3  in  Fig.  3.  Suppose  the  preamplifier  has  a 
noise  bandwidth  of  5,  =  100  MHz  and  a  noise  figure  (NF) 
of  10  dB.  Let  us  assume  that  there  is  a  further  10-dB  degrada¬ 
tion  in  S/N  ratio  due  to  noise  contribution  along  the  signal 
path  from  the  inputs  of  SW1  and  SW2  to  the  inputs  of  the 
100-kHz  crystal  filters  F3  and  F4.  The  effective  NF  referred 
to  the  input  of  the  preamplifier  is  therefore  20  dB.  At  the 
preamplifier  input,  typically  the  signal  level  A  is  about  5  mV. 
The  effective  input  S/N  ratio  is  given  by 


1l  =  A  2/2  III 

N,  Fpx4kTB,R  Tp  ’ 


(13) 


where  Fp  is  the  total  noise  figure  (power)  and 
k=  1.38x10' 23  J/K  is  Boltzmann’s  constant.  At  room 
temperature  ( T  =  300  K),  and  for  a  50-fl  system,  S,/N,  is 
about  +  4  dB.  The  corresponding  S/N  ratio  at  the  crystal 
filter  output  is  +  62.7  dB  from  Eq.  ( 8 ) ,  which  translates  to  a 
root-mean-square  phase  noise  of  A#RMS  =  0.04°  according 
to  Eq.  ( 10).  Suppose  that  we  have  the  same  amount  of  phase 
noise  at  the  outputs  of  F3  and  F4.  The  theoretical  phase 
resolution  limit  as  given  by  Eq.  (12)  is  0.06°.  Thus,  under 
typical  measurement  conditions,  we  are  operating  fairly 
close  to  the  phase  resolution  limit  of  the  lock-in  amplifier, 
which  is  0.1°.  In  certain  situations  where  the  input  S/N  ratio 
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can  be  improved  by  increasing  the  signal  duty  cycle  Tw/Tp, 
as  in  V(z)  measurements  where  signal  pulses  are  several  mi¬ 
croseconds  long,  the  lock-in  amplifier  becomes  the  limiting 
factor  in  the  phase  resolution  of  the  measurement  system. 
The  lock-in  amplifier  has  to  be  replaced  by  more  precise 
phase  measurement  instrumentations  in  such  cases  in  order 
to  take  full  advantage  of  the  corresponding  reduction  in 
phase  noise. 


V.  DISCUSSION 

We  have  described  the  basic  principle  and  demonstrated 
the  operation  of  a  precision  phase  measurement  system.  Al¬ 
though  the  system  is  specifically  designed  for  acoustic  mi¬ 
croscopy,  its  application  can  readily  be  extended  to  situa¬ 
tions  where  precision  phase  measurement  is  required  for 
synchronous  periodic  rf  pulses  as  short  as  two  or  three  cy¬ 
cles. 
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An  electronically  scanned  optical  microscope  which  quantitatively  measures  amplitude  and 
phase  is  described.  The  system  is  insensitive  to  mechanical  vibrations.  The  phase  information 
makes  it  possible  to  measure  surface  height  variations  with  an  accuracy  of  better  than  100  A  and 
can  also  be  used  to  improve  the  lateral  resolution. 


Optical  microscopy  is  widely  used  as  a  method  of  evalu¬ 
ating  small  scale  surface  features.  In  examining  these  fea¬ 
tures,  both  the  lateral  dimensions  as  well  as  the  height  are  of 
interest.  Conventional  optical  microscopy  gives  little  quanti¬ 
tative  information  on  surface  height.  Nomarski  differential 
phase  contrast,1  as  well  as  conventional  phase  contrast  mi¬ 
croscopy,2  provides  qualitative  contrast  for  objects  differing 
from  their  surroundings  in  optical  phase  rather  than  ampli¬ 
tude.  However,  it  is  very  difficult  to  extract  precise  height 
information  from  these  images.3  Interference  microscopy4 
involves  the  counting  of  fringes  to  evaluate  surface  height, 
but  again  the  fringe  maps  are  difficult  to  interpret. 

We  present  here  an  electronically  scanned  optical  mi¬ 
croscope  which  directly  measures  the  optical  amplitude  and 
phase  of  light  reflected  from  a  sample.  The  system  incorpo¬ 
rates  internal  references  which  make  it  insensitive  to  envir¬ 
onmental  vibrations  which  are  a  common  problem  in  inter¬ 
ferometers.  Information  from  the  linear  scan  is  obtained  in  a 
form  which  is  directly  compatible  with  digital  processing. 
Since  both  amplitude  and  phase  are  measured,  it  is  possible 
to  take  a  Fourier  transform  of  the  complex  spatial  variation. 
This  can  be  used  to  remove  aberrating  effect  in  the  optical 
system  by  Weiner  filtering  with  the  measured  point  spread 
function  of  the  system.  This  should  make  it  possible  to  more 
accurately  determine  the  lateral  dimensions  of  surface  fea¬ 
tures,5  such  as  linewidths  on  intefrated  circuits. 

On  a  sample  of  uniform  composition,  the  optical  phase 
can  be  simply  calibrated  to  yield  the  height  of  surface  fea¬ 
tures.  Thus,  the  scanning  microscope  can  be  used  as  a  non¬ 
contacting  optical  profilometer  which  provides  height  infor¬ 
mation  similar  to  that  obtained  using  a  mechanical  stylus 
profilometer.  On  inhomogeneous  samples,  compensation  is 
required  for  phase  and  polarization  changes  on  reflection. 

Earlier  work  on  optical  profilometry  using  a  differential 
technique,6  or  by  heterodyning,  using  a  Zeeman  laser7  does 
not  make  use  of  electronic  scanning.  In  our  system,  a  Bragg 
cell  is  used  to  frequency  shift  as  well  as  to  scan  the  optical 
beam.  Line  scan  rates  of  over  10  kHz  are  possible  with  the 
present  Bragg  cell,  provided  the  detection  electronics  can  be 
designed  to  be  sufficiently  fast.  This  makes  rapid  image  pro¬ 
cessing  a  possibility. 

Two  schematics  of  the  optical  system  are  shown  in  Fig. 
1.  A  simplified  version  of  the  basic  system,  shown  in  Fig. 
1(a),  illustrates  how  a  Bragg  cell  can  be  used  to  split  and 
frequency  shift  the  optical  beam.  Similar  optical  heterodyn¬ 
ing  with  a  Bragg  cell  has  been  used  to  detect  acoustic  surface 
vibrations.*  The  acoustic  intensity  is  adjusted  so  that  the 
input  beam  is  split  equally  between  an  undiffracted  optical 


beam  of  frequency and  a  Bragg  diffracted  beam  upshift ed 
to  frequency  f0  + fB,  where  f B  is  the  frequency  of  the  input 
signal  to  the  Bragg  cell.  A  combination  of  a  field  lens  to 
image  the  center  place  of  the  Bragg  cell  to  the  back  focal 
plane  of  a  microscope  objective,  and  the  objective  itself,  con¬ 
verts  the  angular  deviation  produced  by  the  Bragg  cell  into  a 
displacement  of  one  of  the  optical  spots  on  the  sample.  After 
reflection  from  the  sample,  the  beams  retrace  their  paths  and 
reenter  the  Bragg  cell.  Half  of  the  previously  undiffracted 
beam  is  now  down  shifted  by  the  Bragg  cell  to  f0  —  fB  and 
falls  on  a  detector.  In  addition,  half  of  the  previously  de¬ 
flected  beam  passes  through  the  Bragg  cell  with  no  further 
frequency  shift  and  follows  the  same  path  to  the  detector. 
Thus,  two  signals  at  frequencies /„  +  fB  and  f0  —  fB  interfere 


FIG.  1,  Scanning  optical  microscope.  The  Bragg  cell  splits  the  input  beam 
and  scans  the  deflected  and  frequency  shifted  beam  across  the  sample.  The 
undiffracted  beam  is  used  as  a  reference,  {al  Simplified  schematic,  (b)  Secon¬ 
dary  reference  on  sample  to  compensate  for  surface  tilt. 
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at  the  detector  to  produce  a  difference  frequency  of  2/„ .  The 
phase  of  the  signal  at  2  f„  depends  on  the  optical  path  length 
difference  between  the  two  beams  impinging  on  the  sample, 
since  the  optical  signals  at/,  4-  /„  and/,  —  fH  originate  from 
the  scanned  and  stationary  spots  on  the  sample,  respectively. 

For  plane  wave  illumination,  a  surface  height  variation 
of  h  will  introduce  an  optical  phase  change  of  2kh  as  the 
beam  is  scanned  over  a  feature  of  height  h.  For  a  tightly 
focused  beam,  the  observed  phase  change  is  less  since  much 
of  the  incident  light  arrives  at  the  sample  at  off-normal  inci¬ 
dence.  By  integrating  over  a  uniformly  illuminated  spherical 
lens  aperture,  it  can  be  shown  that  the  feature  introduces  a 
phase  change  of 

Atf>  =  kh(\  +  cos  0„),  (1) 

where  0„  is  the  half-angle  of  the  objective  aperture,  and 
k  =  2ir/A.  As  0O  approaches  zero,  this  formula  reduces  to 
the  phase  change  for  a  plane  wave.  Changes  in  the  distance 
from  the  lens  to  the  sample  do  not  affect  the  measurement 
since  only  path  length  differences  between  the  two  beams  are 
important.  Thus,  mechanical  vibrations  in  the  lens  spacing 
do  not  disturb  the  measurements  so  long  as  they  are  less  than 
the  depth  of  focus  of  the  objective.  In  the  simplest  configura¬ 
tion  [Fig.  1(a)],  the  measurement  is  sensitive  to  surface  tilt. 

To  measure  the  phase  shift  in  the  interference  signat  at 
2/b  as  the  Bragg  cell  frequency  fB  is  varied  over  the  scan 
range,  an  electronic  reference  signal  is  required.  Simply  fre¬ 
quency  doubling  the  Bragg  cell  drive  signal  is  not  sufficient 
since  the  finite  acoustic  delay  in  the  Bragg  cell  will  produce  a 
large  linear  phase  shift  as  the  Bragg  frequency  is  changed. 
Instead,  a  second  optical  signal  is  generated  [Fig.  1(b)]  by 
splitting  the  beam  incident  on  the  Bragg  cell  in  a  direction 
perpendicular  to  the  direction  of  acoustic  wave  propagation. 
Thus,  four  spots  are  focused  on  the  sample.  If  the  sample 
geometry  is  chosen  so  that  the  second  set  of  spots  is  located 
on  a  uniform  flat  section  of  the  sample,  then  they  can  be  used 
for  a  reference.  Alternatively,  the  second  set  of  spots  can  be 
magnified  to  be  much  larger  than  the  size  of  surface  features 
of  interest  and  placed  anywhere  on  the  sample. 

The  Bragg  cell  is  imaged  on  a  photodiode  detector  to 
compensate  for  the  variation  of  the  Bragg  angle  with  fre¬ 
quency.  Due  to  the  use  of  a  coherent  detection  scheme,  the 
system  is  a  type  II  confocal  microscope.'*  The  Rayleigh  reso¬ 
lution  of  0.57  A  /NA  is  slightly  better  than  the  resolution  of 
0.61  A  /NA  of  the  incoherent  microscope. 10  The  ripples  ob¬ 
served  while  scanning  across  a  discontinuity  are  reduced 
over  that  of  an  incoherent  type  I  system. 10 

Figure  2  shows  the  signal  processing  used  to  detect  the 
amplitude  and  phase  of  the  optically  derived  signals  at  2 fB . 
Since  it  is  advantageous  to  detect  phase  at  a  fixed  frequency, 
the  reference  signal  is  first  limited  and  then  single  sideband 
modulated  at  an  intermediate  frequency  of  60  MHz.  Mixing 
this  reference  with  the  signal  photodetector  output  gives  a 
fixed  frequency  signal  at  60  MHz  which  contains  all  the  am¬ 
plitude  and  phase  information.  In  our  initial  experiments,  a 
vector  voltmeter  is  used  to  compare  this  signal  to  an  elec¬ 
tronic  reference. 

Calibration  scans  are  performed  on  an  optical  flat  and 
amplitude  and  phase  variations  due  to  the  electronics  as  well 
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FIG.  2.  Phase  measurement  system.  Phase  detection  is  performed  at  a  fixed 
frequency  of  60  MHz. 

as  optical  aberrations  are  divided  out  and  subtracted  out, 
respectively.  After  this  calibration,  scans  across  a  flat  uni¬ 
form  sample  show  amplitude  variations  less  than  3%  and 
phase  variations  of  less  than  5°,  suggesting  a  minimum  sur¬ 
face  height  sensitivity  of  better  than  100  A. 

In  our  experiments,  5  mW  from  an  etalon  controlled 
single-frequency  argon-ion  laser  at  510  nm  was  used.  The 
Te02  Bragg  cell10  has  a  bandwidth  of  60-1 10  MHz  and  de¬ 
flection  angles  of  3°-6l>.  Microscope  objectives  with  numeri¬ 
cal  apertures  in  the  range  of  0. 1-0.86  were  used.  The  length 
of  the  scan  on  the  sample  is  proportional  to  the  frequency 
change  times  the  objective  focal  length.  With  the  higher  nu¬ 
merical  aperture,  a  1/e2  spot  radius  of  1  /jm  was  measured  by 
examining  the  reflection  from  a  metallized  edge. 

Figure  3  demonstrates  an  initial  experiment  on  the  ap¬ 
plication  of  the  microscope  to  optical  profilometry.  In  order 
to  check  the  system  performance,  our  aim  was  to  eliminate 
phase  changes  due  to  the  use  of  a  wide  aperture  lens  caused 
by  polarization  effects  and  changes  in  phase  of  the  reflection 
coefficient.  A  step  of  aluminum  900  A  thick  (as  measured 
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FIG.  3.  Amplitude  and  phase  variation  across  an  aluminum  step  on  a  metal¬ 
lized  substrate.  The  phase  step  indicates  a  height  of  820  A. 
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with  a  mechanical  stylus)  was  deposited  on  a  glass  substrate 
metallized  with  500  A  of  aluminum.  To  simulate  plane  wave 
illumination,  a  long  focal  length  microscope  objective  (16- 
mm  focal  length)  was  illuminated  with  a  narrow  beam  less 
than  1  mm  in  diameter.  This  yields  an  effective  numerical 
aperture  of  less  than  0.03  and  ^„=;0  in  Eq.  ( 1 ).  The  Bragg  cell 
frequency  was  scanned  over  the  range  of  82.0-87.8  MHz, 
causing  the  focused  beam  to  scan  across  90fim  on  the  sam¬ 
ple.  In  Fig.  3,  the  amplitude  and  phase  of  the  scanned  spot 
are  shown.  After  calibration,  5’  of  phase  offset  has  been  sub¬ 
tracted  from  the  experimental  phase  curve.  The  theoretical 
curves  are  evaluated  numerically  for  a  Gaussian  beam  with  a 
1/e2  radius  of  9  fj. m,  which  is  consistent  with  the  small  nu¬ 
merical  aperture  of  the  lens  used.  The  beam  is  assumed  to 
have  constant  optical  phase  and  is  convolved  across  the 
phase  step.  Due  to  aberrations  and  misfocusing,  there  may 
be  phase  variations  across  the  beam  which  could  account  for 
the  discrepancy  between  the  experimental  and  theoretical 
amplitude  curves.12 

The  phase  curve  indicates  a  step  height  of  820  A,  which 
is  in  fair  agreement  with  the  mechanical  stylus  measurement 
of 900  A,  and  demonstrates  the  ability  to  obtain  good  quanti¬ 
tative  contrast  on  samples  with  features  differing  only  in  op¬ 
tical  phase. 

In  conclusion,  an  electronically  scanned  optical  micro¬ 
scope  which  can  quantitatively  measure  amplitude  and 
phase  is  demonstrated.  Surface  height  variations  can  be  ac¬ 
curately  measured  with  the  system.  In  addition,  the  ability  to 


measure  optical  phase  promises  to  make  possible  digital  fil¬ 
tering  to  remove  the  effect  of  aberrations  and  to  accurately 
determine  lateral  dimensions  on  the  sample. 

The  authors  wish  to  thank  B.T.  Khuri-Yakub,  S.D. 
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their  informative  discussions.  E.  Roos  provided  the  Bragg 
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Air  Force  Office  of  Scientific  Research  under  contract  No. 
F49620-79-C-02 1 7. 
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The  on-axis  intensity  response  of  a  confocal  scanning  optical  microscope  was  measured  for  an  objective  of  numerical 
aperture  0.9.  The  data  compare  favorably  with  theoretical  calculations  obtained  by  numerical  integration  of  the 
standard  theory,  provided  that  lens  aberrations  are  taken  into  account.  The  invariance  of  the  shape  of  the  central 
lobe  to  surface  roughness  and  tilt  is  also  demonstrated. 


The  confocal  scanning  optical  microscope  (CSOM)1 
has  been  the  subject  of  much  recent  study.  One  of  the 
potentially  most  useful  properties  of  this  type  of  mi¬ 
croscope  for  making  film-thickness  and-profile  mea¬ 
surements  is  its  shallow  depth  of  Held.  We  have  mea¬ 
sured  the  intensity  response,  |V(z)l2,  of  a  CSOM  as  a 
function  of  the  amount  of  defocus  of  a  plane  reflector 
from  the  focal  plane.  We  have  found  good  agreement 
between  the  experimental  results  and  the  theory.  We 
have  also  found  that  the  shape  of  the  central  lobe  is 
invariant  to  surface  roughness  and  tilt  in  the  sample. 
This  property  is  potentially  useful  for  remote  distance 
sensing  and  other  applications. 

The  amplitude  V(z)  of  the  wave  reaching  a  coherent 
detector  in  a  confocal  system  has  been  determined  by 
Liang  et  al 2  and  Sheppard  et  al.3  for  the  acoustic 
microscope  and  in  scalar  form  by  Wilson  and  Shep¬ 
pard4  for  the  scanning  optical  microscope.  We  have 
derived  a  nonparaxial  form  of  the  vector  field  theory. 
This  can  be  summarized  as  follows. 

We  assume  a  lens  with  a  pupil  function  P(6),  where  0 
is  the  angle  between  the  ray  from  the  pupil  plane  to  the 
focal  point  and  the  axis.  By  following  the  derivation 
of  Richards  and  Wolf,5  with  some  minor  changes,  it 
can  be  shown  that  the  transverse  electric  field  associ¬ 
ated  with  a  plane  wave  focused  by  the  lens  to  a  point 
on  axis  a  distance  z  beyond  the  focal  plane  is  of  the 
form 

E(8,  z)  =  [  —  +  C°8  g  exp(ifez  cos  0)P(0)d0. 

Jo  (cos  8yn 

(1) 

It  will  be  noted  that  there  is  a  (cos  0)1/2  term  in  the 
denominator  rather  than  in  the  numerator,  as  in  the 
theory  of  Richard  and  Wolf,  in  order  to  conserve  power 
at  a  flat  exit  plane  from  the  lens. 

The  integrand  of  Eq.  (1)  expresses  the  amplitude 
and  phase  of  the  plane  wave  components  of  the  E  field 
at  an  angle  8  to  the  axis.  When  a  focused  beam  is 
reflected  from  a  plane  mirror,  its  image  will  be  a  dis¬ 
tance  2 z  away  from  the  focus.  In  our  apparatus,  the 
reflected  image  is  refocused  onto  the  pinhole  in  front 
of  the  detector.  Thus  it  is  the  field  on  axis  at  the  focal 
plane  of  the  objective  lens  that  is  imaged  at  the  pin¬ 
hole.  This  field  is  of  the  form 


'9°  (1  +  cos0)sin  8 
.o  (cos  6) 1/2 


X  exp(— 2ikz  cos  0)P(0)R(8)A8.  (2) 

In  these  expressions,  R(8)  is  the  reflection  coeffi¬ 
cient  of  the  plane  reflector,  k  is  the  wave  number  (2jt/ 
X),  f  is  the  focal  length  of  the  lens,  and  sin  do  is  the 
numerical  aperture  of  the  lens.  Given  a  lens  of  radius 
a  with  uniform  excitation,  [P(0)  *  1  for  8  <  0O]  for  large 
ka  and  a  plane  reflector  [i?(0)  =  1].  An  approximate 
expression  for  |  V(z)\  has  been  derived  by  Liang  et  al.2: 


IV(z)|  = 


sin  kz(  1  -  cos0o) 
fez(l  -  COS0O) 


This  expression  accurately  predicts  the  shape  of  the 
central  lobe;  however,  it  fails  to  account  for  asymme¬ 
tries  in  the  sidelobes.  The  depth  of  focus  is  given  by 
the  3-dB  points  of  the  central  lobe  in  Eq.  (3)  (Ref.  2): 


(Az)3dB  = 


0.45X 
1  -  cos  0O 


(4) 


The  basic  technique  for  measuring  optical  IV(z)l2 
was  outlined  by  Cox  et  al.6  In  Fig.  1,  the  laser  (X  =  633 
nm),  the  beam  expander,  and  the  iris  diaphragm  pro- 


Fig.  1.  Schematic  diagram  of  apparatus.  H.V.,  high-volt¬ 
age  power  supply  to  drive  scanning  stage;  D/A,  digital-to- 
analog  converter;  A/D,  analog-to-digital  converter. 
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Fig.  2.  I  V(z)\2  curve  for  objective  of  N.A.  0.9.  Experimental  data,  dotted  curve.  Theoretical  calculations,  solid  and  dashed 
curves. 


duce  plane  waves.  The  Bragg  cell  modulates  the  sig¬ 
nal  so  that  ac  detection  techniques  may  be  employed. 
The  plane  waves  are  focused  by  a  microscope  objective 
onto  an  aluminum  front-surface  mirror.  The  light 
reflected  from  the  mirror  is  then  focused  by  another 
lens  onto  a  pinhole  and  eventually  onto  a  photodetec¬ 
tor.  The  output  of  the  photodetector  is  measured  by  a 
lock-in  amplifier.  The  mirror  is  mounted  on  a  piezo¬ 
electric  transducer  so  that  it  can  be  scanned  in  the 
axial  (z)  direction.  Its  relative  position  is  measured 
by  using  a  linear  variable  differential  transformer 
(LVDT). 

The  apparatus  is  computer  controlled.  The  com¬ 
puter  reads  the  mirror  position  and  signal  intensity. 
It  then  increments  the  voltage  to  the  piezoelectric 
transducer,  thereby  moving  the  mirror  and  once  again 
taking  a  data  set.  The  computer  repeats  this  opera¬ 
tion  for  the  full  range  of  the  piezoelectric  transducer 
(about  12  m),  taking  data  points  approximately  every 
0.02  nm. 

We  present  data  for  a  microscope  objective  with  a 
numerical  aperture  (N.A.)  of  0.9  in  air  (dotted  curve  in 
Fig.  2).  The  z  axis  (position)  is  taken  to  be  zero  at  the 
position  of  best  focus,  with  positive  values  of  z  repre¬ 
senting  increased  objective-sample  separation.  The 
0.9-aperture  data  are  similar  to  data  published  by  Cox 
et  al.6 

The  theoretical  data  were  calculated  by  numerical 
integration  of  Eq.  (2).  In  the  integration,  a  uniform 
R(9)  was  assumed.  If  the  nonuniform  R(6)  for  alumi¬ 
num  is  used,  a  more  complex  and  time-consuming 
integral  is  generated.  When  evaluated,  this  integral 
does  not  noticeably  change  the  results  of  Fig.  2.  If  a 
uniform  pupil  function  P(9)  is  assumed,  the  calculated 
I  V(z)|2  curves  are  nearly  symmetrical  about  the  posi¬ 
tion  of  best  focus,  as  shown  by  the  dashed  lines  in  Fig. 


2.6  This  assumption,  however,  is  not  valid  in  our  ap¬ 
paratus.  Optical  reflections  due  to  antireflection¬ 
coating  mismatches  at  633  nm  within  a  thick  lens 
should  generate  some  amplitude  apodization.  In  ad¬ 
dition,  we  expect  that  with  a  system  focused  on  axis, 
there  will  be  some  spherical  aberration.  We  therefor e 
used  a  trial  apodization  function  of  the  form 

P(9)  =  exp(A  sin2  0)exp(—jB  sin4  6)  ( 6  <  0O), 

P(0)  =  0  (0  >  90),  (5) 

where  A  is  the  amplitude  coefficient  and  B  is  related  to 
the  Seidel  coefficient  for  spherical  aberration.  We 
used  sin  9  rather  than  9  in  these  approximations  since 
the  aperture  is  so  large.  We  tried  different  values  of  A 
and  B  to  obtain  the  best  fit  with  the  experimental 
results.  The  calculated  result  for  A  =  0.2  and  B  =  6  is 
shown  by  the  solid  curve  in  Fig.  2. 

The  phase  error  at  the  edge  of  the  lens  corresponds 
to  1.25ir,  well  within  the  expected  values  for  primary 
spherical  aberration.7  The  value  of  A  indicates  that 
the  amplitude  is  increased  by  17%  at  the  edge  of  the 
lens.  These  curves  illustrate  that  a  relatively  small 
amount  of  aberration  can  give  rise  to  asymmetry  in  the 
I  V(z)|2  curves.  With  no  aberrations,  the  W(z)\2  curves 
are  symmetric.  However,  with  aberrations,  the  I  V(z)|2 
curves  become  asymmetric  around  z  =  0,  and  the 
asymmetry  reverses  its  shape  with  the  sign  of  the  aber¬ 
ration. 

An  interesting  and  important  feature  of  the  I  V(z)|2 
curves,  for  sensing  applications  with  low-numerical- 
aperture  long-focal-length  lenses,  is  their  insensitivity 
to  surface  roughness  and  tilt.  This  insensitivity  arises 
because  the  point  illuminated  on  the  reflector  is  still  a 
point,  even  if  the  reflector  is  a  rough  surface.  In  this 
case,  the  reflected  light  is  scattered  over  a  wide  angu- 
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Fig.  3.  V(z)  curve  for  a  rough  substrate  tilted  at  an  angle  of  5°  (solid  curve).  Reflections  from  a  mirror,  dashed  curve. 


lar  range,  but  an  image  of  the  focal  point,  albeit  re¬ 
duced  in  intensity,  is  obtained  at  the  pinhole.  In  our 
apparatus,  the  mirror  in  Fig.  1  was  replaced  with  a 
quartz  flat,  which  had  been  roughened  using  5-pm  grit 
and  coated  with  150  nm  of  gold.  This  sample  is  typi¬ 
cal  of  a  machined  surface.  The  sample  was  then  tilted 
at  an  angle  of  5°  from  a  plane  normal  to  the  z  axis  of 
our  system. 

A  |V(z) I2  trace  was  made  using  a  N.A.  =  0.07  lens 
(solid  curve  in  Fig.  3).  A  0.07  lens  was  used  because  it 
had  a  spot  size  considerably  larger  than  the  average 
surface  roughness  and  hence  would  generate  a  diffuse 
reflection.  The  central  lobe  of  this  curve  is  unchanged 
from  the  central  lobe  obtained  using  an  aluminum 
mirror,  as  shown  by  the  dashed  curve  in  Fig.  3.  The 
peak  value  is,  however,  66  dB  below  the  value  obtained 
with  a  mirror,  indicating  that  this  signal  is  generated 
primarily  by  diffuse  reflection.  The  result  in  Fig.  3 
also  gives  a  signal-to-noise  ratio  of  approximately  90 
dB  for  this  system.  The  stability  and  sharpness  of  the 
|V(z)|2  curves  make  these  curves  useful  in  remote  dis¬ 
tance  sensing  applications  on  a  variety  of  surfaces. 

A  depth-of-focus  curve,  I  V(z)|2,  for  a  CSOM  has 
been  presented  for  a  numerical  aperture  of  0.9.  Theo¬ 
retical  calculations  based  on  reasonable  assumptions 
are  presented  as  well  and  are  shown  to  be  in  good 
agreement  with  the  experimental  data.  The  |V(z)|2 


curves  are  invariant  to  surface  roughness  and  tilt. 
This  feature  has  potentially  far-reaching  applications 
in  the  areas  of  optical  range  finding  for  distances  from 
micrometers  to  meters. 

This  research  was  supported  by  the  U.S.  Air  Force 
Office  of  Scientific  Research  under  contract  no. 
AFOSR-84-0063C.  The  research  of  T.  R.  Corle  was 
supported  by  an  IBM  Graduate  Predoctoral  Fellow¬ 
ship  for  the  duration  of  the  research. 
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INTRODUCTION 

Recently,  a  great  deal  of  interest  has  been  shown  in  making  accur¬ 
ate  range  measurements  with  good  transverse  definition.  This  capability 
makes  it  possible,  in  machine  vision  systems,  to  extract  geometrical 
shape  information  from  the  Images.  In  robot  position  sensing,  It  is  Im¬ 
portant  to  determine  the  absolute  distance  Instead  of  distance  change  so 
that  noncontlnuous  measurements  can  be  made  without  the  need  for  cali¬ 
bration  at  start-up.  A  third  application  of  great  importance  Is  to 
measure  the  shape  and  size  of  machined  parts  with  a  noncontacting 
sensor. 

Research  in  the  area  of  distance  measurements  has  been  carried  out 
for  many  years  using  several  different  techniques.1  One  technique  is  to 
use  interferometry,  which  requires  a  high-quality  laser.  The  method 
yields  excellent  range  accuracy,  but  suffers  from  phase  wraparound; 
thus,  without  additional  complexity,  such  as  the  use  of  two  laser  fre¬ 
quencies  or  the  employment  of  extremely  high-frequency  modulation,  it  is 
not  possible  to  make  noncontlnuous  distance  measurements.  An  alterna¬ 
tive  technique  is  to  use  a  triangulation;  this  method  does  not  usually 
provide  very  good  transverse  resolution  and  is  attractive  only  when  the 
accuracy  requirement  is  on  the  order  of  100  ym  or  greater. 

In  this  paper,  we  describe  a  new  optical  technique  which  is  based 
upon  the  type  II  microscope.  With  this  system,  we  have  made  absolute 
distance  measurements  with  a  range  accuracy  of  2  ym  and  a  transverse 
definition  of  10  ym  at  a  working  distance  of  15  cm  .  The  system  is 
stable,  easy  to  align,  and  largely  insensitive  to  the  tilt  and  roughness 
of  the  object. 


THEORY  AND  EXPERIMENTAL  SET-UP 

A  schematic  diagram  of  the  optical  range  finder  is  shown  in  Fig.  1. 
The  objective  lens  focuses  the  laser  beam  to  a  spot  of  the  order  of 
5  ym  diameter,  so  that  most  of  the  beam  passes  through  a  10  ym 
diameter  pinhole.  A  camera  lens  placed  beyond  the  pinhole  focuses  the 
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Fig.  1.  Experimental  arrangement  of  optical  range  finder. 


light  onto  a  point  on  the  object.  An  image  of  this  point  is,  in  turn, 
focused  on  the  pinhole.  The  reflected  light  passes  back  through  the 
pinhole  and  is  detected  by  the  photodetector.  When  the  object  is  lo¬ 
cated  exactly  at  the  focus  of  the  transmitted  light,  the  maximum  amount 
of  light  passes  through  the  pinhole.  Otherwise,  when  the  system  is  par¬ 
tially  defocused,  the  amount  of  light  passing  through  the  pinhole  is 
reduced.  We  therefore  expect  the  signal  output  from  the  detector  to 
have  a  strong  dependence  on  the  position  of  the  object. 

Because  a  focused  beam  is  employed  and  the  phase  changes  between 
the  different  rays  comprising  the  beam  are  very  small,  the  phase  fluc¬ 
tuations  of  the  laser  are  unimportant1,  therefore,  a  semiconductor  laser, 
or  even  an  incoherent  source,  can  be  used.  Furthermore,  the  employment 
of  a  focused  beam  implies  that  the  transverse  definition  of  the  system 
is  excellent. 

It  will  be  noted  that  the  system  shown  is  deceptively  simple.  In 
confocal  scanning  microscopy,  it  is  common  to  employ  two  pinholes  so  as 
to  eliminate  the  reflected  light  from  the  transmitting  pinhole,  or  to 
use  a  collimated  input  beam  to  the  objective  lens  instead  of  a  pin¬ 
hole. In  both  cases,  the  light  passes  through  the  pinhole  only  once. 
The  optical  alignment  procedure  required  for  an  optical  range  finder, 
which  has  a  tmich  longer  working  distance  than  the  microscope,  is  ex¬ 
tremely  difficult  with  such  a  two-pinhole  system.  Thus,  our  system  uses 
the  single  pinhole  technique  in  which  the  light  beam  passes  through  the 
pinhole  twice.  It  is  easy  to  align  and  very  stable.  In  this  case,  the 
light  reflected  back  from  the  10  pm  pinhole  is  eliminated  by  focusing 
the  beam  to  a  spot  size  of  5-6  um  . 


We  assume  a  lens  with  a  pupil  function  P(9)  where  9  is  the 


angle  between  the  ray  from  the  pupil  plane  to  the  focal  point  and  the 
axis.  By  following  the  derivation  of  Richards  and  Wolf"’  for  vector 
fields,  with  some  minor  changes,  it  can  be  shown  that  the  transverse 


electric  field  associated  with  a  plane  wave  focused  by  the  lens  to  a 
point  on  axis  a  distance  z  beyond  the  focal  plane  is  of  the  form 


■/  — 


cos  9)  sin  9  „  . 

_ ' _  jkz  cos  9 


cos  9) 


P(9)  d9 


1/2 

It  will  be  noted  that  there  is  a  (cos  9)  '  terra  in  the  denominator 
rather  than  in  the  numerator,  as  it  is  in  Wolf's  theory,  in  order  to 
conserve  power  at  a  flat  exit  plane  from  the  lens. 


? 


The  Integrand  of  Eq.  (1)  expresses  the  amplitude  and  phase  of  the 
plane  wave  components  of  the  E  field  at  an  angle  9  to  the  axis.  When 
a  focused  beam  is  reflected  from  a  plane  mirror,  its  image  will  be  a 
distance  2z  away  from  the  focus.  In  our  apparatus,  the  reflected  im¬ 
age  is  refocused  onto  the  pinhole  in  front  of  the  detector.  Thus,  it  is 
the  field  on  axis  at  the  focal  plane  of  the  objective  lens  which  is  im¬ 
aged  at  tbe  pinhole.  This  field  is  of  the  form 


V(z) 


f  (1  +  cos  9)  sin  9 
J  ( cos  9) 


ITT 


:2jkz  cos  p(e)  R(0) 


d9 


(2) 


where  the  output  signal  from  the  detector  Is  proportional  to 


I(z)  -  |v (z)|  2 


(3) 


In  these  expressions,  R(0)  is  the  reflection  coefficient  of  the 
plane  reflector,  k  is  the  wave  number  (2u/A)  ,  f  is  the  focal  length  of 
the  lens,  and  sin  0Q  is  its  numerical  aperture,  where  sin  9q  »  a/f  , 
and  the  radius  of  the  lens  is  a  .  With  uniform  excitation,  [p(9)  “  1 
for  9  <  0q ]  ,  ka  »  100  ,  and  a  plane  reflector  [R(9)  -  1]  •  Ag 
approximate  expression  for  |V(z)|  has  been  deriv-  4  by  Liang  et  al 


|v(z)| 


sin  kz(l  -  cos  9q) 
kz(l  -  cos  9q) 


(4) 


This  expression  accurately  predicts  the  shape  of  the  central  lobe; 
however,  it  falls  to  account  for  asymmetries  in  the  sidelobes.  The 
depth  of  focus  is  given  by  the  3  dB  points  of  the  central  lobe  in 
Eq.  (4). 


(Az) 


3  dB 


0.443  A 


1  -  cos  9q 


(5) 


For  small  ,  the  numerical  aperture  (N.A.)  of  the  lens  is2given  by 

the  relation  [N.A.  -  sin  (0Q)  -  0Q]  .  With  cos  eQ  -  1  -  9Q/2  ,  we  can 
write  Eq.  (5)  in  the  form 


(iz)3  dB 


0.886  A 

(N.A.)2 


(6) 


We  used  a  4  mW  He-Ne  laser  with  a  wavelength  of  6328  A  in  these 
experiments.  A  semiconductor  laser  could  have  been  employed  equally 
well;  we  chose  to  work  with  the  gas  laser  Initially  only  because  of  its 
availability  and  the  greater  ease  of  working  with  visible  light  in  the 
early  experiments.  We  found  it  convenient  to  use  a  Bragg  cell  for  am¬ 
plitude  modulation  of  the  incident  light;  with  a  semiconductor  laser,  we 
could  directly  modulate  the  laser  itself.  Since  the  objective  lens  is 
not  specially  coated  for  the  He-Ne  wavelength,  there  is  still  some  light 
reflected  back  from  the  objective  lens.  To  reduce  the  amount  of  the  re¬ 
flected  light  reaching  the  detector,  we  simply  place  an  iris  diaphragm 
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Fig.  2.  The  dependence  of  the  detected  signal  amplitude  on 
object  distance  for  a  mirrored  surface. 

before  the  photodetector.  An  iris  diaphragm  placed  between  the  pinhole 
and  the  objective  lens  reduced  the  numerical  aperture  of  the  reflected 
light.  This  Increased  the  spot  size  of  the  returning  one-way  spatial 
filter.  The  object  was  moved  by  two  actuators  In  the  axial  direction 
z  and  the  transverse  direction  x  .  The  whole  system  was  controlled  by 
a  computer. 


RESULTS 

To  demonstrate  the  usefulness  of  this  system  for  obtaining  accurate 
distance  measurements,  we  used  a  mirror  as  the  reflecting  object  and 
scanned  it  in  the  axial  direction  z  .  Figure  2  shows  the  dependence  of 
the  signal  on  the  position  of  the  object.  The  numerical  aperture  of  the 
optical  system  is  0.054  ,  the  corresponding  depth  of  focus  is  188  pm  , 
and  the  spot  size  is  about  10  pm  .  The  experimental  curve  fits  the 
theoretical  curve  very  well,  except  for  some  discrepancy  in  the  side- 
lobes,  which  is  believed  to  be  caused  by  aberrations  in  the  optical  sys¬ 
tem.  By  determining  the  position  of  the  peak  of  the  curve,  we  can  de¬ 
termine  the  position  of  the  object  very  accurately. 
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Fig.  3.  One-dimensional  scan  over  a  tilted  mirror. 
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Fig.  One-dimensional  scan  over  a  220  pm  mirror  step. 


To  determine  the  sensitivity  of  our  measurement  technique,  we 
scanned  the  focused  beam  over  a  tilted  mirror,  as  shown  in  Fig.  3.  The 
straight  line  represents  the  mirrored  surface  and  the  dots  are  exper¬ 
imental  results.  The  estimated  accuracy  is  about  2  vim  .  As  we  tilted 
the  mirror,  the  height  of  the  central  lobe  decreased  and  the  depth  of 
focus  increased.  This  is  because  less  light  could  be  collected  by  the 
camera  lens,  so  the  effective  numerical  aperture  decreased. 

The  results  of  a  subsequent  scan  over  a  220  pm  step  on  a  specular 
reflector  are  shown  in  Fig.  4.  The  transverse  spacing  between  adjacent 
points  is  5  urn  .  Except  for  the  overshoot  near  the  step,  the  accuracy 
of  the  range  measurement  is  about  2  pm  .  The  overshoot  is  due  to  sha¬ 
dowing  and  Interference  effects  which  are  illustrated  in  Fig.  5.  The 
focused  light  is  reflected  by  both  top  and  bottom  surfaces.  Most  of  the 
light  reflected  from  the  bottom  surface  is  blocked  by  the  step,  except 
that  reflected  back  at  a  small  angle.  The  interference  between  the  two 
beams  shifts  the  center  peak  and  causes  the  overshoot  effect. 

This  system  is  not  only  suitable  for  measurements  on  smooth  sur¬ 
faces,  but  also  for  measurements  on  rough  surfaces.  This  is  because 
when  the  light  is  focused  to  a  point  on  a  rough  surface,  the  reflected 
light  will  be  scattered  in  all  directions,  but  an  image  of  this  point 
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Fig.  5.  Illustration  of  shadowing  and  interference  effects. 
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Fig.  6. 


Dependence  of  the  detected  signal  amplitude  on  the  object 
distance  for  a  raw  aluminum  sample. 


can  still  be  obtained  at  the  pinhole.  We  pay  a  price  for  this  conven¬ 
ient  result  -  the  light  reaching  the  pinhole  is  much  reduced  in  inten¬ 
sity;  the  advantage  is  that  the  alignment  of  the  rough  surface  is  rela¬ 
tively  uncritical.  To  demonstrate  this  effect,  we  scanned  a  sample  of 
sheet  aluminum  over  the  surface,  as  supplied  by  the  manufacturer.  The 
results  we  obtained  are  shown  in  Fig.  6.  The  central  lobe  of  this  curve 
is  unchanged  from  the  central  lobe  obtained  using  a  mirror.  As  expect¬ 
ed,  the  detected  signal  from  this  rough  surface  is  much  weaker  than  that 
for  a  smoother  surface,  but  it  is  still  strong  enough  to  be  easily  de¬ 
tected  with  a  narrowband  receiving  system.  The  results  obtained  were 
insensitive  to  the  tilt  of  the  object,  which  is  essential  for  the  mea¬ 
surement  of  ordinary  machined  parts. 


CONCLUSIONS 

We  have  demonstrated  an  optical  range  finder  with  a  range  accuracy 
of  2  pm  and  a  transverse  definition  of  10  urn  at  a  working  distance 
of  15  cm  .  The  device,  based  upon  the  type  II  microscope  principle,  Is 
very  stable,  easy  to  align,  and  suitable  for  measurement  over  both 
smooth  and  rough  surfaces. 
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Distance  measurements  by  differential  confocal  optical  ranging 


Timothy  R.  Corle,  Jeffrey  T.  Fanton,  and  Gordon  S.  Kino 


A  new  technique  is  described  for  measuring  the  distance  between  a  lens  and  reflecting  surface  extremely 
accurately.  It  is  based  on  the  sharply  peaked  depth  response  of  type  II  confocal  systems.  By  dithering  either 
the  sample  or  the  optical  system,  a  differential  measurement  is  generated,  placing  a  zero-crossing  at  the  peak 
of  the  depth  response  and  improving  the  ranging  accuracy.  The  technique  is  independent  of  surface 
roughness  or  tilt  and  hence  is  useful  for  robotics  or  machining  applications.  Sensitivities  to  surface  vibrations 
of  0.01  nm  and  thin  film  measurements  to  0.04  fim  demonstrated.  Signal-to-noise  calculations  are  presented, 
and  the  procedure  for  measuring  the  thickness  of  transparent  films  is  outlined. 


Type  II  confocal  scanning  optical  microscopes 
(CSOM)  differ  from  standard  (type  I)  confocal  micro¬ 
scopes  in  that  they  use  a  columnated  beam  to  illumi¬ 
nate  the  objective  lens,  thus  producing  a  small  spot  on 
the  sample.  The  wave  reflected  from  the  sample  is 
then  focused  to  a  point  detector. 12  The  net  effect  is  to 
square  the  lens  response,  thus  sharpening  the  point 
spread  function  compared  with  an  ordinary  micro¬ 
scope.  The  trade-off  comes  from  having  to  take  data 
one  point  at  a  time  and  build  up  the  image  using  a 
raster  scan.  CSOMs  are  also  known  for  their  shallow 
depth  of  field.3-4  By  exploiting  this  property  the  dis¬ 
tance  between  the  object  and  the  lens  can  be  accurately 
determined  by  looking  for  the  peak  in  the  amplitude 
response  V(z)  as  a  function  of  distance  z  from  a  reflec¬ 
tor  to  the  focal  point  of  a  lens.  We  demonstrate  that, 
by  vibrating  the  sample,  we  can  obtain  an  rf  output 
from  the  detector  at  the  vibration  frequency  fi,  which 
is  effectively  the  differential  of  |V(z)|2.  Thus  we  can 
obtain  a  zero  rf  output  when  I  V(z)|2  is  maximum.  This 
zero  point  can  be  used  to  measure  small  variations  in 
the  lens  object  separation  with  great  accuracy.  As 
with  direct  measurements  of  I  V(z)(2,  the  technique  is 
independent  of  surface  roughness  or  tilt.  We  also 
present  signal-to-noise  calculations  and  results  of 
measurements  for  thin  films  made  with  photoresist  on 
silicon. 


It  can  be  shown  that,  as  a  plane  reflector  is  scanned 
axially  through  the  focal  plane,  the  amplitude  response 
(unnormalized)  of  a  confocal  microscope  is  given  by 
the  relation4 


exp(2i7iz  cos 3)  P{B)R(0)d6. 


In  these  expressions  R(0)  is  the  reflection  coefficient  of 
the  plane  reflector,  P(0)  is  the  pupil  function  of  the 
lens,  k  is  the  wave  number  (2ir/X),  z  is  the  distance 
between  the  reflector  and  the  focal  plane,  and  sin0o  is 
the  numerical  aperture  of  the  lens.  Given  a  uniformly 
excited  lens  of  radius  a,  with  ka  »  1  and  (1  +  cos0)/ 
cos01/2  ~2,  an  approximate  expression  for  V(z)  has 
been  derived5: 


V(z)  =  exp[[fcz(l  +  cos0o)) 


sinfcz(l  -  cos#0) 


fcz(l  -  costf0) 


This  function  is  characterized  by  a  sharp  central  lobe, 
with  3-dB  depth  of  field  given  by5 


(1  -  cosfl0) 
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For  a  lens  with  numerical  aperture  of  0.9  in  air  and  X 
=  633  nm,  the  3-dB  width  can  be  as  small  as  500  nm, 
hence  distance  measurements  with  this  accuracy  are 
easily  achieved.  For  finer  resolution  it  is  desirable  to 
measure  the  position  of  the  maximum  of  the  reflectiv¬ 
ity  function  as  accurately  as  possible.  To  measure  this 
position  we  vibrated  the  sample  using  a  vibration  am¬ 
plitude  comparable  to  or  less  than  the  depth  of  focus. 
This  vibration  generated  a  differential  measurement 
as  shown  in  Fig.  1. 

Physically  the  effect  can  be  understood  by  referring 
to  Fig.  2.  When  the  object  is  in  the  focal  plane  of  the 
lens  (z  =  0),  the  return  rays  are  focused  at  the  pinhole 
and  a  maximum  in  the  output  signal  is  obtained.  If 
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Fig.  1.  Experimental  I  V(z)l2  and  V*(z)  V'(z)  curves.  These  results 
are  for  a  plane  reflector  and  a  len  of  numerical  aperture  0.8  in  air. 


the  object  is  then  vibrated  at  a  frequency  fi,  two  maxi¬ 
ma  will  be  obtained  in  each  cycle  corresponding  to  the 
two  passes  through  the  focal  plane.  Thus  an  output 
signal  is  obtained  at  the  second  harmonic  (211),  and  a 
zero  is  obtained  at  the  output  frequency  12.  When  the 
system  is  defocused,  however,  a  signal  at  frequency  12  is 
obtained  whose  magnitude  depends  on  both  the  aver¬ 
age  distance  from  the  reflector  to  the  focus  zq  and  the 
vibration  amplitude  Z\ 

If  geometry  or  size  considerations  make  it  inconve¬ 
nient  to  vibrate  the  sample,  the  lens  or  the  pinhole  may 
be  moved  to  achieve  the  same  effect.  Alternatively 
other  techniques  such  as  the  use  of  a  flexible  mirror 
whose  focal  length  can  be  changed  electronically,  or  an 
electrooptic  cell,  may  be  used  to  periodically  defocus 
the  image.6 

There  are  two  types  of  measurement  of  interest:  (1) 
the  determination  of  the  absolute  position  of  a  surface, 
or  by  looking  at  two  surfaces,  the  thickness  of  a  thin 
film;  and  (2)  the  measurement  of  vibration  ampli¬ 
tudes.  To  determine  absolute  distances  we  use  the 
zero-crossing  of  the  d|  V(z)|2/dz  curve,  while  in  the  lat¬ 
ter  case  it  is  apparent  that  we  should  choose  a  value  of 
z,  where  the  slope  of  I  V(z)l2  is  a  maximum. 

We  may  put  these  ideas  on  a  more  mathematical 
basis  as  follows.  The  small  sample  vibrations  can  be 
expressed  as  Zj  cosfit.  We  may  then  write 


z  =  z0  +  Zj  coaflt  =  z0  +  Az.  (4) 

Expanding  I  V(z)|2  in  a  Taylor  Series  about  Zo  we  ob¬ 
tain,  to  the  third  order  in  z, 


IV(z)l2  -  |  V(z0)|2  +  AZ 


alv(z)P, 


dz 

2  a2l  v(z>!2 1 


+  . . . 


dz 

I  V(z0)|2  +  Z[  ccraQf  2Re[y*(z0)V'(z0)| 


+  z?  cos2«f  Re[V"*(z0)V"(z0)  +  V*(z0)V'(z0)l 
z3 

+  y  cosJ«f  Re(3V"*(z0) V"(z0)+  V*(z0 V“(z0)]  +. . .. 

(5) 


It  will  be  observed  from  Eq.  (5)  that  the  leading 


Fig.  2.  Illustration  of  the  effects  of  defocus  on  a  type  II  microscope. 


contribution  to  the  rf  term  at  12  is  of  the  form 
V*  (z0)  V'(zo).  This  derivative  term  is  zero  at  zo  =  0. 

It  follows  from  Eq.  (2)  that  we  can  also  expand  V(z0) 
in  a  series  of  the  form 

K^z2  KAzA 

V(z0)  =  exp(ifcz0)  [1 - -±  +  —±  +  ...],  (6) 

where  K  =  k(l  —  cos 80).  Substituting  this  expansion 
and  its  derivatives  into  Eq.  (5),  to  the  fourth  order  in  z0 
and  z\  we  find  that 


3  45  6  15 


+  cosfit 


/  2K2z0z,  8ff4ZpZ,  2/f,z0z3\ 

\  3  +  45  +  15  /  ’ 


(7) 


Equation  (7)  gives  us  a  series  approximation  for  the 
D.C.  and  first  harmonic  (cosfit)  terms  in  Eq.  (5).  It 
may  be  used  to  calculate  the  optimum  values  of  zo  and 
z\  for  the  measurement  of  absolute  distances  and  vi¬ 
bration  amplitudes.  To  measure  absolute  distance 
accurately  we  work  near  the  focal  point,  where  zo  is 
small.  Keeping  liner  terms  in  z0  we  find  that  the  rf 
output  varies  with  the  vibration  amplitude  and  is  max¬ 
imum  when  K2z  x2  =  5/3,  or 

0-20* 

z'  (1  -  cos 0o) 


a  movement  of  only  0.35X  for  a  0.9  numerical  aperture 
lens.  With  this  vibration  amplitude  Eq  (7)  reduces  to 

I  t/(z)|2  =  0.72  +  cosSif  ( —0.57  Kz0).  (9) 


The  output  signal  at  frequency  fi  becomes  zero  when 
zo  =  0,  i.e.,  when  I  V(z0)l2  is  maximum.  The  sign  of  the 
rf  signal  is  determined  by  the  term  in  the  brackets,  so 
the  rf  signal  is  positive  for  negative  zo  and  negative  for 
positive  zo  as  in  Fig.  1.  Hence  a  true  zero-crossing  is 
obtained  at  z0  =  0. 
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The  best  sensitivity  for  the  detection  of  acoustic 
waves  and  other  small  rf  vibrations  is  obtained  when 
we  operate  at  the  maximum  slope  of  the  I  V(z) |2  curve. 
The  maximum  slope  is  obtained  when  K?zo2  =  15/8,  or 


0.22A 

0  (1  —  cosfl0) 


(10) 


If  we  include  only  the  zeroth  and  first  harmonic  terms, 
the  magnitude  of  I  V(z) I2  when  z0  is  fixed  at  this  value  is 

I  V(z)|2  =  0.63  +  cosfif(-0.49ifz,).  (H) 

We  are  now  in  a  position  to  determine  the  SNR  of 
this  technique.  We  shall  carry  out  the  derivation  for 
small  displacements  and  determine  from  Eq.  (9)  the 
minimum  value  of  zq  that  can  be  detected. 

The  power  entering  the  detector  is  given  by 

P  =  P0\V(z)\2,  (12) 

where  Po  is  the  power  reaching  the  detector  at  zo  =  0. 
The  current  I  from  the  detector  is 


qyP  _  qvP0\VU)\2 
hv  hv 


(13) 


where  q  is  the  charge  of  an  electron,  q  is  the  quantum 
efficiency  of  the  detector,  h  is  Planck’s  constant,  and  v 
is  the  optical  frequency. 

The  mean-square  noise  current  from  the  detector 
(in2)  is 

<0  =  2  qloB  +  4kTBF/R0,  (14) 

where  Iq  is  the  average  current  from  the  detector,  k  is 
Boltzmann’s  constant,  T  is  the  temperature,  B  is  the 
bandwidth,  /?o  is  the  input  impedance,  and  F  is  the 
noise  figure  of  the  preamplifier  when  supplied  from  a 
high  impedence  detector. 

We  shall  assume  that  either  the  input  impedence  R0, 
the  input  power,  or  both  quantities,  are  large  enough  so 
that  the  last  term  in  Eq.  (14)  can  be  ignored  and  the 
system  is  shot-noise  limited.  In  this  case  we  can  write 

,  2  q2yPB 

Hi)  =  2 qloB  =  ■  (15) 


-  It  follows  from  Eqs.  (9)  and  (13)  that  the  mean- 
square  signal  current  from  the  detector  is 


, 2  _  1  0.57  Kz0qyP0l? 
'  •'  2 1_0.72  hv  J 

So  the  SNR  is 


(16) 


SNR  = 


<*n> 


=  0.11 


Pol(Kz0)2 

hvB 


(17) 


For  a  signal-to-noise  of  unity,  the  minimum  detectable 
value  of  zq  is  z0min: 


=  0.40A  \hvB-\'n 

°"i"  (l-cos@0)  P0 17 


(18) 


As  an  example  with  X  =  633  nm,  hv  =  3.13  X  10-19  J, 
and  using  a  lens  with  an  angle  aperture  of  0.9  in  air,  a 
detector  with  quantum  efficiency  q  =  0.8,  a  bandwidth 
of  1  kHz,  and  Po  =  0.1  mW,  we  find  zomi„  =  8.9  X  10-4 
nm.  A  similar  calculation  based  on  Eq.  (11),  using  the 


i 

i 


Fig.  3.  Schematic  of  the  electronically  controlled  type  II  micro¬ 
scope  used  to  make  the  measurements. 


nm. 


same  parameters,  indicates  that  the  minimum  detect¬ 
able  value  of  zx  is,  zlmin  =  8.6  X  10-4  nm.  Hence,  with 
our  technique  in  a  shot-noise  limited  system,  intensity 
measurements  can  yield  interferometric  accuracies. 

To  measure  I  V(z)|2  and  V*(z)V'(z)  we  used  the  elec¬ 
tronically  controlled  type  II  microscope  shown  in  Fig. 
3.  The  sensitivity  of  our  technique  to  small  surface 
movements  was  determined  by  attaching  a  reflecting 
surface  to  a  PZT  bulk  wave  transducer.  Scans  were 
made  while  the  transducer  was  driven  at  150  kHz. 
Figure  4  shows  the  results  obtained  for  oscillations 
with  an  amplitude  of  Z\  =  0.01  nm.  As  can  be  seen,  the 
signal  is  still  quite  strong,  indicating  surface  move¬ 
ments  in  this  range  can  easily  be  detected.  The  results 
in  Fig.  4  also  imply  that  our  experimental  apparatus  is 
not  shot-noise  limited,  primarily  because  the  input 
impedance  of  our  amplifier  was  only  50  ft. 

In  addition  to  excellent  distance  resolution,  this 
technique  has  several  other  characteristics  which 
make  it  useful  for  remote  distance  sensing.  First,  the 
system  has  excellent  transverse  resolution.  Because 
the  objective  lens  is  used  twice,  once  to  focus  the  in¬ 
coming  plane  waves  to  a  point  and  a  second  time  to 
image  this  point  onto  the  pinhole,  the  transverse  reso¬ 
lution  is  given  by  the  square  of  the  lens  response.  For 
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Fig.  5.  |V(z)l2  and  V’(z)V'(z)  curves  for  1.8  gm  of  photoresist  on 
silicon. 


Fig.  6.  I  V(z)|2  and  V*(z)V'(z)  curves  for  4.86  of  photoresist  on 
silicon. 

a  0.9  numerical  aperture  lens  and  X  =  633  nm,  the  3-dB 
width  of  the  point  spread  function  is  only  0.35  pm,  a 
27%  improvement  over  conventional  microscopes. 
Also,  it  has  been  demonstrated  that  the  shape  of  the 
central  lobe  of  I  V(z)\ 2  is  independent  of  surface  rough¬ 
ness  or  tilt  (although  the  amplitude  does  vary  with  the 
reflectivity  of  the  sample).4  Hence  V*(z)V'{z)  mea¬ 
surements  are  also  independent  of  surface  conditions. 
Finally,  since  this  technique  only  depends  on  the  fo¬ 
cusing  properties  of  the  objective  lens,  phase  stabiliza¬ 
tion  of  the  system  is  not  necessary.  A  multimode 
semiconductor  laser,  or  even  an  incoherent  LED 
source,  can  be  used. 

We  can  apply  differential  techniques  to  the  mea¬ 
surement  of  thin  Films.  With  a  thin  film  more  than 
one  reflecting  plane  is  present,  and  the  |V(z)|2  curve 
has  two  peaks  corresponding  to  the  focusing  on  at  the 
top  and  bottom  surfaces  of  film.  In  this  case  differen¬ 
tial  techniques  may  be  used  to  improve  the  accuracy  of 
the  Film  thickness  measurements.  Results  are  shown 
in  Fig.  5  for  1.8  gm,  and  in  Fig.  6  for  4.86-pm  films.  For 
thick  films  such  as  the  4.86-pm  Film,  the  results  ob¬ 
tained  are  very  accurate.  For  thinner  Films,  where  the 
two  peaks  tend  to  overlap,  the  measurement  of  film 
thickness  is  accurate  only  of  the  order  of  half  a  wave¬ 
length.  This  inaccuracy  is  primarily  due  to  interfer¬ 
ence  with  an  unknown  phase  difference  between  the 
reflected  waves  from  the  top  and  bottom  surface  of  the 


Film.  This  problem  can  be  overcome  by  using  a  broad¬ 
band  source  with  a  coherence  length  less  than  the  film 
thickness,  or  by  direct  comparison  with  calibrated 
samples.  Alternatively  the  equations  for  V*(z)V'(z) 
could  be  fit  to  the  data  in  an  iterative  manner,  but  this 
approach  is  often  difficult  and  time-consuming.  For 
thicker  films,  where  the  peaks  do  not  overlap,  interfer¬ 
ence  does  not  occur,  so  the  film  thickness  may  be 
directly  measured. 

To  determine  the  film  thickness  from  the  data  one 
needs  to  account  for  the  index  of  refraction  of  the  film 
material.  One  approach  is  to  modify  the  form  of  Eq. 
(1).  Let  the  reflection  coefficients  from  the  top  and 
bottom  surfaces  be  Ri  and  R?.  We  shall  denote  all 
quantities  outside  and  inside  the  film  by  the  subscripts 
1  and  2,  respectively,  and  for  simplicity  we  shall  take  R  \ 
and  R<i  to  be  independent  of  0. 

In  this  case  Eq.  (1)  can  be  rewritten  as 


i  (1  —  cosfl^sin#! 
COS0, 


exp(2i>,z  cosSj) 


X[R,  +  R2  exp(2ifc2/i  cos02)]  P(0i)dU{,  (20) 


where  h  is  the  film  thickness,  and  k\  and  are  the 
wave  numbers  outside  and  inside  the  film.  The  angle 
02  can  be  found  by  using  Snell’s  law, 


n\ 

sin02  =  —  sinO, 
n2 


(21) 


where  n i  is  the  index  of  refraction  of  the  material 
surrounding  the  objective  lens  (in  our  case,  air),  and  ri2 
is  the  index  of  refraction  of  the  film  material. 

Direct  integration  of  Eq.  (20)  for  h  =  4.86  pm,  n  = 
1.63,  P(0)  =  1,  Ri  =  0.39,  and  i?2  =  1.0  reproduces  the 
data  fairly  accurately  as  shown  in  Fig.  7.  Past  experi¬ 
ence  has  shown  us  that  we  can  improve  upon  this  result 
by  introducing  a  pupil  function  of  the  form 


P(0)  =  exp(A  sin20)  exp(i'B  sin40).  (22) 


Values  of  A  =  0.2  and  B  =  6  have  been  shown 
previously  to  account  for  aberrations  in  our  system.4 
In  this  case,  however,  the  introduction  of  phase  aberra¬ 
tions  makes  things  much  worse  as  shown  in  Fig.  7. 
This  behavior  could  be  the  result  of  unknown  phase 
aberrations  within  the  film,  the  fact  that  the  aberra¬ 
tions  of  Eq.  (22)  are  an  oversimplification,  or  other 
system  aberrations  not  yet  understood. 

From  the  integration  of  Eq.  (20)  we  measure  the 
distance  between  the  two  maxima  (hc).  The  ratio  of 
the  actual  film  thickness  to  the  calculated  distance  (h/ 
hc)  gives  us  the  multiplicative  factor  needed  to  account 
for  the  index  of  refraction  of  the  film.  This  multiplica¬ 
tive  factor  is,  to  the  first  order,  independent  of  the  film 
thickness.  Thus  the  ratio  h/hc  need  only  be  calculated 
once  for  each  index  material.  For  a  photoresist  film 
with  n  =  1.63  at  X  =  633  nm,  h/hc  =  2.08  yielding  a 
measured  film  thickness  of  4.86  pm.  Lack  of  detailed 
understanding  of  the  system  aberrations  remains  the 
limiting  factor  in  improving  the  accuracy  of  the  thin 
film  measurements.  With  the  current  level  of  under¬ 
standing  we  feel  that  accuracies  of  0.04  pm  (1%)  are 
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Fig.  7.  Comparison  of  experimental  data  and  theoretical  calcula¬ 
tions  for  a  4.86-jim  film  of  photoresist  on  Silicon. 


where  n  is  the  index  of  refraction  of  the  film  material, 
and  the  material  outside  the  film  is  assumed  to  be  air 
( n  -  1).  To  determine  the  focal  shift  of  the  beam,  and 
hence  the  multiplicative  factor,  we  integrate  Eq.  (23) 
over  all  angles  in  the  focused  beam: 

h  Jo  cos0 


For  photoresist  with  n  =  1.63  at  X  =  630  nm,  h/hc  = 
2.12  giving  a  thickness  of  4.91  pm.  While  this  method 
is  not  strictly  correct  since  it  does  not  account  for 
diffraction  within  the  film,  it  offers  a  quick  way  of 
calculating  the  correction  factor  to  within  a  couple  of 
percent. 

We  have  demonstrated  that,  by  vibrating  the  sample 
in  a  confocal  scanning  microscope,  the  distance  be¬ 
tween  the  lens  and  the  sample  can  be  accurately  deter¬ 


mined.  The  technique  is  based  on  the  shallow  depth 
of  focus  properties  of  confocal  optical  microscopes.  It 
exhibits  excellent  transverse  resolution  and  can  be 
used  to  image  rough  samples.  It  also  has  the  advan¬ 
tage  of  being  phase  independent,  hence  cheap  diode 
lasers  or  incoherent  sources  may  be  used.  Signal-to- 
noise  calculations  are  presented  which  illustrate  that 
our  intensity  technique  is  comparable  in  sensitivity  to 
phase-dependent  interferometric  techniques.  The 
thicknesses  of  photoresist  films  have  been  measured  to 
within  0.04  pm  and  theoretical  calculations  presented 
which  correct  for  the  index  of  refraction. 


Fig.  8.  Refraction  of  a  ray  incident  upon  a  dielectric  resulting  in  a 
change  in  the  focal  plane  for  a  focused  beam. 


reliably  obtained.  Ellipsometry  readings  taken  of  the 
same  sample  give  a  thickness  of  4.86  pm  at  index  1.64. 

A  somewhat  simpler  method  of  computing  h/hc  is  by 
calculating  the  distance  the  focal  point  is  refracted  into 
the  material  by  the  presence  of  the  dielectric.  A  given 
ray  impinging  at  an  angle  0  on  a  dielectric  interface  will 
cross  the  optical  axis  of  the  system  at  a  distance  (/) 
from  the  surface  of  the  film  (Fig.  8).  The  ratio  be¬ 
tween  the  distance  (/)  and  the  point  where  the  ray 
would  intersect  the  optical  axis  if  no  film  was  present 
(do)  ia  given  by 

,  {23) 

d0  cos0 
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Abstract 


We  describe  two  ways  which  we  have  developed  for  obtaining  phase  information  in 
Type  2  (confocal)  microscopes:  heterodyne  interferometry  using  a  Bragg  cell  for  both 
scanning  and  frequency  shifting,  and  a  novel  "heterodyne  phase  contrast"  technique  using  a 
t i me- va ri ab 1 e  el ect ro -opt i c  phase  plate,  somewhat  similar  to  a  Zernike  phase  plate  with 
mechanical  scanning.  The  heterodyne  interferometer  can  measure  amplitude  to  0.2%  and 
phase  to  the  0.1°  level  at  a  rate  of  30,000  to  50,000  points  per  second,  with  exel- 
lent  vibration  immunity.  The  heterodyne  phase  contrast  system  is  slower,  but  is  a  single¬ 
beam  system  and  could  be  retrofitted  to  a  standard  Type  2  system  without  affecting  normal 
operation.  If  scan  stage  vibration  can  be  neglected,  for  example  in  acoustic  wave  detec¬ 
tion,  it  should  be  able  to  measure  opti cal  phase  to  around  0.1  millidegree;  however,  this 
accuracy  is  not  yet  demonstrated. 


I nt  roduct i on 


The  purpose  of  phase  measurements  in  microscopy  is  well  known:  they  provide  addi¬ 
tional  information  about  the  sample,  independent  of  the  intensity  information,  which  may 
be  obtainable  in  no  other  way.  Conventional  phase  contrast,  as  well  as  shearing  and  No- 
marski  interference  microscopes,  perform  very  well  for  optically  thin  or  highly  reflective 
samples;  the  lack  of  amplitude  contrast  with  such  samples  means  that  the  image  is  simply  a 
phase  map.  These  instruments  are  incapable  of  producing  quantitative  phase  data  for  more 
general  classes  of  samples  because  they  lack  any  provision  to  separate  the  amplitude  and 
phase  contributions  to  the  contrast.  Our  systems  remove  this  ambiguity  by  using  AC  tech¬ 
niques  to  sample  the  optical  amplitude  and  phase  simultaneously  and  independently. 

Heterodyne  I nt e rf e romet ry 

The  idea  of  this  system  has  been  presented  before,1,2  and  so  will  be  described  only 
briefly.  Figure  1,  re printed  from  reference  1,  shows  a  simplified  version  of  the  optical 
system  of  the  heterodyne  interference  microscope.  A  single-frequency  laser  beam  at 
514.5  nm  comes  into  an  acousto-optic  Bragg  cell  which  is  placed  at  the  pupil  of  a 
factory-modified  Leitz  Orthoplan  microscope.  An  undiffracted  beam  passes  through  the 
Bragg  cell,  along  with  a  steerable  first-order  beam  which  has  been  downshifted  by  the 
acoustic  frequency  f b  .  These  beams  are  focused  to  d i f f ra c t i on -1 i mi t ed  spots  on  the 
sample.  The  undiffracted  beam  forms  a  stationary  reference  spot,  the  steerable  beam  a 
spot  which  scans  as  the  acoustic  frequency  is  varied. 


BRAGG  CEIL 
60-110  MHz 

I 


STATIONARY 

BEAM 


MICROSCOPE  f° 


Fig.  1.  Simplified  diagram  of  the  heterodyne  interference  microscope. 

The  beams  reflected  from  the  sample  retrace  their  paths  and  pass  through  the  cell 
again,  in  the  other  direction.  Of  the  resulting  four  beams,  we  select  the  two  (coaxial) 
singly-diffracted  ones,  one  from  each  spot.  The  beam  from  the  scanned  spot  has  been 
downshifted;  the  one  from  the  stationary  spot  has  been  upshifted.  These  two  beams  inter¬ 
fere  on  a  photodiode,  producing  a  product  signal  at  twice  the  acoustic  frequency.  The 
normalized  complex  amplitude  of  this  signal,  an  RF  tone  in  the  120-220  MHz  band,  is 


92  /  SPIE  Vo!  749  Metrology  Figure  and  Finish  11987) 


equal  to  the  complex  reflection  coefficient  of  the  sample  at  the  position  of  the  moving 
spot.  Now  the  problem  of  measuring  the  optical  amplitude  and  phase  of  the  scanned  beam  is 
reduced  to  an  RF  measurement,  which  can  be  done  very  quickly  and  accu rat e 1 y  .  > 3  We  have 
demonstrated  scan  speeds  of  35,000  points  per  second  with  accuracies  of  a  few  tenths  of 
a  degree  in  phase,  and  0.2%  in  amplitude  with  a  dynamic  range  of  70  dB  .  This  level  of 
phase  error  corresponds  to  a  height  resolution  of  a  few  Angstroms  . 

The  signal  from  the  photodetector  is  the  product  of  the  amplitude  illumination  pat¬ 
tern  (the  scanned  spot),  the  sample  reflection  coefficient  as  a  function  of  position,  and 
the  sensitivity  pattern  (the  stationary  spot  superimposed  by  the  Bragg  cell).  Since  the 
illumination  and  sensitivity  patterns  are  identical  jinc  functions,  the  microscope's  point 
spread  function  is  their  square;  thus,  the  microscope  is  of  Type  2.  We  refer  to  this  pro- 
perty  as  a  “virtual  pinhole." 

Since  the  reference  beam  is  on  the  sample,  the  system  has  excellent  immunity  to 
sample  vibrations;  any  vertical  vibration  of  the  sample  causes  equal  phase  shifts  in  the 
two  beams,  which  then  cancel  each  other  when  the  two  beams  interfere.  Horizontal  vibra¬ 
tions  are  no  more  of  a  problem  than  in  conventional  microscopes;  twisting  motions  can 
introduce  differential  phase  shifts,  but  as  the  two  spots  are  usually  100  urn  apart  (or 
less),  this  is  not  very  serious. 

This  system  has  one  disadvantage,  as  compared  to  single-beam  microscopes,  which  it 
shares  with  shearing  interference  microscopes:  it  demands  a  smooth  area  for  the  reference 
spot.  If  this  condition  is  not  met,  the  sensitivity  pattern  will  no  longer  be  identical 
to  the  illumination  pattern,  and  the  microscope  will  no  longer  be  of  Type  2.  This  condi¬ 
tion  is  usually  accompanied  by  a  large  decrease  in  signal  level.  The  smooth  area  need  be 
only  a  micron  or  so  in  diameter,  so  this  isn't  a  particularly  serious  fault. 

Since  both  amplitude  and  phase  are  available,  we  exploit  the  extra  spatial  frequency 
bandwidth  of  the  Type  2  transfer  function4  with  an  inverse  filter  to  achieve  a  factor  of 
two  sharper  resolution^  at  any  given  wavelength.  The  filtering  is  done  digitally  and 
takes  less  than  <C  seconds  for  1024  points  on  a  personal  computer. 

Heterodyne  Phase  Contrast 

The  heterodyne  phase  contrast  system  is  shown  in  Fig.  2;  it  is  a  standard  Type  2 
microscope  with  an  electro-optic  phase  plate  at  the  pupil.  The  E-0  plate,  shown  in 
Fig.  3,  is  a  2  mm  thick  block  of  PLZT  with  a  transparent  indium  tin  oxide  electrode 
covering  a  central  disk.  Since  PLZT  is  an  el ect rost ri ct i ve  material  with  large  off- 
diagonal  coupling  coefficients,  an  applied  voltage  causes  the  optical  path  length  through 
the  central  region  to  increase  or  decrease;  thus  the  relative  phase  of  the  pupil's  central 
region  and  outer  annulus  can  be  modulated.  There  is  an  obvious  similarity  between  our  E-0 
plate  and  the  Zernike  phase  plate;  however  the  principle  is  quite  different,  as  will 
presently  be  seen.  A  complete  analysis  of  the  imaging  behavior  of  this  system  in  both  the 
horizontal  and  the  vertical  directions  is  under  way,  to  be  published  at  a  later  date,  but 
the  main  features  of  the  phase  imaging  may  be  shown  by  the  following  heuristic 
presentation. 


ELECTRO-OPTIC 
PHASE  SHIFTER 


MICROSCOPE 


Fig.  2.  Schematic  of  the  mechanically-scanned  heterodyne  phase  contrast  microscope 
showing  the  placement  of  the  el ect ro -opt i c  phase  plate. 


Consider  an  optical  beam  focused  just  inside  a  dielectric  film  (nj),  which  covers  a 
reflecting  plane  (Fig.  4).  The  top  surface  is  considered  to  be  z  =  0  .  If,  in  the  ab¬ 
sence  of  the  film,  the  beam  would  be  focused  at  a  point  z  ,  then,  because  of  refraction, 
the  outermost  rays  will  actually  cross  the  axis  at  a  point  z'  ,  where: 
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Fig.  3.  Illustration  of  the  electro-optic  plate  used  to  phase  shift  a  portion  of  the 
i 1 1 umi nati ng  beam. 


(b) 


Fig.  4.  (a)  An  optical  beam  focused  into  a  dielectric  ni  illustrating  the  path  dif¬ 

ference  between  the  center  and  outermost  rays.  (b)  The  same  beam  focused  into  a 
dielectric  ng  . 


z(nf 


sin2  n)1'2 


COS  8 


and  8  is  the  angle  the  rays  make  with  the  axis.  The  path  difference  between  the  axial 
and  outermost  rays  is: 


AZ  ' 

X 

/z'2  +  r2  -  z ' 

=  [z  2  ( n  2 

-  sin2 

0  )  / c  o  s  2  8  +  z2  tan2  e]1^2 

-  z  { n  2 

-  sin2 

e)1/2/cos  9 

so  that  the  phase  difference  is: 
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Hence,  the  phase  difference  between  the  rays  at  any  given  plane  depends  on  the  film's 
index  of  refraction.  If  we  use  the  e 1 ect ro -op t i c  plate  to  modulate  the  phase  of  the  axial 
ray,  that  phase  will  be  exp[i(ut  +  a  cos  at)]  ;  the  phase  of  the  outermost  ray  is 
exp[i(u>t  +  A  4> )  1  .  The  two  rays  will  interfere  on  the  detector  (thanks  to  the  pinhole), 
and  the  resulting  photocurrent  will  have  the  form: 


Ip  =  cos(a  cos  fit  -  &4>)  +  (other  DC  terms) 

For  small  modulation  index  (a  <<  1),  the  AC  term  is  I_  =  a  sin  a  ;>  cos  at  .  We  see  that  a 
signal  at  the  modulation  frequency  is  obtained  whose  amplitude  gives  the  sine  of  the  op¬ 
tical  phase;  if  it  is  multiplied  by  the  modulating  signal  in  a  balanced  mixer  or  synchron¬ 
ous  rectifier,  the  resulting  bipolar  output  voltage  will  faithfully  reproduce  this 
characteristic.  The  DC  term  is  just  the  standard  Type  2  microscope  response,  which  is 
essentially  unaltered  by  the  phase  plate's  action.  The  AC  signal  still  contains  the 
amplitude  informa'  on  as  well  as  the  phase.  To  obtain  truly  independent  data,  we  must 
divide  the  AC  voltage  by  the  DC  term,  an  elementary  task  whether  done  in  analog  circuitry 
or  in  so f twa re . 


Fig.  5.  Amplitude  point  spread  function  of  the  heterodyne  phase  contrast  system  as 
compared  with  Type  1,  Type  2,  and  coherent  Zernike  microscopes. 

We  have  calculated  the  in-focus  point  spread  function  ( P  S  F )  of  the  microscope  from 
Fouri'er  optics  for  both  the  AC  and  DC  components.  The  DC  PSF  is  only  very  slightly  af¬ 
fected  by  the  presence  of  the  E-0  plate;  it  is  equivalent  to  the  conventional  Type  2  PSF 
shown  in  Fig.  5.  The  AC  PSF  has  a  slightly  narrower  main  lobe  and  somewhat  higher  side- 
lobes  than  the  DC  response;  the  difference,  however,  is  by  no  means  as  pronounced  as  that 
between  the  Type  1  and  coherent  Zernike  systems.  In  dramatic  contrast  to  the  coherent 
Zernike  system,  the  AC  PSF  of  our  microscope  exhibts  a  very  compact  structure  with  tightly 
confined  sidelobes;  there  is  none  of  the  halo  or  shading  off  which  bedevils  edge  measure¬ 
ments  with  Zernike  systems. 

Although  we  are  still  in  the  early  stages  of  this  work,  by  analogy  with  other  sys¬ 
tems  the  ultimate  height  resolution  of  the  system  should  be  shot  noise  limited  at  around 
10  A  .  As  it  is  a  one-beam  system,  however,  vertical  vibrations  will  limit  the  resol¬ 
ution.  In  practice  the  objective  to  sample  distance  is  generally  only  a  few  millimeters, 
and  so  vibrations  can  be  minimized  by  careful  mechanical  design.  Some  applications,  such 
as  detection  of  acoustic  waves,  are  insensitive  to  these  low-frequency  vibrations  and  so 
could  take  advantage  of  the  instrument's  full  capabilities. 

In  this  early  ersion,  the  inner  region  of  the  E-0  cell  was  driven  with  a  10  kHz 
sine  wave  of  300  V  (p-p)  amplitude,  producino  a  phase  shift  of  about  6°  in  the  inner 
region  of  the  pupil  (later  versions  will  drive  the  outer  annulus  in  order  to  eliminate  the 
need  for  a  conducting  path  to  the  center  electrode).  Figure  6  shows  the  bipolar  AC  res¬ 
ponse  observed  while  a  plane  reflector  was  scanned  through  the  focal  region,  the  curve 


was  obtained  by  measuring  the  AC  output  at  10  kHz  with  a  single-channel  lock-in  ampli¬ 
fier.  We  see  that,  as  the  plane  reflector  moves  away  from  the  focal  plane,  a  signal  is 
generated  proportional  to  the  distance  off  focus,  as  the  simple  theory  predicts. 


Fiq.  6.  Variation  of  the  AC  signal  with  depth  in  the  heterodyne  phase  contrast  system. 
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